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AUTOM/'ted  calculation  of  protection  factors 
FOR  the  sodium  CHLORIDE  RESPIRATOR 
QUANTITATIVE  FIT  TEST  INSTRUMENT 


INTRODUCTION 

The  purpose  of  this  report  is  to  present  an  automated  procedure  for  cal¬ 
culating  a  respirator's  protection  factor  (PF)  afforded  to  the  respiratory 
tract  and  eyes  against  chemical  warfare  (CW)  agents  in  particulate,  aerosol, 
or  vapor  form.  The  participating  nations  of  the  Air  Standardization  Coordi¬ 
nating  Conwittee  (ASCC)  have  drafted  a  North  Atlantic  Treaty  Organization 
(NATO)  Standardization  Agreement  (STANAG)  recommending  a  sodium  chloride  quan 
titative  fit  test  scheme  for  this  procedure  [1].  The  Chemical  Defense  Estab¬ 
lishment  (CDE)  in  the  Ministry  of  Defense  (MOO)  of  the  United  Kingdom  is 
credited  with  developing  the  basic  sodium  chloride  (NaCI)  quantitative  fit 
test  technology  [2-11].  The  United  States  Air  Force  School  of  Aerospace  Medi 
cine  (USAFSAM),  which  has  had  several  years  of  laboratory  experience  with  the 
sodium  chloride  quantitative  fit  test  instrumentation  [1],  has  reviewed  and 
endorsed  this  draft  agreement. 

Our  experience  with  the  sodium  chloride  quantitative  fit  test  instrumen¬ 
tation  has  shown  one  area  to  be  of  particular  concern;  namely,  the  method  of 
reducing  the  collected  data  (respirator  sodiun  chloride  penetration  concentra 
tion  for  a  particular  exercise  ijrotocol)  and  calculating  a  protection  factor. 
Tht.  STANag  agreement  draft  proposes  that: 

"'ILen  the  output  of  the  flame  photometer  fluctuates  during 
o  measunment  of  penetration,  the  maximum  output  1$  to  be 
used  to  calculate  the  Protection  Factor.  For  this  pur¬ 
pose,  occasional  transient  increases  in  output  the  dura¬ 
tion  of  each  of  which  does  not  exceed  2  sec  may  be 
ignored."  [1] 

In  addition,  various  interested  organizations  have  developed  and  reported  on 
the  following  data  reduction  schemes  for  respirator  quantitative  fit  testing 
[2,13-43]: 

a.  selection  of  the  overall  maximum  output  peak 

b.  arithmetic  average  of  the  maximum  output  peaks 

c.  arithmetic  average  of  the  maximum  output  peaks  and  minimum 
valleys  (midpoint) 

d.  visual  estimation  of  the  midpoint  between  the  maximum  output 
peaks  and  minimum  valleys 

e.  time-averaged  or  integrated  value. 


This  report  develops,  with  some  rigor,  an  automated  procedure  to  reduce 
the  NaCl  leak  test  data  and  calculate  a  protection  factor.  An  overview  of  the 
USAFSAM  sodium  chloride  respirator  quantitative  fit  test  instrument  is  fol¬ 
lowed  by:  a  discussion  on  conventional  protection  factor  calculations;  the 
use  of  a  voltage-to- frequency  (V/F)  converter  circuit  to  do  time-averaged 
integration;  and,  finally,  a  least-squares  curve  fit  computer  program  to  cal¬ 
culate  a  protection  factor. 


USAFSAM  SODIUM  CHLORIDE  RESPIRATOR 
QUANTITATIVE  FIT  TEST  INSTRUMENT 

The  COE  sodium  chloride  respirator  quantitative  fit  test  method,  de¬ 
scribed  in  British  Standards  4400  and  2091,  has  been  adapted  and  modified  by 
USAFSAM  to  measure  the  protection  factor  in  the  respiratory  and  eye  compart¬ 
ments  of  aircrew  chemical  defense  respirators  This  instrument  gener¬ 

ates  a  solid  aerosol  of  sodium  chloride  crystals  as  the  challenge  atmosphere. 
The  concentration  of  the  challenge  atmosphere  in  an  aircrew  respirator  is 
measured  using  a  hydrogen  flame  photometer,  and  the  result  is  displayed  on  a 
strip-chart  recorder.  This  technique  allows  protection  factors  as  high  as 
10°  to  be  calculated. 


Instrument  Description 

Illustrated  in  Figure  1  are  the  primary  components  used  in  the  USAFSAM 
sodiun  chloride  respirator  quantitative  fit  test  instrument  [1,3,11,33].  The 
sodium  chloride  solid  aerosol  challenge  atmosphere  is  generated  by  atomizing  a 
sooiint  chloride  solution,  drying  the  liquid  aerosol  in  a  drying  tube,  and 
delivering  the  dry  cloud  to  the  top  of  a  transparent  plastic  hood.  A  subject, 
having  donned  a  respirator,  enters  the  hood  and  performs  a  series  of  breathing 
and  head  movement  exercises.  A  sampling  pump  is  used  to  draw  a  portion  of  the 
atmosphere  from  the  interior  compartment  of  the  respirator.  This  sample  is 
vaporized  in  a  hydrogen  flame  photometer,  and  the  output  signal  is  displayed 
on  a  strip-chart  recorder. 

Production  of  the  Sodium  Chloride  Challenge  Atmosphere 

The  challenge  atmosphere  is  generated  by  atomizing  a  5*  aqueous  sodium 
chloride  solution,  prepared  by  dissolving  50.00  grams  of  sodium  chloride  in 
lOOQ  ml  of  distilled  water.  Operating  a  Oautrebande  atomizer  at  70  psi 
(482.3  kPa)  (clean  dry  compressed  air)  yields  a  flow  through  the  atomizer  of 
35  liters/m^n  (STP).  The  liquid  aerosol  from  the  atomizer  is  then  injected 
perpendicularly  into  the  air  stream  ficswing  through  a  mixing  and  drying  tube. 
The  50  liter/min  (STP)  source  of  clean  dry  air  evaporates  water  from  the 
liquid  sodium  chloride  aerosol  and  produces  a  solid  sodium  chloride  aerosol 
challenge  atmosphere.  The  mass  median  aerodynamic  diameter  (NNAD)  of  the  dry 
sodtunj  chloride  crystals  range  from  0.4  -  0.6  wm.  The  sodium  chloride  chal¬ 
lenge  atmosphere  is  delivered  through  a  short  length  of  Tygon  tubing  to  the 
top  center  of  a  transparent  plastic  hood  that  covers  the  subject  from  head  to 
waist.  The  hood's  vertical  displacement  is  controlled  by  releasing  a  locking 
pin  and  turning  a  cable- connected  crank.  The  diameter  of  the  hood  permits  a 
subject  to  move  freely  his  arms,  shoulders,  and  head.  An  adjustable  cloth  col¬ 
lar  attached  to  the  bottom  edge  of  the  hood  is  drawn  snugly  around  the  sub¬ 
ject's  waist  and  serves  to  contain  the  challenge  atmosphere  [1,3,11,33]. 
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Figure  1,  USAFSAM  sodium  chloride  respirator  quani’Uative  fit  test  instrument. 


Measurement  of  the  Sodium  Chloride  Respirator  Leakage 

When  evaluating  a  full-face  aircrew  chemical  defense  respirator,  a  pri¬ 
mary  concern  is  the  penetration  of  the  sodium  chloride  challenge  atmosphere 
into  the  visual  compartment.  In  order  to  make  this  measurement,  an  aluminum 
tube,  approximately-!  in.  (2.54  cm)  long  and  0.25  in.  (0.635  cm)  i.d.,  is  so 
fitted  through  and  sealed  to  the  protective  device's  visor  that  the  distance 
from  the  cornea  to  the  open  end  of  the  aluminum  sampling  tube  (interior  to  the 
respirator's  visor)  is  not  greater  than  0.8  in.  (2  cm)  [1]. 

The  concentration  of  the  sodium  chloride  challenge  atmosphere  that  has 
leaked  into  the  visual  compartment  is  determined  by  continuously  sampling  gas 
from  this  site  and  vaporizing  it  in  the  flame  photometer.  Sampling  is  accom¬ 
plished  by  attaching  one  end  of  a  short  length  of  Tygon  tubing  to  the  open  end 
of  the  aluminum  sampling  tube  (exterior  to  the  respirator's  visor),  and  then 
passing  the  opposite  end  of  the  plastic  tubing  through  a  sealed  port  in  the 
top  of  the  hood.  A  pump--of  the  metal  bellows  positive  displacement,  continu¬ 
ously  sampling  type--is  connected  to  the  open  end  of  the  hood's  port,  and  is 
used  to  draw  a  gas  sample  from  the  respirator's  visual  compartment  (constant 
7-1  iter /min  flow).  Before  the  7-1  iter /min  sample  is  injected  into  the  hydro¬ 
gen  flame  photometer,  a  6-liter /min  amount  is  bled  to  the  cmbient  atmosphere 
through  a  calibrated  orifice.  The  resulting  1-liter/min  sample  is  diluted 
with  a  14-liter /min  flow  of  clean  dry  compressed  air,  and  this  mixture  is  then 
injected  into  the  photometer  for  analysis.  The  determination  of  the  7-liter/ 
min  mask  sampling  rate,  6-1  iter /min  bleed-off,  and  subsequent  14-1  iter /min 
dilution  of  the  1-liter/min  portion  of  the  penetration  sample  was  based  on  two 
experimental  observations:  First,  the  7-liter /min  sampling  rate  was  selected 
to  minimize  the  negative  pressure  within  the  respirator's  visual  compartment; 
i.e.,  a  greater  sampling  rate  was  observed  to  distort  the  penetration  measure¬ 
ment.  Second,  the  sodium  chloride-air  mixture  composition  was  metered  to 
yield  optimum  performance  of  the  hydrogen  flame  photometer  for  PF's  ranging 
from  10^  to  10°  [33]. 

The  flame  photometer  used  to  analyze  the  sampled  sodium  chloride  aerosol 
from  the  interior  of  a  respirator  consisted  of  two  primary  components:  a 
burner  to  vaporize  the  sodium  chloride  crystals,  and  an  electronic  package  to 
quantify  the  actual  leakage  concentration. 

Hydrogen  is  used  as  the  combustion  gas  for  the  burner  because  of  its 
pure,  almost  colorless  (pale  blue)  flame.  With  the  particular  gas  jet  being 
used,  the  hydrogen  source  is  regulated  to  flow  at  a  nominal  0.05  -  0.07  psi 
(0.34  -  0.48  kPa);  these  conditions  produce  a  1-in.  (2.54  cm)  vertical  flame 
[33]. 


The  electro-optic  scheme  used  to  quantify  a  respirator's  sodium  chloride 
leakage  is  shown  in  Figure  2  [33].  In  operation,  the  sampled  sodium  chloride 
crystals  are  vaporized  in  the  hydrogen  flame  photometer.  The  589  nanometer 
(nm)  wavelength  optical  bandpass  filter  is  used  to  detect  the  yellow  sodium 
emission  lines  and  to  reject  all  other  undesirable  light  energy  emissions. 
Detection  and  quantification  of  the  intensity  of  the  filtered  yellow  light  is 
accomplished  with  a  photomultiplier  tube  (PMT).  In  principle,  the  PMT  output 
current  is  directly  proportional  to  the  intensity  of  the  yellow  light 
impinging  upon  Its  cathode.  Thus,  since  the  Intensity  of  yellow  light  pro¬ 
duced  by  the  hydrogen  flame  Is  proportional  to  the  concentration  of  sodium 
chloride  in  the  medium  surrounding  the  flame,  the  PMT  output  current  Is 
also  proportional  to  the  sodium  chloride  concentration  around  the  flame* 
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processing  scheme  used  in  the  USAFSAM  sodium  chloride  respirator 


Electronic  processing  of  the  PMT  output  current  is  necessary  for  three 
reasons: 

a.  The  PMT  output  current  must  be  converted  to  a  voltage  of  suffi¬ 
cient  magnitude  so  that  the  leakage  can  be  displayed  on  a  strip- 
chart  recorder. 

b.  The  PMT  output  current  must  be  filtered  to  reduce  the  random 
high-frequency  noise  components. 

c.  The  PMT  output  current  must  be  scaled  logarithmically  to  accom¬ 
modate  a  wide  dynamic  range  (four  orders  of  current  magnitude) 
on  a  single  span  strip-chart  recorder. 


The  electronic  processing  scheme  used  in  the  USAFSAM  sodium  chloride  respira¬ 
tor  quantitative  fit  test  instrument  is  shown  in  Figure  3  [33],  A  Harris 
Semiconductor  Corporation  (HA2-2905-5)  integrated  circuit  is  used  as  the 
active  element  in  the  current-to-voltage  converter;  an  Analog  Deviaes  Corpora¬ 
tion  (Model  755N)  integrated  circuit  is  used  as  the  logarithmic  amplifier;  and 
a  Signetics  Corporation  general  purpose  741  operational  amplifier  integrated 
circuit  is  used  in  the  interface  circuit  to  drive  the  strip-chart  recorder. 

The  low- pass  filter  and  logarithmic  amplifier  are  shown  in  Figure  4  [33]. 
Illustrated  in  Figure  5  are  the  electronics  associated  with  the  strip-chart 
recorder  interface  circuit  [33].  The  PMT  detector  is  an  International  Tele¬ 
phone  and  Telegraph  (ITT)  special  purpose  16-stage,  electrostatically  focused, 
model  FW130  tube.  With  an  operating  anode-to-cathode  potential  of  1300-2200 
volts,  the  PMT's  dark  current  is  100  times  less  than  the  PMT  background  cur¬ 
rent  (PMT  background  current  is  defined  to  be  that  current  produced  by  the 
tube  when  the  hydrogen  flame  is  surrounded  by  a  medium  free  of  sodium  chlo¬ 
ride).  With  the  associated  electronics  package  (Figs.  4  and  5),  the  following 
transfer  function  relates  the  PMT  output  current  to  the  strip-chart  recorder 
voltage: 


V  =  2  +  log,^ 


I 


_i _ 

1.1  yA 


} 


(1) 


where 


V  =  output  voltage  to  strip-chart  recorder  (in  volts),  and 

i  =  PMT  output  current  (in  microamperes) 

In  addition,  this  transfer  function  is  illustrated  graphically  in  Figure  6 
[33].  Of  particular  attention  is  the  fact  that  the  PMT  output  current  spans 
11  nanoamperes  (nA)  to  110  microamperes  (yA),  and  the  associated  strip-chart 
recorder  voltage  spans  0  to  4  volts;  i.e.,  a  1-volt  output  change  occurs  per 
decade  of  input  current  change. 

In  order  to  calculate  the  protection  factor  for  a  subject’s  chemical 
defense  respirator,  the  foregoing  information  concerning  the  dynamics  of  the 
instrument  must  be  coupled  with  one  more  factor— the  calibration  procedure. 
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Figure  6.  Transfer  function  plot  relating  PMT  output  current  to  strip-chart  recorder  voltage. 
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Calibration  Procedure 

The  calibration  procedure  being  used  for  the  USAFSAM  sodium  chloride 
respirator  quantitative  fit  test  instrument  utilizes  .the  serial  dilution  tech¬ 
nique  [1].  This  technique  is  based  on  the  fact  that  the  concentration  of 
sodiim  chloride  in  the  challenge  atmosphere  is  directly  proportional  to  the 
concentration  of  the  sodium  chloride  in  the  atomizer  solution.  Thus,  since  a 
5%  sodiun  chloride  solution  standard  is  used  to  produce  the  challenge  atmo¬ 
sphere,  a  series  of  dilutions  of  this  standard  can  be  used  to  produce  known 
PMT  voltage  responses.  If  the  serial  dilutions  are  carefully  selected,  a 
series  of  PMT  voltages,  directly  proportional  to  these  concentrations,  can  be 
measured.  In  the  end,  it  will  be  possible  to  interpolate  between  these  known 
PMT-voltage  response  vs.  serial-dilution  concentrations,  and  to  calculate  the 
soditm  chloride  leakage  (concentration)  (and  thus  PF)  based  on  the  associated 
strip-chart  recorder  voltage  fluctuations.  The  judicious  selection  of  serial 
dilutions  was  made  to  produce  calibration  standards  ranging  from  10"®  to  10"^ 
of  the  basic  5%  sodium  chloride  challenge  atmosphere  solution.  Table  1 
depicts  the  relative  concentrations  and  component  amounts  [33]  for  the  chal¬ 
lenge/calibration  sodium  chloride  dilutions  (see  “Author's  Note,"  below). 


TABLE  1.  SOOIIW  CHLORIDE  CHALLENGE/CALIBRATION  SOLUTIONS  [Ref.  33] 

Relative  Mass  of  sodium  Volune  of  distilled 

concentration  chloride  (grams)  water  (liters) 


10®  (5X  challenge) 

10*1 

50.00 

1.0 

5.000 

1.0 

10*2 

0.5000 

1.0 

10*5 

0.05000 

1.0 

10*4 

0.005000 

1.0 

10-5 

0.0005000 

1.0 

10-® 

0.00005000 

1.0 

When  calibrating  the  instrument,  to  avoid  contaminating  a  weaker  solution 
by  a  stronger  one,  the  operator  begins  with  an  atomizer  containing  the  10"® 
solution  and  advances  to  the  10^  solution.  The  strip-chart  voltage  amplitude 
for  each  calibration  sample  is  annotated  and  retained  for  subsequent  analy¬ 
sis.  A  typical  plot  (semilogarithmic)  of  the  PMT  output  voltage  response  for 
each  of  the  sodium  chloride  concentration  standards  is  illustrated  in  Fig¬ 
ure  7;  and  Table  2  depicts  this  relationship  with  actual  data. 


AUTHOR'S  NOTE:  The  dual  function  of  the  10®  NaCl  solution  is  to  serve  as 
a  calibration  solution  and  the  challenge  solution. 


17 


Sodfutn  chlorfde  cal i brat fon/chal lenge  concentration  vs.  PMT  output  voltage  (semilogarithmic 


TABLE  2.  SODIUM  CHLORIDE  CHALLENGE/CALIBRATION  SOLUTION  CONCENTRATION  VS. 
PMT  OUTPUT  VOLTAGE 


Sodium  chloride  challenge/ 
calibration  solution  concentration 


PMT  output  voltage 
(volts) 


10^  (5%  challenge) 

3.355 

10-1 

2.915 

10-2 

2.310 

10-3 

1.500 

10“^ 

0.545 

10-5 

0.i65 

10-6 

0.105 

CONVENTIONAL  PROTECTION  FACTOR  CALCULATIONS 

A  general  discussion  of  a  protection  factor  is  presented  here,  along  with 
a  description  of  the  conventional  method  (hand  calculation)  used  to  calculate 
the  PF's  associated  with  aircrew  chemical  defense  respirator  quantitative  fit 
testing.  Through  this  information,  the  reader  can  evaluate  the  respective 
advantages  and  disadvantages  of  the  conventional  and  automated  methods  of 
calculating  PF's. 


Protection  Factor 

A  respiratory  protection  factor  is  defined  as  the  ratio  of  the  ambient 
challenge  atmosphere  concentration  external  to  the  respiratory  protective 
device  to  tnat  of  the  sampled  leakage  concentration  drawn  from  the  interior  of 
the  device  [1,2,8,11-20,26,32,34,37,41-43].  Formally,  this  relationship  can 
be  expressed  in  mathematical  terms: 

PF  =  I-  (2) 

where 

PF  =  protection  factor 

Cj  »  ambient  challenge  atmosphere  concentration 
Cj  =  sampled  leakage  concentration 

Note  that  a  PF  is  a  dimensionless  quantity.  In  the  ratio,  the  units  of  con¬ 
centration  in  the  numerator  and  denominator  cancel  (assuming  that  and 
Cj  were  measured  and  appropriately  converted  to  a  consistent  set  of 
concentration  units;  e.g.,  parts  per  million,  micrograms  per  liter,  percent, 
etc.). 


Also  important  in  respirator  quantitative  fit  testing  is  the  calculation 
of  an  average  protection  factor  (PF).  This  calculation  becomes  important  when 
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the  subject  being  evaluated  performs  a  series  of  bre^thinn  anw  u  j 

face-ti?a1eSl^crS:!"‘„ 


I 

i=l  PFi 


where 


PF  =  average  protection  factor  for  n  exercises 
i  =  the  i^H  exercise,  i  =  1.  2.  3.  ...,  n 

PF  =  protection  factor  associated  with  a  particular  exercise 

9JeaJlr"or 

exercise  PF's.  The  most  common  examnip  assigned  to  individual 

exercise  protocol  Is  pfrfoS  fo^7jtfi!l\^^  an 

time  would  become  the  weiqhtina  facfnr  length  of  time;  in  this  case, 

expression  for  an  average  weighted  Jf  is:  ^  mathematical 

n 

_  I  “1  PFl 

DC  _  1=1 


/ 


I  Wi 
i=l 


where 


PFw  “  weighted  average  protection  factor  for  n  exercises 
i  =  the  it^i  exercise,  i  =  1,  2.  3,  ....  n 
Wj  =  weighting  factor  for  the  i^h  exercise 
PF  =  protection  factor  associated  with  a  particular  exercise 

CoTOntioral  Method  of  Colcolatino  o  Protection  Factor 

displlj!  rlco^,  or"“,™"a‘teVV>1"  FyeteMs,  doea  not 

display  the  relative  penetration  fleakJLl  record  and 

calculation  of  PF's  for  the  USAFSAM^tnef^  challenge  atmosphere.  The 

example.  Shown  n  Fiqlrrs  is  I  J^n  r!?  explained  through  an 

recording  that  included  the  Jrelil^  strip-chart 

tion  for  a  set  of  six  exercises;  ^  calibration  and  penetration  informa- 


v'v,.  -.W- 


a.  normal  breathing  (NB) 

b.  deep  breathing  (DB) 

c.  turning  head  side-to-side  vrith  deep  breathing  (TH) 

d.  moving  head  up-and-down  with  deep  breathing  (UD) 

e.  talking  (T) 

f.  facial  grimacing  (FG) 

The  analysis  of  Figure  8  begins  at  the  bottom  of  the  strip-chart  record¬ 
ing.  The  first  section  of  information  uniquely  identifies  the  particular  sub¬ 
ject  and  type  of  respirator.  The  next  section  contains  the  instrument  cali¬ 
bration  data  (a  steady-state  response  measurement  is  recorded  for  each  of  the 
serial  dilutions  and  the  5%  or  10^  sodium  chloride  challenge  concentration). 
Before  the  breathing  exercises  begin,  a  zero  calibration  or  washout  measure¬ 
ment  is  taken  to  establish  a  baseline.  The  six  exercises  follow  in  sequence, 
each  being  performed  for  a  predetermined  time  period. 

The  cyclic  nature  of  the  recorder's  trace  during  the  exercises  is  a 
direct  function  of  the  subject's  breathing  cycle.  Figure  8,  for  instance, 
reveals  that  the  slight  negative  pressure  created  in  the  facepiece  during 
inhalation  increases  the  penetration  of  the  challenge  atmosphere.  Exhalation, 
on  the  other  hand,  creates  a  slightly  positive  pressure,  and  acts  to  reduce 
the  penetration  of  the  challenge  atmosphere.  Because  samples  are  drawn  from 
the  visual  cavity  of  chemical  defense  respirators,  absorption  of  the  sodium 
chloride  aerosol  by  the  lungs  is  negligible  [1].  Therefore,  respirator  per¬ 
formance  is  based  on  the  average  of  the  penetration  peaks  and  valleys  for  each 
of  the  exercises.  Finally,  the  overall  respirator  performance  is  based  on  the 
arithmetic  average  of  the  six  exercise  PF's. 

The  average  of  the  penetration  peaks  and  valleys  (location  of  their  mid¬ 
point)  is  generally  deduced  visually.  By  considering  each  exercise  separate¬ 
ly,  a  line  can  be  drawn  through  the  "visual  average"  of  the  peaks  and  val¬ 
leys.  This  midpoint  line  is  then  extended  until  it  Intersects  the  calibration 
curve  (Fig.  8).  A  plot  of  the  sodiun  chloride  calibration  concentrations  vs. 
their  strip-chart  recorder  displacements  reveals  an  approximate  logarithmic 
relationship  (Fig.  9).  Thus,  logarithmic  interpolation  is  used  to  identify 
the  average  ^)enetrat1on  value  between  adjacent  calibration  decades.  Presented 
in  Table  3  is  a  suninary  of  the  calculations  for  the  example  in  Figure  8. 


TABLE  3.  QUANTITATIVE  FIT  TEST  PENETRATION  RECORD 


j- 


Exercise 


Average  penetration 


Normal  breathing  5.5  x  10'^ 
Deep  breathing  7,/  x  10"^ 
Turning  head  side-to-side  (deep  breathing)  7.2  x  10"® 
Moving  head  up-and-down  (deep  breathing)  5.5  x  10"| 
Talking  5.9  x  10"® 
Facial  grimacing  8.3  x  10"^ 


21 


Figure  9.  Sodium  chloride  logarithmically  scaled  calibration  data  set  used  to  interpolate  mask  leakage 
concentration. 


With  this  information,  the  PF  for  each  exercise  can  he  calculated  using 
the  following  relationship: 


PF 


10°  1 
Cs  Cg 


(5) 


That  is,  since  the  challenge  atmosphere  is  generated  from  the  5%  or  10°  sodium 
chloride  solution  concentration,  and  10°  1,  U  follows  that  a  PF  is  simply 

the  reciprocal  of  the  exercise  average  penetration.  To  complete  this  example. 
Table  4  depicts  the  calculated  exercise  PF's  and  the  overall  arithmetic  aver¬ 
age  PF. 


TABLE  4.  QUANTITATIVE  FIT  TEST  PF  PFCORD 


Exercise  '  PF 


Normal  breathing  1.8  x  10** 
Deep  breathing  1.3  x  10** 
Tuning  head  side-to-side  (deep  breathing)  1.4  x  10** 
Moving  head  up-and-down  (deep  breathing)  1.8  x  10** 
Talking  1.7  x  10* 
Facial  grimacing-  1 .2  x  10^ 


Overall  average  PF  =  1;5  x  lO'* 


Although  the  strip-chart  recorder  data  can  be  interpreted  without 
significant  mathemat'cal  rigor,  this  exercise  can  be  exasperating  when  more 
than  a  dozen  subjects  are  involved.  Having  analyzed  a  series  of  13  tests, 
this  author  has  developed  an  alternative  method  that  utilizes  the  USAFSAM 
POP-11/70  computer  to  perform  these  calculations.  This  automated  scheme 
yields  a  data  reduction  turn-around  time  of  approximately  4  min  per  subject 
vs.  40  min  per  subject  by  manual  calculation. 


USING  A  VOLTAGE-TO-FREQUENCY  CONVERTER 
CIRCUIT  TO  DO  TIME-AVERAGED  INTEGRATION 

Various  techniques  have  been  applied  to  integrate  electronic  signals 
produced  in  the  laboratory.  These  techniques  include  the:  ball  and  disk 
mechanical;  low-inertia  motor;  electrochemical;  analog-to-digital  conversion 
followed  by  counting;  operational  amplifier;  and  voltage- to- frequency  (V/F) 
conversion  followed  by  counting  [63],  The  first  four  techniques  have  been 
evaluated  by  numerous  investigators  for  processing  laboratory- type  recorded 
signals.  Such  Integrators,  however,  possess  certain  Inherent  disadvantages; 
e.g.,  high  cost,  Indirect  readout.  Insufficient  accuracy  and  precision,  and 
Insufficient  dynamic  range. 


Operational  Amplifier  Integrators 


Operational  amplifier  integrators  have  been  used  to  solve  the  short-term 
laboratory  signal  integration  problem  [44-50,53,56,68,62-64],  A  typical  oper¬ 
ational  amplifier  integrator  circuit  is  shown  in  Figure  10  [64J.  The  output 
voltage  of  the  circuit  (Tout)  is  related  to  the  input  voltage  (Ei^)  by  Equa¬ 
tion  6. 


-out 


-  -  -  I 

RC  'o 


in 


dt  + 


^  /  ib  dt 
0  0 


C 


ip  dt 


1_ 

RC 


t 

/  E  dt 
0  os 


(6) 


The  first  tenn  of  Equation  6  represents  the  desired  integrated  value;  and  the 
second  through  fourth  terms  represent  the  output  error  generated  by  inte¬ 
grating  the  input  bias  current  (ib).  the  current  leakage  through  the  inte¬ 
grating  capacitor  (ip^),  and  the  offset  voltage  (Eqs),  respectively.  These 

errors  can  be  minimized  by  restricting  the  period  of  integration  from  1  msec 
to  100  sec  [51 ,52, 54, 55.59, 70],  However,  when  the  integrating  period 
approaches  100  sec  and  an  overall  \t  accuracy  is  desired,  a  large  value,  hiqh- 
perfomance,  expensive  polystyrene  capacitor  is  required,  as  well  as  an 
expensive  operational  amplifier  whose  input  offset  current  and  drift  are 
negligible.  In  order  to  optimize  accuracy,  to  integrate  signals  lasting  from 
milliseconds  to  hours,  and  to  keep  the  cost  of  components  within  manageable 
limits,  the  V/F  integrator  scheme  becomes  an  extremely  attractive  alternative 
[69,71-73], 


Voltage- to- Frequency  Converter  Integrators 

The  basic  function  of  a  V/F  converter  is  to  transform  a  variable  direct- 
current  voltage  (usually  0  to  10  volts)  into  a  pulse  train  whose  re|>etition 
rate  (frequency)  is  a  direct  linear  function  of  the  direct- current  input  volt¬ 
age,  An  excellent  technique  for  precisely  integrating  an  analog  voltage  sig¬ 
nal  is  simply  to  add  a  counter  stage  to  the  output  of  a  V/F  converter  and 
accumulate  the  pulse  count.  By  accimtulating  the  V/F  converter  output  pulses, 
the  "area  under  the  input  voltage  curve,"  or  integral,  is  calculated.  This 
operating  principle  [68,69]  is  illustrated  in  Figure  U  In  addition,  digital 
integrators  using  V/F's  have  three  primary  advantages  over  their  operational 
amplifier  counterparts--for  it  is  very  easy  to: 

a.  simply  switch  the  counter  ^.o  a  "hold"  iwsition,  and 
the  accumulated  count  (integrated  value)  is  held  in¬ 
definitely  with  absolutely  no  drift; 

b.  preset  the  digital  counter  to  any  desired  level  and 
integrate  up  (or  down)  from  that  initial  condition; 
and 

c.  adapt  standard  electro-optical  readouts  (such  as  fiixie 
tubes,  light-emitting  diodes,  etc.)  to  display  the 
integrated  count  value. 

The  attractive  features  of  the  V/F  integration  technique  stimulated  the 
development  of  the  integrator  used  with  the  USAFSAM  R^T  sodiun  chloride 
instrument. 
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operating  principle  of  a  V/F  integrator 


USAFSAM  Sodium  Chloride  RQFT  Voltage- to-Frequency 
Integrator  Circuit  Design 

Shown  in  the  respective  figures  are  the  schematic  diagrams  for:  the 
USAFSAM  sodium  chloride  RQFT  instrument  (Fig.  12);  the  digital  display 
(Fig.  13);  and  the  voltmeter  (Fig.  14).  The  major  components  [66,67,74-82] 
are  listed  in  Table  5. 


TABLE  5.  MAJOR  COMPONENTS  OF  THE  V/F  CONVERTER  ELECTRONIC  INTEGRATOR 


Schematic  Diagram  Ref¬ 
erences  (Figs.  12  -  14) 

Description 

R1 

27K  ohms,  1/4  W, 

5% 

R2 

IK  ohms,  1/4  W, 

5% 

R3 

20K  ohms,  1/4  W, 

5% 

R4 

500-ohm  variable 

potentiometer 

R5 

50K-ohm  variable 

potentiometer 

R6 

lOOK  ohms,  1/4  W, 

5X 

R7 

lOOK  ohms,  1/4  W, 

S% 

R8 

50K  ohms,  1/4  W, 

5X 

R9 

lOOK  ohms,  1/4  W, 

5X 

RIO 

50K  ohms,  1/4  W, 

5t 

Rll 

lOOK  ohms,  1/4  W, 

5» 

RI2 

lOOK  ohms,  1/4  W, 

5X 

R13 

15K  ohms,  1/4  W, 

51 

R14 

lOK  ohms,  1/4  W, 

51 

R15 

IK  ohms,  1/4  W, 

51 

(ContM.  on  facing  page) 


TABLE  5  (Cont'd.) 


Schematic  diagram  ref- 
erences  (Figs.  12  -  14) 

R16 

R17 

Cl 

T1 

ICl 

IC2 

IC3 


IC4 


1C5 


IC6 


IC7 


Digital  Display  Module  - 
Integrator  Count 

Digital  Display  Module  - 
Integration  Time 


Description 


lOOK-ohm  variable  potentiometer 

lOK  ohms,  1/4  watt,  5% 

0.1-uF  capacitor 

2N956  NPN  transistor 

Precision  Monolithics  Incorporated, 
operational  amplifier,  OP- 7 

Analog  Devices  Incorporated,  high  per¬ 
formance  V/F  converter,  AD450J 

Motorola  Semiconductor  Products  Incorpo¬ 
rated,  decade  counter/divider, 

MC14017B 

Motorola  Semiconductor  Products  Incorpo¬ 
rated,  noninverting  hex  buffers,  MC 
14050B,  Vcc  -  Pin  1.  Vjs  -  Pin  8, 
Ground  Pins  11  and  14,  Vqc  “  +5  volts 

Motorola  Semiconductor  Products  Incorpo¬ 
rated,  quad  2- input  OR  gate,  MC14071B, 
Vcc  ■  Pi'’  Vss  ”  Pin  7,  Ground  Pins 
5,  6,  7,  8,  9,  12,  and  13,  V^c  ®  +5 
volts 

Texas  Instruments  Incorporated,  hex 
inverter,  SN7404,  Vrr  -  Pin  14,  Vcc  - 
Pin  7,  Vcc  =  +5  volts 

Motorola  Semiconductor  Products  Incorpo¬ 
rated,  industrial  time  base  generator, 
MC14566B 

Dialight  LEO  display  module  749-1706 


Oialight  LED  display  module  749-1704 
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Figure  12:  Part  2  (of  3).  Sodium  chloride  RQFT  instrument  integrator. 


instrument  integrator. 


Figure  13.  Sodium  chloride  RQFT  instruii’°nt  integrator  digital  display. 

(♦Delete  these  two  digits'from  the  module  for  the  4-digit  LED  time  display.) 


Figure  14.  Sodium  chloride  RQFT  instrument  integrator  voltmeter. 


Description  of  the  Analog  Devices  AD450J 
V/F  Converter  Integrated  Circuit 

The  heart  of  the  USAFSAM  sodium  chloride  RQFT  instrument  integrator  design 
is  the  Analog  Devices  AD4500  V/F  converter.  The  electronic  specifications  for 
this  device  [66,67,78]  are  summarized  in  Table  6,  and  a  block  diagram  of  the 
converter  is  shown  in  Figure  15  [66,67,78]. 

The  AD450J  is  a  low-cost  high-performance  V/F  converter  that  provides 
exceptional  linearity  and  temperature  stability  over  a  wide  dynamic  input  sig¬ 
nal  range.  The  key  to  the  AD450J‘s  precise  operation  is  the  unique  charge 
balance  conversion  technique.  This  feature  means  that  an  analog  signal  is 
accurately  converted  to  a  train  of  pulses  (each  of  which  has  a  constant  width 
and  amplitude)  at  a  rate  directly  proportional  to  the  analog  signal  ampli¬ 
tude.  The  output  continuously  tracks  and  responds  directly  to  changes  in  the 
input  signal. 

The  versatile  operational  amplifier  buffer  serves  as  the  input  stage.  Its 
purpose  is  to  convert  the  applied  input  voltage  to  a  control  current  for  the 
charge  balance  conversion  circuitry. 

External  clock  synchronization  is  not  required  for  the  AD450J.  The  inter¬ 
nal  temperature  compensated  timing  reference  yields  an  accurate  square-wave 
output  pulse  train  void  of  low-frequency  cycle-to-cycle  jitter. 

Finally,  the  output  drive  circuitry  is  responsible  for  conditioning  the 
stable  output  pulse  train.  The  high  current-handling  capacity  of  this  cir¬ 
cuitry  allows  the  designer  to  interface  directly  with  any  low-cost  digital 
processing  logic  family.  These  attractive  features  motivated  the  selection  of 
the  AD450J  as  the  V/F  converter  circuit  for  the  long-term  precision  integrator 
design. 


Operation  of  the  USAFSAM  Sodium  Chloride  RQFT 
V/F  Integrator 

Operation  of  the  USAFSAM  sodium  chloride  RQFT  V/F  integrator  is  quite  sim¬ 
ple.  The  integrator's  capability  and  function  can  be  appreciated  by  analyzing 
Figure  12.  Only  one  connection  to  the  integrator  is  required.  The  analog 
voltage  signal  to  be  integrated  (the  input  signal  to  the  strip-chart  recorder) 
is  connected  directly  to  the  0-10  volt  integrator  input  terminal. 

The  integrator  is  energized  with  the  on/off  switch,  and  the  front  panel 
neon  lamp  yields  a  positive  indication  for  applied  power.  Two  direct-current 
power  supplies  (±15  volt  and  +5  volt)  provide  the  operating  voltages  for  the 
integrator's  digital  circuits. 

Because  the  output  signal  of  the  sodium  chloride  RQFT  instrument  is  biased 
with  a  low-level  noise  component,  a  Precision  Monolithics  Corporation  OP-7 
operational  amplifier  buffer  is  used  to  compensate  for  this  unwanted  signal 
[74].  The  degree  of  compensation  is  accomplished  by  rotating  the  PMT  noise 
offset  voltage  adjustment  potentiometer  (R16).  This  adjustment  must  be 
checked  each  time  the  RQFT  instrument  is  turned  on. 


TABLE  6.  ELECTRONIC  SPECIFICATIONS  OF  THE  ANALOG  DEVICES  AD450J  V/F  CONVERTER 


Characteri Stic 

Value 

TRANSFER  FUNCTION 

Voltage  Input 

fout  =  (10’ 

ANALOG  INPUT 

Voltage  Signal  Range  (ein) 

0  to  +10V  min 

Overrange 

50%  min 

Impedance  (eif,) 

20kfl 

Max  Safe  Input  Voltage  (eip) 

+25V  (-Vg) 

ACCURACY 

Warmup  Time 

1  min 

Nonl inearity 

e^n  =  +lmV  to  +15V 

±0.01%  max 

Full  Scale  Error 

(+0.5  +  1.5)%  max 

Gain 

vs.  Temperature  (0  to  70®C) 

±50ppni/°C  max 

vs.  Supply  Voltage 

±200ppm/%max 

vs.  Time 

±100ppm/day 

Input  Offset  Voltage 

±5mV  max 

vs.  Temperature  (0  to  70°C) 

+50wV/»C 

vs.  Supply  Voltage 

±10ppm/Xmax 

vs.  Time 

±10uV/day 

(Cont'd.  on  facing  page) 


TABLE  6  (Cont'd.) 
Characteristic 


Value 


RESPONSE 

Settling  Time  for  +10V  Step  Input 
Overload  Recovery  Time 

OUTPUT 
Waveform 
Pulse  Width 
Rise/Fall  Time 
Pulse  Polarity 
Logic  "1"  (High)  Level 
Logic  "0"  (Low)  Level 
Capacitive  Loading 
Fan  Out  Loading 
Impedance 

POWER  SUPPLY 

Voltage,  Rated  Performance 
Voltage,  Operating 
Current,  Quiescent 

TEMPERATURE  RANGE 
Rated  Performance 
Operating 
Storage 

CASE  SIZE 


120yS 

15ms 

Train  of  TTL/DTL  compatible  pulses 
50us 
200ns 
positive 
+2.4V  min 
+0.4V  max 
lOOOpF  max 
10  TTL  loads  min 
3.3kn 

tl5V  dc 

±(12  to  18)V  dc 
(+15,  -9)mA 

0  to  +70'‘C 
-aS^C  to  +80‘’C 
-55“C  to  ♦85'’C 

(1.5  x  1.5  X  0.4)  inches 
transistor- transistor  logic. 


OTL  “  diode- transistor  logic;  and  TTL 
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g  Devices  A0450J  V/F  converter 


The  overall  accuracy  and  dynamic  range  of  the  AD450J  V/F  converter  is 
established  with  an  initial  calibration  of  the  two  trim  potentiometers  (R4  and 
R5)  shown  in  Figure  12.  This  calibration  is  accomplished  after  the  integrator 
has  reached  an  operating  equilibrium  {5-min  warmup).  A  precision  regulated 
voltage  source  of  +1.00  mV  is  connect^  to  pin  6  ofICl,  and  kn  oscilloscope 
is  connected  between  ground  and  pin  6  of  IC2.  The  zero  offset  adjustment 
potentiometer  (R5)  is  rotated  so  that  a  1  Hz  frequency  pulse  train  is  visible 
on  the  oscilloscope.  The  precision  regulated  voltage  source  is  then  set  to 
+10.000  V,  and  the  full  scale  adjustment  potentiometer  (R4)  is  rotated  so  that 
a  10  kHz  frequency  pulse  train  is  visible  on  the  oscilloscope. 

The  analog  input  voltage-to-frequency  conversion  process  is  accomplished 
with  the  Analog  Devices  AD450vJ  V/F  converter.  Its  output  pulse  train  is  pro¬ 
cessed  by  the  Motorola  Semiconductor  Products  MC14017B  decade  counter-divider 
integrated  circuit.  This  integrated  circuit  conditions  and  accumulates  the 
V/F  converter  output  signal  for  the  6-digit  LED  integrator  count  display.  The 
user  can  select  a  signal  dividing  constant  via  the  panel -mounted  integrator 
count  scale  switch.  Selection  of  the  divide-by-one  position  means  a  one-co- 
one  correspondence  exists  between  the  actual  integrator  count  (area  under  the 
analog  voltage  input  signal)  and  the  magnitude  on  the  digital  display.  Selec¬ 
tion  of  any  switch  position  2  through  9  means  that  the  magnitude  on  the  digi¬ 
tal  display  must  be  multiplied  by  the  switch  position  number  to  yield  the 
actual  integrator  count.  This  feature  permits  the  user  to  integrate  input 
signals  for  several  hours  without  overloading  the  integrator  count  display. 

An  internal  time  base  generator  is  an  inherent  part  of  the  integrator 
design.  A  Motorola  Semiconductor  Products  MC14566  integrated  circuit  is  used 
to  provide  a  clocked  output  pulse  each  second.  The  accumulated  time  (in  sec¬ 
onds)  is  available  via  the  4-digit  LEO  display.  The  user  controls  the  length 
of  time  an  analog  input  signal  is  integrated.  Two  front-panel  mounted 
switches  give  the  user  the  flexibility  to  clear  both  LED  displays  to  read  zero 
and  halt  the  integration  of  a  signal  at  any  point  in  time.  To  clear  the  dis¬ 
plays,  the  user  simply  toggles  the  display  reset  switch  and  trips  the  integra¬ 
tor  functior\  switch  to  the  hold  position.  When  the  operator  is  ready  to  inte¬ 
grate  a  signal,  he  toggles  the  integrator  function  switch  to  the  run  position. 
Thus,  the  signal  will  be  integrated  and  the  integrator  count  and  elapsed  time 
will  be  automatically  recorded  on  the  appropriate  LEO  displays.  When  the 
integration  is  to  be  halted,  the  usa'  toggles  the  integration  function  switch 
to  the  hold  position.  This  last  operation  gWes  the  user  the  opportunity  to 
record  the  magnitude  of  the  integrator  count  and  elapsed  time.  When  another 
signal  is  to  be  integrated,  the  user  resets  the  displays  to  zero  and  repeats 
the  sequence  described  above. 


Data  Collection  with  the  US/tfSAM  Sodium 
'  Chloride  RQFT  Instrument  V/F  Integrator 

The  collection  of  RQFT  data  for  subsequent  PF  calculations  is  a  simple 
process  when  the  integrator  is  used.  Shown  in  Figures  16  and  17  are  the  data 
sheets  used  for  this  purpose.  After  the  RQFT  Instrument  has  reached  its  oper¬ 
ating  equilibrium,  and  the  integrator's  PHT  noise  offset  voltage  has  been 
adjusted,  the  sodium  chloride  calibration  and  challenge  atmosphere  atomizer 
solutions  are  processed. 
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RESPIRATOR  QUANTITATIVE  FIT  TESTING 
--USAFSAM  SALT  FOG  INSTRUMENTATION- 


Concentratlon 

10  to  the  Zero 
10  to  the  Minus  One 
10  to  the  Minus  Two 
10  to  the  Minus  Three 
10  to  the  Minus  Four 
10  to  the  Minus  Five 
10  to  the  Minus  Six 


SUBJECT  NAME: 
TYPE  OF  MASK: 
DATE  TESTED: 
TIME  TESTED: 


Exercise 


Normal  Breathing  Straight  Ahead 
Deep  Breathing  Straight  Ahead 
Talking 

Side-to-Side  Head  Movements  (Deep  Breathing) 
Up-and-Down  Head  Movements  (Deep  Breathing) 
Facial  Grimacing 


Count  Time  Period 

(In  Seconds) 


Figure  16.  Sodium  chloride  RQFT  data  collection  form  No.  1. 
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RESPIRATOR  QUANTITATIVE  FIT  TESTING 
-USAFSAM  SALT  FOG  INSTRUMENTATION- 


Concentratlon 

10  to  the  Zero 
10  to  the  Minus  One 
10  to  the  Minus  Two 
10  to  the  Minus  Three 
10  to  the  Minus  Four 
10  to  the  Minus  Five 
10  to  the  Minus  Six 


Voltage  (In  Volts) 


SUBJECT  NAME: 
TYPE  OF  MASK: 
date  TESTED: 
TIME  TESTED: 


Exercise 


Count  Time  Period 

(in  Seconds) 


Normal  Breathing  Straight  Ahead 

Normal  Breathing  Left 

Normal  Breathing  Right 

Normal  Breathing  Down 

Normal  Breathing  Up 

Deep  Breathing  Straight  Ahead 

Deep  Breathing  Left 

Deep  Breathing  Right 

Deep  Breathing  Down 

Deep  Breathing  Up 

Talking 

Fa.ial  Grimacing 

Side-to-Side  Head  Movements  (Normal  Breathing) 
Up-and“0own  Movements  (Normal  Breathing) 
Side-to-Side  Head  Movements  (Deep  Breathing) 
Up-and-Down  Head  Movanents  (Deep  Breathing) 


Figure  17.  Sodium  chloride  RC^T  data  collection  form  No.  2. 
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The  data  collected  for  each  of  the  standards  is  the  steady-state  output  voltage 
i  displayed  on  the  integrator's  digital  voltmeter.  These  voltages  will  range 

:  from  a  maximum  of  4.0  volts  to  a  minimum  of  0.09  volts;  each  reading  is  re- 

1  corded  opposite  the  sodium  chloride  concentration  (Fig.  16).  The  next  block  of 

I  information  to  be  recorded  is  the  subject  and  respirator  identification  infor- 

I  mation.  After  each  exercise  is  accomplished,  the  integrator  count  on  the 

\  6-digit  LED  display  and  time  on  the  4-digit  LED  display  are  recorded.  The 

;  elapsed-time  display  is  used  to  initiate  and  terminate  each  exercise.  Before 

proceeding  to  the  next  exercise,  the  integrator  displays  are  reset  to  zero. 

After  all  subjects  have  been  tested,  the  user  proceeds  to  a  computer  term¬ 
inal  and  enters  the  information  from  the  data  collection  forms.  The  interac¬ 
tive  curve-fitting  and  PF  calculation  program  described  in  the  next  section  is 
used  to  process  the  RQFT  data. 


LEAST  SQUARES  CURVE  FITTING  COMPUTER  PROGRAM 
TO  CALCULATE  PROTECTION  FACTORS 

The  fitting  of  empirical  data  by  formulas  or  equations  can  be  accomplished 
by  two  methods.  The  first  is  to  have  a  polynomial  that  is  satisfied  exactly  at 
the  observed  data  points;  this  is  commonly  referred  to  as  "the  polynomial 
interpolation  method"  [83,84,86-92].  The  second  method,  however,  is  a  more 
desirable  way  to  analyze  data  which  have  been  gathered  from  experimental  obser¬ 
vations  that  are  biased  with  various  errors  of  measurement;  this  is  commonly 
referred  to  as  "the  least  squares  approximation  method"  [83-87,93].  Because 
the  measurement  errors  associated  with  the  RQFT  calibration  data  have  been 
empirically  determined  to  be  relatively  small  in  magnitude,  a  least  squares 
curve  fitting  algorithm  is  used  to  determine  an  interpolating  polynomial.  The 
theory  of  fitting  nonlinear  curves  by  the  method  of  least  squares  yields  an 
interpolating  polynomial  when  the  degree  of  the  polynomial  is  one  less  than  the 
number  of  data  points.  Additionally,  this  curve  fitting  polynomial  is  unique 
and  the  data  points  are  fully  parameterized  [96].  A  summary  of  the  mathemati¬ 
cal  theory  for  the  method  of  least  squares  curve  fitting  is  presented  first, 
and  is  followed  by  its  direct  application  to  orocess  the  RQFT  data  and  calcu¬ 
late  PF's. 


The  Method  of  Least  Squares  Curve  Fitting 

The  objective  of  the  method  of  least  squares  curve  fitting  is  to  relate  by 
some  function,  y  «  f(x),  a  set  of  m  points  (xj,  yj),  (j  »  1,  2,  3,  ...,  m), 
which  have  been  gathered  through  some  measuring  process.  The  method  of  least 
squares  curve  fitting  assumes  that  the  function,  y  “  f(x),  can  be  written  as  a 
polyncx'aial  of  degree  n  <  m; 


y  «  +  ajx  +  ajx^  +  ...  anx^  »  j  a^x'* 

i»0  ^ 


(7) 
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The  next  step  in  the  process  is  to  determine  the  value  of  the  coefficients 
a-j,  (i  =  0,  1,  2,  n) ,  such  that  the  polynomial  described  by  Equation  7 
is  a  "good  fit"  to  the  data  (xj,  yj).  Ry  substituting  the  data  points 
into  the  polynomial  described  oy  Equation  7,  a  set  of  m  simultaneous  equations 
are  generated: 


Ri  =  ao  +  a^xi  +  a2Xi^  +  ...  +  apXi^  -  y^ 

o  n 

R2  "  3o  +  ^1^2  ^  ^2^2  ^  •••  ^n^2  "  ^2 


(8) 


Rm  =  ap  ^2^  +  ...  +  anXm  -  ym 


These  equations  are  not  exactly  equal  to  zero,  because  the  polynomial  does  not 
necessarily  pass  through  all  the  points  (except  in  the  case  where  the  degree 
of  the  polynomial  is  one  less  than  the  number  of  data  points).  As  shown  in 
Figure  18,  the  differences  between  a  polynomial  value  and  a  data  value  can  be 
positive,  negative,  or  zero.  This  difference  is  called  a  residual.  Residuals 
are  readily  calculated  by  Equation  9. 


Rj  =  I  (J  =  2.  m)  (9) 
i=0  ^ 

Thus,  the  set  of  equations  described  by  Equation  8  can  be  referred  to  as  "the 
residual  equations."  The  principle  of  least  squares  curve  fitting  states  that 
the  best  representation  of  the  data  is  that  which  makes  the  sum  of  the  squares 
of  the  residuals  a  minimum  [83,86,87,93].  Therefore,  it  is  desirable  to  force 
the  function 


f(ao,  ai,  a2,  ....  an)  “  Ri*  +  Rz^  ♦  +  ...  Rm^  (10) 


to  be  as  close  to  zero  as  possible.  The  condition  which  fulfills  this 
requirement  is  that  the  partial  derivatives  of  Equation  10  be  exactly  zero 
[85,94]. 


Substituting  the  results  of  Equation  12  into  Equation  11  yields: 


Ri  +  R2  +  R3  +  ...  Rm  “  0 


XiRi  +  X2R2  ^3^3  +  ...  +  X[i]Rni  =  0 
Xi^Ri  +  X2^R2  +  X3^R3  +  ...  +  Xm^Rm  =  0 


(13) 


Xi^iRi  +  X2^R2  +  X3'1R3  +  ...  +  Xm^Rn,  =  0 


Replacing  the  Rj's  by  their  values,  defined  in  Equation  8,  and  collecting 
the  coefficients  of  the  (n+1)  unknowns  of  a-j,  (i  =0,  1,  2,  3,  ...»  n), 
yields: 


mao 

+  Ixjai 

* 

+  ...  +  Ixj^an 

-  bi  0 

Ixjao 

*  1^3^32 

+  ...  +  Ixji^+^an 

“  0 

Ixj^ao 

-  Ix/a; 

+  ...  + 

-  Ixj^Yj.  =  0 

+  i;xjf>+‘av  +  +  ...  +  Ixj2%  -  Ixjfiyj  =  0 

where  all  sumnatlons  are  from  1  to  m;  that  is, 

IXj®  «  Xi^  +  X2^  +  Xj^  +  ...  +  x„,^ 

“  xi^yi  +  ♦  ...  +  Xm^Ym 


(14) 


(15) 


The  set  of  equations  defined  in  Equation  14  are  known  as  the  normal  equa¬ 
tions.  All  of  the  associated  summations  are  known;  so  the  system  of  equations 


defined  by  Equation  14  is  a  system  of  (n+1)  linear  equations  in  the  (n+1) 
unknowns  of  a-j  (i  =  0,  1,  2,  ...»  n).  The  solution  of  Equation  14  yields 
the  coefficients,  a^-,  and  thus  the  polynomial  defined  by  Equation  7  is 
determined. 

The  principle  of  least  squares  curve  fitting  is  not  limited  to  polynomi¬ 
als  defined  solely  in  terms  of  x.  The  functional  relationship  between  y  and 
X  can  be  any  known  form  (i.e.,  e^,  logjgX,  sin  x,  1/x,  etc.)  as  long  as  the 
functional  form  is  defined  and  the  resulting  normal  equations  can  be  solved. 

The  simplest  case  by  far,  however,  is  for  y  defined  directly  in  terms  of  x 
[83-87,93]. 

The  coefficients  for  Equation  7  are  found  by  calculating,  from  the  mea¬ 
sured  data,  all  the  sums  defined  by  Equation  14.  These  sums  are  then  substi¬ 
tuted  into  this  system  of  equations,  and  the  coefficients  a^  (i  =  0,  1,  2, 

....  n)  are  thus  determined. 

Application  of  the  Method  of  Least  Squares 
Curve  Fitting  to  Calculate  RQFT  PF's 

The  method  of  least  squares  curve  fitting  was  applied  to  process  the  cal¬ 
ibration  and  exercise  integrator  count  data;  the  final  product  is  the  calcula¬ 
tion  of  PF's.  The  author  wrote  a  Fortran  computer  program  that  accepted  the 
calibration  data  and  calculated  a  least  squares  polynomial  curve  fit  equa¬ 
tion.  The  integrator  count  data  for  each  exercise  was  then  substituted  into 
the  polynomial  equation,  and  a  corresponding  mask  leakage  penetration  (concen¬ 
tration)  was  calculated.  Finally,  PF's  were  calculated  using  Equation  5.  The 
subsequent  paragraphs  of  this  report  outline  the  analysis  used  to  accomplish 
this  method  of  calculating  PF's. 

The  objective  of  writing  a  least  squares  curve  fitting  computer  program 
was  to  take  advantage  of  the  intrinsic  importance  of  the  calibration  data  and 
relate  this  set  of  seven  points  (xj,  yj),  (j  =  1,  2,  3,  ...,  7)  by  some 
function  y  =  f(x).  In  this  particular  case,  the  x  coordinate  is  the  digital 
voltmeter  reading  (in  volts),  and  the  y  coordinate  is  the  corresponding 
sodium  chloride  calibration  or  challenge  concentration.  (This  determination 
was  made  to  avoid  introducing  additional  error  in  the  calculation  of  a  PF  by 
having  to  implement  a  complex  inverse  iteration  algorithm.)  By  referring  to 
Figure  16,  the  reader  can  observe  that  these  seven  ordered  pairs  correspond  to 
the  data  at  the  top  of  the  RQFT  form,  A  typical  set  of  calibration  data  are 
•^hown  in  Table  7. 

After  several  sets  of  calibration  data  had  been  examined,  and  various 
functional  definitions  of  the  variable  had  been  tried,  it  was  determined  that 
a  polynomial  of  the  following  form  yielded  the  best  and  most  stable  fit: 

y  =  ao  +  ajex  +  a2(e>«)2  +  +  a4(ex)‘*  +  agfe^^)®  +  agCe^)®  (16) 

The  "goodness"  of  fit  was  checked  utilizing  the  conditional  Eq.  10  to  verify 
that  the  polynomial  indeed  passed  through  the  RQFT  calibration  data  points. 

The  computer  program  has  calculated  the  "goodness"  of  fit;  and  the  results  for 
a  typical  run,  using  the  data  in  Table  7,  are  presented  in  Table  8.  In 
addition,  a  computer- generated  (Cal comp)  plot  (semi logarithmic)  is  shown  in 
Figure  19. 
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TABLE  7.  TYPICAL  SODIUM  CHLORIDE  RQFT  CALIBRATION  DATA  USED 
FOR  LEAST  SQUARES  CURVE  FITTING 


X-coordinate 
(voltage  reading) 

[in  volts] 

Y-coordinate 
(sodium  chloride 
concentration) 

3.380 

1.0  . 

2.915 

0.1 

2.310 

0.01 

1.500 

0.001 

0.545 

0.0001 

0.165 

0.00001 

0.105 

0.000001 

TABLE  0.  LEAST  SQUARES  CURVE  FIT  CALCULATIONS 
FOR  THE  DATA  CONTAINED  IN  TABLE  7 


The  order  of  the  desired  polynomial  =  6 

The  polynomial  functional  definition  of  the  variable  (X)  is  in  terms  of: 
(Exponential  (X)) 


Coefficient 

Number 

1 

s 

-4.17904E-05 

Coefficient 

Number 

2 

s 

-7.31487E-05 

Coefficient 

Number 

3 

s 

1.25863E-04 

Coefficient 

Number 

4 

s 

-2.68095E-05 

Coefficient 

Number 

5 

= 

3.80138E-06 

Coefficient 

Number 

6 

s 

-1.71541E-07 

Coefficient  Number 

7 

= 

3.88406E-09 

1  residuals  are 

calculated  by  the  following  equation: 

[Y(DATA)-Y(ESTIMATED)] 

SAMPLE  NUMBER 

Y(DATA) 

Y(ESTIMATED) 

RESIDUAL 

1 

1.00000 

1.00000 

-0.483548E-07 

2 

0.100000 

0.100000 

-0.189111E-07 

3 

O.lOOOOOE-01 

O.lOOOOOE-01 

-0.711342E-08 

4 

O.lOOOOOE-02 

O.lOOOOOE-02 

-0.401929E-08 

5 

0.100000E-03 

0.100004E-03 

-0.390294E-08 

6 

O.lOOOOOE-04 

0.100034E-04 

-0.340602E-08 

7 

O.lOOOOOE-05 

0.100433E-05 

-0.432502E-08 

Sum  of  the  square  residuals  =  2.80811E-15 
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From  the  information  in  Table  8,  the  polynomial  curve  fit  equation  is: 

y  =  (-4.17904  x  lO'^)  +  (-7.31487  x  10-^)eX  +  (1.35863  x  10-'' (eX)^  + 
(-2.68095  X  10-5)(eX)^  +  (3.80138  x  10-^)(e>')''  + 

(-1.71541  X  10-’)(eX)^  +  (3.88406  x  10-^)(eX)®  (17) 


The  residuals  were  calculated  to  measure  the  "goodness  of  fit"  for  the  poly¬ 
nomial  curve  fit  Equation  17.  Since  the  RQFT  data  points  are  fully  parameter¬ 
ized  by  Equation  17,  the  residuals  and  sum  of  the  square  residuals  should  be 
exactly  zero.  However,  these  values  are  not  exactly  equal  to  zero  because  of 
computer  calculation  round-off  errors.  The  computer  program  made  these  calcu¬ 
lations  by  substituting  the  set  of  sodium  chloride  calibration/challenge  con¬ 
centration  (Y-coordinate  values  of  Table  7)  values  into  the  polynomial  curve 
fit  Equation  17;  the  corresponding  calculated  Y  values  are  called  the  Y-esti- 
mated  values.  One  can  readily  observe  that  the  sum  of  the  square  residuals  is 
2.80811  X  10"^^,  a  very  small  number  indeed.  Interestingly,  the  worst-case 
residual  occurs  for  the  10"^  sodium  chloride  calibration  standard;  this 
finding  implies  that  the  accuracy  of  predicting  a  particular  exercise  mask 
leakage  concentration  will  not  be  in  error  by  more  than  0.433%  [(0.100433  x 
10-^  -  0.100000  X  10-^)7(0.100000  X  10-^)  =  0.00433  x  100].  The  magnitude  of 
this  error  is  representative  of  that  found  for  several  sets  of  calibration 
data  that  have  been  analyzed. 

The  following  explanation  is  the  key  to  understanding  how  a  best  fit 
polynomial  curve  (relating  sodium  chloride  calibration/challenge  concentration 
to  a  digital  voltmeter  response)  can  be  used  to  calculate  a  PF.  Two  concepts 
are  involved  in  the  analysis. 

First:  the  Y-axis  is,  in  reality,  an  exact  scale  for  the  time-averaged 
mask  leakage  penetration  (concentration).  When  the  user  calibrates  the  RQFT 
instrument,  the  strip-chart  recorder  displays  (on  its  vertical  or  Y-axis)  the 
PMT's  response  for  a  known  diluted  sodium  chloride  concentration  sample.  It 
is  important  to  recognize  that  the  response  for  the  diluted  concentration  is 
equivalent  to  that  which  would  be  measured  if  one  were  to  evaluate  the  fit  of 
a  mask  with  the  10°  sodium  chloride  challenge  concentration,  and  have  the  mask 
fit  be  perfect,  except  for  a  "calibrated"  leak  that  would  permit  a  correspond¬ 
ing  known  diluted  concentration  of  the  10°  challenge  to  penetrate  the  mask. 

The  essence  of  the  polynomial  curve  fit  equation  is  that  it  permits  "interpo¬ 
lation"  between  the  responses  for  adjacent  sodium  chloride  calibration  con¬ 
centrations. 

Second:  the  integrator  count— the  source  of  data  used  to  calculate  a 
mask  leakage  penetration  (concentration)-- is,  in  reality,  a  time-averaged 
voltage  response.  This  fact  can  be  derived  through  the  following  analyses: 

a.  Integrator  sensitivity  Is  1000  counts-per-volt-sec. 

b.  Each  exercise  is  performed  for  a  predetermined  length  of  time;  for 
example,  10  sec. 

c.  The  Integrator  count  (IC)  value  recorded  for  a  particular  exercise 
Is  actually  the  time-averaged  area  under  the  strip-chart  recorder 
response  (refer  to  Figs.  8  and  11). 


Therefore, 


in  which  V  volts  is  the  time- averaged  voltage  for  a  particular  exercise. 

The  computer  program,  written  to  accept  the  integrator  exercise  count 
data  and  the  time  duration  for  each  exercise,  calculates  corresponding  exer¬ 
cise  time-averaged  voltages.  Fach  time-averaged  voltage  is  substituted  into 
the  polynomial  curve  fit  Equation  17,  and  a  corresponding  mask  leakage  pene¬ 
tration  (concentration)  is  calculated.  The  final  set  of  calculations  per¬ 
formed  by  the  program  are  the  individual  exercise  PF‘s,  an  arithmetic  average 
PF,  and  a  time-weighted  average  PF.  Equation  5  is  used  to  calculate  individ¬ 
ual  PR's,  and  Equations  3  and  4  are  used  to  calculate  the  average  PR's.  The 
results  of  the  data  in  Tables  7  and  8  are  given  in  Table  9- 

Discussion  of  the  Computer  Programs  Used  to  Process  the 
USAFSAM  Sodium  Qtloride  RQFT  Integrator  Data 


Two  computer  programs  are  used  to  process  the  integrator  exercise  data. 
The  first  and  primary  program  is  called  NACLRQRT.RTN;  the  second, 
NACLGRAPH.RTN. 


The  purpose  of  NACLRQRT.RTN  is  to  use  the  RQFT  information  collected  or; 
the  RQFT  data  sheet  (Fig.  16)  ar>d  calculate  a  set  of  PR's.  The  results  of 
this  program  are  stored  on  three  disk  files: 

1.  DATA. XXX  contains  the  initial  calibration  data,  the 
test  identification  data,  the  time  period  for  each 
exercise,  and  a  listing  of  the  exercises  performed  and 
their  associated  integrator  count  values. 


CALCX.XXX  contains  the  polynomial  curve  fit  equation 
coefficients,  the  calculated  residuals,  the  sum  of  the 
square  residuals,  a  composite  listing  of  the  identifi¬ 
cation  data,  the  exercises  performed  and  their  corre¬ 
sponding  PF's,  and  the  average  PR's. 

GRPHX.XXX  contains  an  array  of  401  X-axis  values  and 
401  Y-axis  values.  These  values  were  generated  using 
the  polynomial  curve  fit  Equation  15.  This  array  Is 
used  by  the  NACL0RAPH.FTN  program  to  generate  a 
CALCOHP  calibration  curve  for  the  sodiun  chloride  RQFT 
cal ibration  data. 


TABLE  9.  PROTECTION  FACTOR  CALCULATIONS  FOR  THE  DATA  CONTAINED  IN 
TABLE  7 


SUBJECT  NAME:  Captain  Edward  S.  Kclesar,  Jr. 
TYPE  OF  MASK:  USA:  M17  -  Medium  (no  glasses) 
DATE  TESTED:  9  April  1980 
TIME  TESTED:  1330  hours 

EXERCISE  INTEGRATOR  COUNT  DATA 


Exercise 

Integrator 

count 

Time  period 

Protection  factor 

Normal  breathing  straight  ahead 

3904 

10 

1 .8E+04 

Deep  breathing  straight  ahead 

4751 

10 

1.3E+04 

Talking 

4628 

10 

1 .3E+04 

Side-to-side  head  movements 
(deep  breathing) 

3976 

P 

l./:+04 

Up-and-down  head  movements 
(deep  breathing) 

40  i' 6 

10 

1 .7E+04 

Facial  grimacing 

4937 

10 

Overall  arithmetic  average  protection  factor  (PF)  for  all  categories  of 
exercises  actually  performed  <»  1.5E+04. 

Ovcjrall  time-weighted  average  protection  factor  (PF)  for  all  categories  of 
exercises  actually  perfor^  •  l,5£*04. 


The  NACLGRAPH.FTN  program  has  been  written  as  a  separate  program  because 
the  USAFSAM  PDP-11/70  CALCOMP  plotter  is  an  off-line  device.  That  is, 
NACLGRAPH.FTN  uses  the  GRPHX.XXX  disk  file  as  input  data  and  sequences  through 
the  CALCOMP  plot  subroutine  library  files,  and  produces  a  user-named  output 
disk  file  (for  example,  PLOT. SCO).  The  NACLGRAPH.f  "N-generated  output  file  is 
then  transferred  to  a  magnetic  tape  by  one  of  the  computer-room  operators. 

The  magnetic  tape  is  then  mounted  on  the  CALCOMP  terminal  and  the  plot  is 
generated  [95]. 

Each  of  the  programs  discussed  in  this  section  is  documented  with  com¬ 
ments  that  define  the  variables  and  explain  the  operations  performed.  There¬ 
fore,  a  line-by-line  analysis  of  the  code  will  not  be  done.  For  the  inter¬ 
ested  reader,  however,  the  following  information  is  available  in  eight  appen¬ 
dixes  (A  -  H) 

Appendix  A:  NACLRQFT.FTN  Fortran  listing 

Appendix  B:  OATA.XXXX  file  contents  for  data  in  Table  7 

Appendix  C:  CALCX.XXX  file  contents  fcr  information 
in  Tables  8  and  9 

Appendix  0:  GRPHX.XXX  file  contents  for  use  with 
NACLGRAPH.FTN  program 

Appendix  E:  NACLGRAPH.FTN  Fortran  listing 

Appendix  F:  CALCOMP  generated  plot  (semi logarithmic)  of  the 
GRPHX.XXX  data. 

Appendix  G:  User's  guide  for  the  NACLRQFT.FTN  computer  program 

Appendix  H:  User's  guide  for  the  NACLGRAPH.FTN  computer  program 
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C  THIS  PROGRAM  CALCULATES  PROTECTION  FACTORS  FOR  DATA  COLLECTED  ON 
C  THE  USAFSAM/VNL  SODIUM  CHLORIDE  RQFT  INSTRUMENT. 

C 

C  CALCULATIONS  ARE  BASED  ON  INTEGRATING  THE  PHOTOMULTIPLIER 
C  TUBE'S  OUTPUT  OVER  A  VARIABLE  TIME  PERIOD  (10-20  SECONDS)  FOR 

C  EACH  EXERCISE  OF  THE  RQFT  PROT(XOL.  SINCE  THE  INTEGRATOR'S  j 

C  SENSITIVITY  IS  1000  COUNTS  PER  (VOLT-SECOND) ,  AND  THE  TIME  PERIOD 
C  FOR  EACH  EXERCISE  IS  VARIABLE  (10-20  SECONDS),  A  SCALED  OR 

C  'AVERAGE'  VOLTAGE  FOR  THE  PHOTOMULTIPLIER  OUTPUT  CAN  BE 

C  CALCULATED  BY  THE  FOLLOWING  RELATIONSHIP; 

C  I 

C  APTVO=(IC)*(l/IS)*(l/VTPEE) 

^  i 

C  WHERE,  APTVO=AVFRAGE  PHOTOMULTIPLIER  TUBE  VOLTAGE  OUTPUT  1 

C 

C  IC=INTEGRATOR  COUNT 

C 

C  IS=INTEGRATOR  SENSITIVITY 

C 

C  VTPEE=VARIABLE  TIME  PERICD  FOR  EACH  EXERCISE 

C 

C  NEXT,  A  LEAST-SQUARES  CURVE  FITTING  METHOD  IS  USED  TO  GENERATE  A 
C  POLYNOMIAL  FUNCTION  OF  THE  FORM:  Y»F(X).  THE  DATA  USED  TO 
C  GENERATE  THIS  FUNCTION  ARE  THE  ORDERED  PAIRS  OF  VALUES  OF  THE  FORM 
C  (X,Y),  TABULATED  WHEN  CALIBRATING  THE  RQFT  INSTRUMENT  WITH  THE 

C  SODIUM  CHLORIDE  CALIBRATION  STANDARDS.  IN  THIS  CASE; 

C 

C  X=PHOTOMOLTIPLIER  DETECTOR  TUBE  VOLTAGE  OUTPUT  (IN  VOLTS) 

C 

C  Y-SODIUM  CHLORIDE  CAIIBRATION  STANDARD  CONCENTRATION 

C 

C  THE  RESULT  OP  THE  CURVE  FITTING  PROCEDURE  WILL  BE  THE  CAPABILITY 

C  TO  DIRECTLY  CALCULATE  THE  MASK  LEAK  CONCENTRATION  FOR  AN  EXERCISE. 

C  THIS  IS  SO  BECAUSE  THE  POLYNOMIAL  CURVE  PITTING  RELATIONSHIP  OP 
C  THE  FORM;  Y-P(X) ,  DIRECTLY  RELATES  THE  PHOTOMULTIPLIER  TUBE'S 
C  AVERAGE  OUTPUT  VOLTAGE  TO  CONCENTRATION  (THAT  IS,  THE  INTEGRATED 
C  COUNT  IS  DIRECTLY  RELATQ)  TO  MASK  LEAKAGE  CONCENTRATION  FOR  A 

C  PARTICULAR  EXERCISE) .  j 

C  i 

C  THUS,  FOR  A  PARTICULAR  EXERCISE,  A  MASK  LEAKAGE  CONCENTRATION 

C  IS  CALCULATED  BY  SUBSTITUTING  THE  'ADJUSTED'  INTEGRATOR  COUNT,  ; 

C  (THAT  IS,  THE  APTVO  VALUE),  INTO  THE  POLYNOMIAL  CURVE  FITTING  I 

C  FUNCTION.  ' 

C  i 

C  FINALLY,  A  MASK  LEAK  CONCENTRATION  IS  CONVERTED  TO  A  PROTECTION 
C  FACTOR  BY  THE  POLUWING  REUTIONSHI  Pj 

C  j 

C  PF-(CC)/(HL) 

C  WHERE,  PF-PROTECTION  FACTOR  | 

C  I 

C  CO-CHALLENGE  CONCENTRATION  (FOR  THIS  SYSTEM,  10  TO  THE  | 

C  ZERO  OR  1.0)  I 

C  I 

C  NL-KASK  LEAK  CONCENTRATION  FOR  A  PARTICUUR  EXERCISE  | 

C 


I 

j 
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C  SINCE  THE  LEAST-SQUARES  CURVE  FITTING  PROCEDURE  IS  THE  KEY  TO 

C  THIS  METHOD  OF  PROCESSING  THE  RQFT  DATA,  THE  FOLLOWING  COMMENTS 

C  WILL  DESCRIBE  ITS  IMPLEMENTATION  IN  THE  PROGRAM. 

C 

C  THIS  IS  A  PROGRAM  WHICH  UTILIZES  A  LEAST-SQUARES  METHOD  OF 
C  APPROXIMATION  TO  DO  CURVE  FITTING. 

C 

C  THIS  PROGRAM  WILL  TAKE  A  SET  OF  (N)  NUMBER  OF  DATA  POINTS  AND 

C  WILL  FIT  A  POLYNOMIAL  fUP  TO  DEGREE  (N-1) ]  TO  THE  SET  OF  POINTS. 

C 

C  THIS  PROGRAM  DETERMINES  THE  COEFFICIENTS  OF  THE  POLYNOMIAL  BY 
C  SOLVING  A  SYSTEM  OF  NORMAL  EQUATIONS  WHICH  ARE  DERIVED  BY 

C  CONSIDERING  SOME  QUANTITY  (Y)  AS  A  POLYNOMIAL  FUNCTION  OF  ANOTHER 

C  QUANTITY  (X). 

C 

C  TO  FIND  THE  BEST  FITTING  CURVE  FOR  A  GIVEN  SET  OF  DATA,  THE 

C  PROGRAM  USER  IS  ONLY  REQUIRED  TO  ENTER  THE  OBSERVED  DATA  POINTS 

C  AMD  THE  DEGREE  OF  THE  POLYNOMIAL  APPROXIMATION. 

C 

C  THE  FIRST  PART  OF  THE  PROGRAM  (STATEMENT  LABEL  NUMBERS:  5101  THRU 
C  403)  ACCEPTS  THE  DATA,  STORES  IT,  AND  ORGANIZES  IT  FOR  USE  IN 
C  THE  REMAINDER  OF  THE  PROGRAM. 

C 

C  THE  SECOND  PART  OF  THE  PROGRAM  (STATEMENT  LABEL  NUMBERS:  400  THRU 

C  500)  ESTABLISH  THE  NORMAL  EQUATIONS.  STATEMENT  LABEL  NUMBERS:  30 

C  THRU  90  SOLVES  THE  NORMAL  EQUATION  SYSTEM.  IN  THE  FINAL  SECTION 

C  OF  THE  PROGRAM  (STATEMENT  NUMBERS:  16  THRU  11)  ,  THE  ERRORS 

C  (DIFFERENCES  BETWEEN  THE  ACTUAL  AND  THE  ESTIMATED  VALUES)  ARE 
C  CALCULATED.  IN  ADDITION  THE  SUM  OF  THE  SQUARE  ERRORS  IS 
C  CALCULATED.  THE  POLYNOMIAL  FUNCTION  PRODUCING  THE  SMALLEST 

C  SUM  OF  SQUARE  ERRORS  YIELDS  THE  BEST  CURVE  FIT.  OBSERVATION 
C  OF  THE  INDIVIDUAL  ERRORS,  ON  THE  OTHER  HAND,  GIVES  A  CLUE 

C  AS  TO  THEIR  RELATIVE  DISTRIBUTION  ABOUT  THE  BEST  FIT 

C  POLYNOMIAL  CURVE. 

C 

C  THE  USER  IS  GIVEN  THE  OPTION  OF  SELECTING  A  FUNCTIONAL 
C  DEFINITION  OF  THE  VARIABLE  (X)  SO  THAT  APPROXIMATIONS 
C  CAN  BE  HADE  TO  CURVES  THAT  ARE  NOT  NECESSARILY 

C  POLYNOMIALS  IN  (X).  FOR  EXAMPLE:  Y»EXP(X)  OR  Y-l/ (X) . 

C 

***************  ********, 


Q*****  ***** 

Q*****  IF  YOU  HAVE  ANY  QUESTIONS  CONCERNING  THIS  PROGRAM  CALL  *’“** 

***** 

C*****  CAPTAIN  EDWARD  S.  KOLESAR,  JR.  ***** 

C*****  ***** 

C*****  USAFSAM/VNL  BROOKS  AFB,  TX  ***** 

C*****  ***** 

*****  AUTOVON  240-2154  COMMERCIAL  (512)536-2154  ***** 

*****  ***** 


ARRAYS  ARE  DOUBLE  PRECISION  TO  INSURE  ACCURACY  OP  THE 
CALCULATIONS. 

IMPLICIT  INTEGBR*4  (I-N) 
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INTEGER  FILTK 

DOUBLE  PRECISION  D( 8, 9) ,W( 0 ; 13) , 2 ( 0 : 13 ) ,B ( 8, 9) ,E ( 8, 9) 

DOUBLE  PRECISION  XAS(7) ,A(14) ,X<13) ,Y(13) .THEORY , RESID,SRESD, XXX 
DOUBLE  PRECISION  DEL.XMIN.XMAX ,XEST ,BEGIN,XXMIN,XXMAS,XXAS(7 ) 
DOUBLE  PRECISION  PF (17) ,IC (17) , ICTP(17) ,XXXX 
DIMENSION  GRPH1(14) ,GRPH2(802) ,ICCBS(2) 

DIMENSION  XX(401) ,YX(401) 

BYTE  SECNO) 

BYTE  P(9)  ,C(10)  .GdO) 

BYTE  ANS.COMMEN.REP.COM 
BYTE  SELECT  (7)  .GROUP  1  .GROUP  2 
BYTE  YES. NO 

BYTE  NAME(45) .MASK(45) ,DATE(45) .TIHE(45) 

DATA  SECN/'1'.'2','3’. U'.'S', '6', d*. '8' . '9'/ 

YES»'Y' 

NO='N' 


ESTABLISH  A  FILE  COOMTEE  AND  DECLARE  THE  FILE  NAMES. 
5101  FILTK-1 

VARIABLES  ARE  SET  EQUAL  TO  ZERO  SO  THAT  ITERATIVE  RUNS 
OH  THE  SAME  DATA  CAN  BE  READILY  PERFORMQ). 


H-0 

NP-0 

ICCBS(l)«i72*256+27 
ICCBS (2) -74*256+27 
4137  FORMATdH  ,2A2) 
ICl-0.0 
1C2-0.0 
JC3-0.0 
IC4-0.0 
lCS-0.0 
IC6-0.0 
IC7-0.0 
IC8-0.0 
IC9-0.0 
leio-o.tS) 
tcu-o.o 

ICl 2-0.0 
IC13-0.0 
IC14-0.0 
IC15-0.0 
1C16-0.0 

xci-o.e 

XC2>0,0 

XC3-0.0 

XC4-0.0 

XCS-0.0 

XC6-0.0 

XC7-0.0 

XC8-0.0 

XC9-0.0 

XClO-O.tt 

xcn-0.0 

XC12-C.0 
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XC13-0.0 

XC14»0.0 

XC15«>0.0 

XC16-O.0 

CALCPF=0.0 

XEST-0.0 

PFEST-0.0 

DEL=0.0 

THEORY=0.0 

RESID-0.0 

SRESD-0.0 

XXX-0.0 

BEGIN-0.0 

DO  9333  1-1,17 

IC(U-0.0 

PF(I)»0.0 

ICTP(I)-0.0 

9333  COWTINUE 

DO  4153  1-1.7 
XAS(I)»0.0 
XXAStn-O.O 
4153  CONTINUE 

THE  FILES  THAT  HOLD  VARIOUS  SECMENTS  OF  DATA  ARE  NAMED. 

DO  9334  1-1,13 
X(I)-0.0 
Ydl-O.O 

9334  CONTINUE 
P<l)-'D’ 

Pt2)-'A' 

P(3) -'T' 

P(4I-*A' 

PtS)-' 

Pt9l-0 

eti)«'c' 

et2)-'A' 

C(3!-'L’ 

Ct4)-'C' 

CtSI-SBCNIFlLTKI 
C«lQ!»0 
G(2)»'R' 

Ot4)-'tt‘ 

C{51«SECN(F1LTX» 

. ' 

OdOI-O 

THE  FILE  CALLED  DATA. XXX,  CONTAINS  THE  SODIUM  C1IIA)RIDE 
CALinRATION  STAKOARD  CONCENTRATION  DATA  TO  d£  PITTED 
WSTM  THE  WLVNOMIAL  FUNCTION. 

THE  PILE  CALtfiD  CALCX.XXX,  CONTAINS  TNE  CALCOUTED  WLVNOMIAL 
COEFFICieiiTS,  RESIDUALS,  AND  OTUER  DCSCRmiVE  INFORMAtlON. 
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THE  FILE  CALLED  GRPHX.XXX  CONTAINS  THE  FOLLOWING  ARRAYS: 

(GRPHl]  SODIUM  CHLORIDE  CALIBRATION  STAJJDARD 
CONCENTRATION  DATA. 

[GRPH2)  ESTIt4ATED  VALUES  DERIVED  FROM  EXERCISING 
THE  POLYNOMIAL  FUNCTION. 

IF  A  PLOT  OF  THE  ACTUAL  DATA  AND  ESTIHATO)  VALUES  IS  DESIRED, 
THE  USER  CAN  USE  THE  FILE  CALLED  GRPHX.XXX  FOR  THIS  PURPOSE. 


NUMBER  THE  FILES  SEQUENTIALLY  SO  THAT  THEY  CAN  BE  EASILY  RETRIEVED 
FOP  PRINTING  AND  PLOTTING. 

TYPE  2006 

2006  FORMATClX,'  •///) 

1002  FORMAT  (I  All 

PROGRAM  USER  WILL  UNIQUELY  IDEITTIFY  EACH  DATA  FILE  WITH  A 
SEQUENTIAL  NUMBERING  SYSTEM  FOR  EASE  OF  RECALL. 

TYPE  4137,ICC8S(11 .ICCBS(2) 

TYPE  2006 
TYPE  2002 

2002  FORHATIIX, 'USER  ATTENTION:  IN  ORDER  TO  KEEP  TRACK  OF  THE  DATA'/ 

C  SETS  BEING  ANALYZED,  IT  IS  RECOMMENDED  THAT  THEY  BE  '/ 

C  SEQUENTIALLY  NUMBERED. •/> 

TYPE  306 
TYPE  2003 

2003  FORMAT (IX, 'ENTER  THE  FQLUOMINGt  001  FOR  THE  FIRST  DATA  SETj  '/ 

C  002  FOR  THE  SECOND  DATA  SBTj  003  FOR  THE  THIRD  DATA  SET.  ETC. '/I 
TYPE  306 
TYPE  2004 

2004  FORMAT! IX, ‘ ENTRY*  ',SJ 
ACCEPT  200S,P(6» ,P(7) ,Pt$) 

2006  FORHATi3Alt 
TYPE  2006 

TYPE  4n7,tceas(n,!ce6st2» 

TYPE  2006 

typf;  6100 

6100  FO»(AttlX, 'USER  ATTEWTIONi  IN  ORDER  TO  .KESP  TRACK  OF  THE'/ 

C  POLYNOMIAL  CURVE  FITTING  COWFICIEOTS  AHC  SBSIOUALS,  IT'/ 

C  58  HECOMKENOeO  THAT  THEY  Ot  S6<J«SMtIALLY  NUMSERBD.*/) 

TYPE  306 
TYM;  6200 

6200  rOW(ATilX.'Kf{TE»  THE  FOLLOWING:  001  FOR  THE  FIRST  JiEStOUAL  S£T{'/ 
C  002  I'OR  THE  SECOND  RESIDUAL  SET:  003  FOR  TUB  THIRD,  ETC.'/) 

TYPE  306 
*n*PE  6300 

6300  FOWATt  IX, 'ENTRY-  ',tj 

ACCEPT  640Q,C(7» ,Ct6),C<9> 

6400  FORMATtSAll 
TYPE  2006 

TYPE  413t,tCC&3(H ,ICCOS(2) 

TYPE  2006 
TYPE  4001 

4001  FORMAT!  IX, 'USER  ATTEWriONs  IN  ORDER  TO  KEEP  TRACK  OF  THE  GRAW/ 
C  sets  BEING  AtlALYSED,  IT  IS  RECOMNEKDEO  THAT  TtlEY'  OE'/ 
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C  SEQUENTIALLY  NUMBERED*/) 

TYPE  306 
TYPE  4004 

1004  FORMAT< IX, 'ENTER  THE  FOLLOWING:  001  FOR  THE  FIRST  GRAPH  SET;'/ 
C  002  FOR  THE  SECOND  GRAPH  SET;  003  FOR  THE  THIRD,  ETC.'/) 

TYPE  306 
TYPE  2004 

ACCEPT  2005,G(7) ,G(8) ,G:9) 

TYPE  2006 

TYPE  4137,ICC8S(1) ,ICCBE(2) 

TYPE  2006 
REP»'N' 

GO  TO  3081 
6000  FILTK=FILTK+l 

C(5) =SECN(FILTK> 

G(5)*SECN(FILTK) 

CLOSE  (UNIT- 3) 

CLOSE ( UNIT- 1) 

OPEN  ( UNI T» 3  ,NAKE«>C, DISPOSE- ' SAVE' , TYPE  - ' NEW' ) 

OPEN ( UNI T- 1 , NANE-G, DISPOSE- *  SAVE  * .TYPE- ' NEW • ) 

REWIND  2 
GO  TO  6001 

3081  OPEN(UNlT-2.NAME-P,DlSPOSE-*SAVE' .TYPE-'NEW') 

OPEN ( UNIT- 3 , NAMB-C, DI SPOSE- *  SAVE '.TYPE-'NEW') 

OPEN ( UNIT- 1 . NAME-G, DISPOSE- ' SAVE ' ,TY PE- *  NEW • ) 

C 

C  ENTER  THE  DATA  TO  BE  PROCESSH).  IN  THIS  SITUATION.  THE  VOLTAGE 

C  NEASURENEWrS  ARE  THE  X-AXIS  DATA  AND  THE  SODIUM  CHLORIDE 

C  STANDARD  CALIBRATION  CONCEOTRATIONS  THE  Y-AXIS  DATA. 

C 

1003  TYPE  1000 

lUOO  FORMAT (IX, 'ENTER  THE  NUMBER  OP  SODIUM  CHLORIDE  CALIBRATION'/ 

C‘  CONCENTRATION  STANDARDS'./) 

TYPE  2004 
ACCEPT  2000..NP 
2000  PORHAT(H2) 

C 

C  OUTPUT  THE  JNFOBMATIO.N  OSH)  TO  PERFORM  THE  CALCULATIONS  IN 
C  THE  OATAX.XXX  FILE. 

C 

TYPE  2Q06 

TYPE  *ni,icc»s(i>,icctts(a» 

T¥I>E  2006 
WHtTKt2,2006) 

200?  FORMATtiX. //////) 

K«lTE(2.30ai)NP 

3001  F'OI!MAttlX.24X,'THE  NUMBER  OF  SODIUM  CHLORIDE 
C  CALIBRATION  CONCEMTHATION  STANDARDS-  *,12) 

WRtTfil2,200«) 

TYPE  1000 

1000  FORMATt  lx, 'DEPRESS  RETURN  XEV  AFTER  EErTERlNC  VOLTAGE*/ 

C*  MEASUREMENT'//) 

TYPE  2006 

C 

C  THE  FOLUMtKG  INFORMATION  IS  AVAILAOU  FROM  THE  DATA  SHEET 
C  USED  OURIMG  AH  ROFT  EVALUATION. 

C 
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TYPE  4000 

4000  FORMAT (IX, 'ENTER  THE  DATA  POINTS'/ 

C'  SAMPLE  NUMBER  VOLTAGE  SODIUM  CHLORIDE  / 

C  MEASUREMENT  CALIBRATION'/ 

IN  VOLTS  CONCENTRATION'/ 

C  (X  DATA)  (I  DATA)  '/) 

WRITE(2,3003) 

3003  FORM AT C IX, 2 4X, 'SAMPLE  NUMBER' ,18X, 'VOLTAGE ' ,23X, 'SODIUM  CHLORIDE') 
WRITE(2,3082) 

3  082  FORMAT  C IX , 2 4X , 1 3X , 18X , ' MEASUREMENT ' , 19 X , ' CALIBRATION  ) 
WRITE(2,3083) 

3083  F0RMAT(1X,24X,13X,18X, ' (IN  VOLTS) ', 20X , 'CONCENTRATION ' ) 
WRITE(2,3C84) 

3084  F0RMAT{1X,24X,13X,18X, * (  X  DATA  )',20X.'(  Y  DATA  )  /) 

LINES=0 

DO  5001  1=1, NP 
TYPE  4500,  I 

4500  FORHAT(lH+,T7 ,1I2,T20,1H  ,S) 

ACCEPT-5000,  X(I) 

CALL  CLEAR(LINES) 

TYPE  4501,  I,X(I) 

4501  FORMAT( '+ ' ,T7,1I2,T21,1G13.6, '  ',?) 

ACCEPT  5000,  Y(I) 

TYPE  305 

5000  F0BHAT(1G13.6) 

XAS(I)=X(I) 

XXAS ( 1 ) =X (I ) 

WRITE(2, 5002)1, X(I) ,V(I) 

5002  F0RHAT{1X,29X,1I2,23X,1G13.6,19X,1G13 .6) 

5001  CONTINUE 

C  ENTER  THE  DESCRIPTIVE  INFORMATION  CONCERNING  THE  SUBJECT,  MASK, 

C  DATE,  AND  TIME  TESTED. 

C 

TYPE  2006 

TYPE  4137,ICCBS(1)  ,ICCBS(2) 

TYPE  2006 
3  077  TYPE  3  085 

3085  FORMAT (IX, 'SUBJECT  NAME!',2X,5) 

ACCEPT  3 08 6  , NAME 

3086  FORMAT  ( 4  5A1) 

TYPE  3087 

3087  FORMAT (IX, 'TYPE  OF  MASK:',2X,S) 

ACCEPT  3086  , MASK 

TYPE  3088 

3088  FORMATdX, 'DATE  TESTED  t ',2X,  5) 

ACCEPT  3089, DATE 

3089  FORMAT  ( 4 5A1) 

TYPE  3  090 

3090  FORMAT  (IX, 'TIME  TESTED  i ', 2X ,  $) 

ACCEPT  3089, TIME 

TYPE  2006 
3092  F0RMAT(I2) 

TYPE  2006 

C  ENTER  ■''HP  INTEGRATION  EXERCISE  COUNT  DATA  FOR  THE  TEST  PROTOCOL. 

C 
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C  THIS  INFORMATION  IS  AVAILABLE  FROM  THE  DATA  SHEET 
C  USED  DURING  AN  RQFT  EVALUATION. 

C 

C 

C  SELECT  THE  EXERCISE  PROTOCOL. 

C 

TYPE  2006 

TYPE  4137,ICCBS(1) ,ICCBS(2) 

TYPE  2006 
4139  TYPE  3127 

3327  FORMATdX, 'THE  USER  IS  FREE  TO  SELECT  ONE  OF  TWO  GROUPS  OF’/ 

C  EXERCISE  PROTOCOLS.'//) 

TYPE  3128 

3128  FORMATdX, 'THE  tGROUP  1)  EXERCISE  PROTOCOL  CONSISTS  OFi '// 

C  111  NORMAL  BREATHING  STRAIGHT  AHEAD'/ 

C  (21  DEEP  BREATHING  STRAIGHT  AHEAD'/ 

C  (31  TALKING’/ 

C  (4)  SIDE-TO-SIDE  HEAD  MOVEMENTS  (DEEP  BREATHING)'/ 

C  (51  UP-AND-DOWN  HEAD  MOVEMENTS  (DEEP  BREATHING) '/ 

C  (61  FACIAL  GRIMACING'///) 

TYPE  2006 

TYPE  4137,ICCBS(1) ,ICCBS(2) 

TYPE  2006 
TYPE  3129 

3129  FORMATdX, 'THE  (GROUP  21  EXERCISE  PROTOCOL  CONSISTS  OFi'// 

C  (11  NORMAL  BREATHING  STRAIGHT  AHEAD'/ 

C  (21  NORMAL  BREATHING  LEFT'/ 

C  (31  NORMAL  BREATHING  RIGHT'/ 

C  (41  NORMAL  BREATHING  DOWN'/ 

C  (51  NORMAL  BREATHING  UP'/ 

C  (61  DEEP  BREATHING  STRAIGHT  AHEAD'/ 

C‘  (71  DEEP  BREATHING  LEFT'/ 

C  (81  DEEP  BREATHING  RIGHT') 

TYPE  3130 

3130  FORMATdX, '(91  DEEP  BREATHING  DOWN'/ 

C  (101  DEEP  BREATHING  OP'/ 

C  (111  TALKING'/ 

C  (121  FACIAL  GRIMACING'/ 

C  (131  SIDE-TO-SIDE  HEAD  MOVEMENTS  (NORMAL  BREATHING)'/ 

C*  (141  UP-AND-DOWN  HEAD  MOVEMENTS  (NORMAI  BREATHING)'/ 

C  (151  SIDE-TO-SIDE  HE,')D  MOVEMENTS  (DEEP  BREATHING)'/ 

C  (161  UP-AND-D(»»N  HEAD  MOVEMENTS  (DEEP  BREATHING) '//) 

4138  TYPE  2006 

TYPE  4137,ICCDS(1) ,ICCBS(2) 

TYPE  2006 
TYPE  3131 

3131  FORMATdX, 'TO  SPECIFY  THF,  EXERCI.SE  PROTWOL  GROUP  OP  INTEREST,'/ 

C  TYPE  eiTHER;  GROUP  1  OR  GROUP  2  ') 

TYPE  306 
TYPE  3199 

3199  FORMATdX, 'ENTRY  -  ',S) 

ACCEPT  3122, SELECT 
3122  FO.RHATt7Al) 

If  (SBLBCTi'7)  .NE. '1*.AND.SBLECT{7)  .NS.  *2')  GOTO  4138 
TYPE  2006 

TYPE  4137,1CCDS(1) ,ICCPS(2) 
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TYPE  2006 
TYPE  3021 

3021  FORMATdX, 'USER  ATTENTTON;  IF  NO  EXERCISE  COUNT  DATA  WAS'/ 

C  COLLECTED  FOR  A  PARTICUIAR  EXERCISE,  TYPE:  000001.  ALSO,'/ 

C  FOR  EACH  TYPED  EXERCISE  COUNT  DATA  ENTRY,  SIX  DIGITS  MUST') 

TYPE  3022 

3022  FORMATdX, 'BE  TYPED,  THAT  IS,  IF  YOU  HAVE  A  SIX  DIGIT  NUMBER,'/ 

C  TYPE  ALL  SIX  DIGITS.  IF  YOU  HAVE  A  FIVE  DIGIT  NUMBER,  TYPE'/ 

C  ONE  LEADING  ZERO  AND  THEN  THE  FIVE  DIGITS.  IF  YOU  HAVE  A') 

TYPE  3033 

3033  FORMATC  FOUR  DIGIT  NUMBER,  T^PE  TWO  LEADING  ZEROS  AND  THEN  THE'/ 
C  FOUR  DIGITS,  ETC.  SEVERAL  EXAMPLES  FOLLOW  AS  AN  ILLUSTRATION') 

TYPE  3034 

3034  FORMATdX, 'FOR  EXAMPLE:  COUNT  DATA=743 182  TYPED  ENTRY-743182') 

TYPE  3035 

3035  FORMATdX, 'FOR  EXAMPLE:  COUNT  DATA-18726  TYPED  ENTRY-018726') 
TYPE  3036 

3036  FORMATdX, 'FOR  EXAMPLE:  COUNT  DATA-6412  TYPED  ENTRy-006412' ) 
TYPE  2006 

TYPE  4137,ICCBSd)  ,ICCBS(2) 

TYPE  2006 
TYPE  3119 

3119  FORMATdX, 'DEPRESS  THE  RETORN  KEY  AFTER  ENTERING  AN'/ 

C  INTEGRATOR  COUNT'/) 

TYPE  2006 
TYPE  3037 

3037  FORMATdX, 'EXERCISE  COUNT  DATA:'/) 

TYPE  3038 

3038  FORMATdX. 'EXERCISE', 29X,' INTEGRATOR', 7X, 'TIME  PERIOD') 

TYPE  3126 

3126  FORMATdX, 40X, 'COUNT' ,9X, '  (IN  SECONDS)'/) 

IF(SELECT(7)  .EQ.  '2')  GO  TO  3133 
TYPE  3039 

3039  FORMATdX, 'NORMA!,  BREATHING  STRAIGHT  AHEAD  ',$) 

ACCEPT  3 040, Id 

CALL  CLEAR (LINES) 

3040  FORItAT(I6) 

TYPE  3112, ICl 

3112  FORMAT(lH+, 'NORMAL  BREATHING  STRAIGHT  AHEAD ' ,T39 , 16 ,T59, IH  ,*) 
ACCEPT  3113,ICTP1 

3113  KORMAT(i2) 

TYPE  3041 

3941  FORMATdX, 'DEIT  BREATHING  STRAIGHT  AHEAD  ',*) 

ACCEPT  3040,IC2 
CALL  CLEAR (LIMBS) 

TYPE  3114, IC2 

3114  FORMAT  dll*, 'DEEP  BREATHING  STRAIGHT  AHEAD' ,T39,I6,TS9, 

•lit  ,S) 

ACCEPT  3113,ICTP2 
TYPE  3042 

3042  FORMATdX, 'TALKING  *,*) 

ACCEPT  3040,IC3 
CALL  CLEAR (LINES) 

TYPE  II1S,1C3 

3US  FOHMATdH-»,'TALKlKG',T39,16,T59,lH  ,  *) 

ACCEPT  3in,lCTP3 
TYPE  3043 
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3043  FORMAT ( IX, 'SIDE-TO-SIDE  HEAD  MOVEMENTS'/ 

C  (DEEP  BREATHING)  ',2 OX, IH  ,  S> 

ACCEPT  3040,IC4 
CALL  CLEAR (LINES) 

TYPE  3116, IC4 

3116  FORMAT  ( 1H+,  '  (DEEP  BREATHING)  '  ,T39 , 16  ,T59 ,  IH  ,  $) 

ACCEPT  3113,ICTP4 

TYPE  3044 

3044  FORMAT (IX, 'UP-AND-DOWN  HEAD  MOVEMENTS'/ 

C  (DEEP  BREATHING)  ',20X,1H  ,$) 

ACCEPT  3 040, ICS 
CALL  CLEAR (LINES) 

TYPE  3116, ICS 
ACCEPT  3113,ICTP5 
TYPE  3045 

3045  FORMATdX, 'FACIAL  GRIMACING ', 20X , IH  ,S) 

ACCEPT  3Q40,IC6 

CALL  CLEAR (LINES) 

TYPE  3117, IC6 

3117  FORMAT (1H+, 'FACIAL  GRIMACING' ,T39, 16, T59,1H  , S) 

ACCEPT  3113,ICTP6 

TYPE  2006 

TYPE  41  37,ICCBS(1)  ,ICCBS(2) 

TYPE  2006 
GO  TO  3153 

3133  CONTINUE 
TYPE  3039 
ACCEPT  3040, ICl 
CALL  CLEAR(LINES) 

TYPE  3112, ICi 
ACCEPT  3113,tCTPl 
TYPE  3134 

3134  FORMATdX, 'NORMAL  BREATHING  LEFT  ,S) 

ACCEPT  3040,IC2 

CALL  CLEAR (LINES) 

TYPE  3135, IC2 

3135  FORMAT (IH*, 'NORMAL  BREATHING  LEFT ‘ ,T39 , 16 ,T59, IH  ,S) 
ACCEPT  3113,IC'TP2 

TYPE  3136 

3136  FORMATdX, 'NORMAL  BREATHING  RIGHT  ',S) 

ACCEPT  3040, IC3 

CALL  CLEAR(LINES) 

TYPE  3137,IC3 

3137  FORMAT dll*, 'NORMAL  BREATHING  RIGHT* ,T39, 16, T59, IH  ,5) 
ACCEPT  3U3,1CTPI 

TYPE  3138 

3138  FORMATdX, 'NORMAL  BREATHING  DOWN  ',$) 

ACCEPT  3040,1C4 

CALL  CLEAR (LINES! 

TYPE  3139, IC4 

3139  TORMATdH*, 'NORMAL  BREATHING  DOWN' ,T39 , 16 ,T59, IH  ,S) 
,\CCEPT  31i3,ICTP4 

TYI’E  3140 

3140  FORMATdX, ‘NORMAL  BREATHING  OP  ',$) 

ACCE1*T  1040,  ICS 

CALL  CLEAR (LINES) 

TYPE  3 141, ICS 
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3141  FORMAT (1H+, 'NORMAL  BREATHING  OP ’ ,T39 , 16 ,T59, IH  , $) 
ACCEPT  3113,ICT?5 

TYPE  3  041 
ACCEPT  3040, IC6 
CALL  CLEAR (LINES) 

TYPE  3114,  IC6 
ACCEPT  3113,ICTP6 
TYPE  3142 

3142  FORMAT (lx, 'DEEP  BREATHING  LEFT  ',$) 

ACCEPT  3040, IC7 

CALL  CLEAR (LINES) 

TYPE  3143,IC7 

3143  F0RMAT(1H+,'DEEP  BREATHING  LEFT  * ,T39 , 16 ,T59 , IH  , $) 
ACCEPT  3113,ICTP7 

TYPE  3144 

3144  FORMATdX, 'DEEP  BREATHING  RIGHT  ',$) 

ACCEPT  3040,IC8 

CALL  CLEAR (LINES) 

TYPE  3 14 5,  ICS 

3145  FORMAT(lH+, 'DEEP  BREATHING  RIGHT ' ,T39 , 16 ,T59, IH  ,$) 
ACCEPT  3113,ICTP8 

TYPE  3146 

3146  FORMATdX. 'DEEP  BREATHING  DOWN  ',$) 

ACCEPT  3040,IC9 

CALL  CLEAR(LINES) 

TYPE  3147,  IC9 

3147  FORMATdH+,'DEEP  BREATHING  DOWN*  ,T39 , 16  ,T59 ,  IH  ,  $) 
ACCEPT  3113,ICTP9 

TYPE  3148 

3148  FORMATdX, 'DEEP  BREATHING  UP  ',$) 

ACCEPT  3040, ICIO 

CALL  CLEAR(LINES) 

TYPE  3149  ICIO 

3149  FORMATdH+, 'DEEP  BREATHING  UP' ,T39 , 16 ,T59 , IH  ,S) 
ACCEPT  3113,ICTP10 

▼yPE  3042 
ACCEPT  3040, ICII 
CALL  CLEAR(LINES) 

TYPE  3115,1011 
ACCEPT  3113,ICTP11 
TYPE  3045 
ACCEPT  3040, 1C12 
CALL  CLEAR (LINES) 

TYPE  3117, IC12 
ACCEPT  3113,ICTP12 
TYPE  3150 

3150  FORMATdX, 'SIDE-TO-SIOE  HEAD  MOVEMENTS'/ 

C  (NORMAL  BREATHING)  ',*) 

ACCEPT  3040, ICl 3 
CALL  CLEAR (LINES) 

TYPE  3151, IC13 

3151  FORMATdll+,' (NORMAL  BREATHING)  '  ,T39 , 16 ,T59,  IH  ,5) 
ACCEPT  3113,1CTPU 

TYPE  3152 

3152  FORMATdX, 'UP-AND-DOWN  MEAD  MOVEMENTS'/ 

C  (NORMAL  BREATHING)  ',$) 

ACCEPT  3040, IC14 
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CALL  CLEAR (LINES) 

TYPE  3151, IC14 
ACCEPT  3113,ICTP14 
TYPE  3043 
ACCEPT  3040, IC15 
CALL  CLEAR (LINES) 

TYPE  3116,  IC15 
ACCEPT  3113,ICTP15 
TYPE  3044 
ACCEPT  3040, IC16 
CALL  CLEAR (LINES) 

TYPE  3116, IC16 
ACCEPT  3113,ICTP16 
TYPE  2006 

TYPE  4137,ICCBS(1) ,ICCBS(2) 

3153  TYPE  2006 

WRITE(2,2006) 

WR1TE(2, 3046) NAME 

3046  POBMAT(6X, 'SUBJECT  NAMEs * ,2X,45A1) 

WRITE (2, 3 047) MASK 

3047  FORMAT (6X, 'TYPE  OF  MASK: ' ,2X,45A1) 

WRITE(2,3048)DATE 

3048  FORMAT(6X, 'DATE  TESTED  :',2X,45A1) 

WRITE(2,3049)TIME 

3049  FORM AT (6X, 'TIME  TESTS)  :',2X,4SA1) 

WRITE(2,2006) 

WRITE(2,9006) 

9006  PORMAT(6X, 'BXERCISF.  INTEGRATOR  COUNT  DATA:'/) 

WRITE(2,3094) 

3094  FORMAT(6X, ' EXERCISE', 29X,* INTEGRATOR  COUNT', 7X, 'TIME  PERIOD') 
WRITE(2,3120) 

3120  PORMAT(6X,8X,29X,16X,7X, ' (IN  SECONDS)',/) 

IF(SELECT(7)  .EQ.  '2')  GO  TO  3154 

XCl-ICl 

XC2-IC2 

XC3-IC3 

XC4-IC4 

XC5-IC5 

XC6-IC6 

IC(1)»XC1 

IC(2)»XC2 

IC(3)«XC3 

IC(4)-XC4 

IC(5)«XC5 

IC(6)-XC6 

lCTP(l)«ICTrl 

ICTP(2)«ICTP2 

ICTP(3)-ICTP3 

ICTP(4)-ICTP4 

ICTP(5)-rCTP5 

ICTP(6)-ICTP6 

WRITE(2,3051) ICl.lCTPl 

3051  F0RKAT(6X, 'NORMAL  BREATHING  STRAIGHT  AHEAD* ,1 IX. 16, 17X, 12) 
WRITB{2,3052)IC2,ICTP2 

WRITE(2, 3053) ICS, ICTP3 

3052  PORNAT(6X,'DEEP  BREATHING  STRAIGHT  AHEAD* ,13X,I6,17X,I2) 

3053  FORMAT(6X, 'TALKING', 35X, 16, 17X, 12) 
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WRITE<2,3054)IC4,ICTP4 

3054  FORMAT (6X, 'SIDE-TO-SIDE  HEAD  MOVEMENTS’/ 

C  (DEEP  BREATHING) • ,26X, 16, 17X, 12) 

WRITE (2, 3055) ICS, ICTP5 

3055  FORMAT (6X, 'UP-AND-DOWN  HEAD  MOVEMENTS'/ 

C  (DEEP  BREATHING) ’ ,26X, 16, 17X, 12) 

WRITE(2,3056) IC6,ICTP6 

3056  FOBMAT(6X, 'FACIAL  GRIMACING ’ ,26X, 16 ,17X , 12//) 
GO  TO  3155 

3154  CONTINUE 
XC1=IC1 
XC2-IC2 
XC3=IC3 
XC4=»IC4 
XC5=IC5 
XC6-IC6 
XC7-IC7 
XC8-IC8 
XC9-IC9 
XCIO-ICIO 
XCll-ICll 
XC12-IC12 
XC13-IC13 
XC14-IC14 
XC15-1C15 
XC16-IC16 
IC(1)*XC1 
1C(2)»XC2 
IC(3)»XC3 
IC(4)»XC4 
IC(5)«XC5 
IC(6)«XC6 
IC(7).XC7 
1C(6)«XC8 
IC(9)«XC9 
IC(10)-XC10 

ic(ii)»xcn 

IC(12)-XC12 

IC(13)«XCU 

IC(14)-XC14 

1C(15)-XC15 

IC(16)-XC16 

ICTPtD-ICTPl 

SCTP(2)»ICTP2 

ICTP(3)»ICTP3 

ICTP(4)«1CTP4 

1CTP(5)-ICTPS 

tCTP(6)-lCTP6 

ICTP(7)«ICTP7 

ICTP(8)»ICTP8 

ICTP(9)«1CTP9 

ICTPtlO)-ICTPlO 

tCTP(U)»ICTPU 

lCTP<l2)«tCTPl2 

ICTP(13)-ICTP13 

1CTPC14)«ICTPU 

ICTP(15)-ICTP1S 


—APPENDIX  A- 


ICTP(16)=ICTP16 
WRITE<2,3051)IC1,ICTP1 
WRITE(2,3156) IC2,ICTP2 

3156  FORMAT (6X, 'NORMAL  BREATHING  LEFT  * ,21X , 16 , 17X, 12) 

WRITE(2,3157) IC3,ICTP3 

3157  FORMAT(6X, 'NORMAL  BREATHING  RIGHT ' ,20X , 16 ,17X, 12) 

WRITE  {2,3158)IC4,ICrP4 

3158  FORMAT(6X, 'NORMAL  BREATHING  DOWN* ,21X, 16 , 17X , 12) 

3159  FORHAT(6X, 'NORMAL  BREATHING  UP' ,23X, I6,17X, 12) 

WRITE ( 2 , 3 159) IC5 , ICTP5 

WRITE(2,3052) IC6,ICTP6 
WPITE(2,3160)IC7,ICTP7 

3160  FORMAT(6X, 'DEEP  BREATHING  LEFT ' , 23X, 16 , 17X , 12) 
l'fRITE(2,3161)  IC8,ICTP0 

3161  FORMAT(6X, 'DEEP  BREATHING  RIGHT ' , 22X , 16 , 17X, 12) 
WR:tE(2,3162)IC9,ICTP9 

3162  F0RMAT(6X, 'DEEP  BREATHING  DOWN' ,23X, I6,17X, 12) 

WR1TE(2,3163) ICIO.ICTPIO 

3163  FORMAT(6X, 'DEEP  BREATHING  OP ' ,25X, I6 ,17X , 12) 

WRITE(2,3Q53) ICll.ICTPll 
WRITE(2,3164)IC12,ICTP12 

3164  F0RMAT(6X, 'FACIAL  GRIMACING' ,26X, 16, 37X, 12) 
WRITE(2,3165)IC13,ICTP13 

3165  FORMAT! 6X, 'SIDE-TO-SIDE  HEAD  MOVEMENTS'/ 

C  (NORMAL  BREATHING) •,24X, 16, 17X, 12) 

WRITE(2,3166)IC14,ICTP14 

3166  F0BMAT(6X, 'UP-AND-OOWN  HEAD  MOVEMENTS'/ 

C  (NORMAL  BREATHING) ',24X, 16, 17X, 12) 

WRITE(2,3054)IC1S,ICTP1S 
WBITB(2,3055)IC16,ICTP16 
3155  REP* ' N ' 

3888  IF(RBP.BQ.NO)  GO  TO  3999 
GO  TO  9983 

C 

C  CALCULATE  THE  AVERAGE  PHOTOMULTIPLIER  TUBE  OUTPUT  VOLTAGE  FOR 
C  EACH  EXERCISE 
C 

3999  CONTINUE 

iF(SBLBCT(7)  .EQ.  *2')  GO  TO  3167 
DO  3093  I*l,6 

IC(I)-IC(I)/(1000.0*ieTP(l)) 

3093  CONTINUE 
GO  TO  3168 

3167  CONTINUE 

DO  3169  1*1,16 
IC(1)*1C(I)/{1000.0»ICTP(I)) 

3169  CONTINUE 

3168  CONTINUE 
END  FILE  2 

9983  CONTINUE 
REWIND  2 

6001  ReAD(2,2006) 

READ (2, 3001 )NP 
3002  F0RMAT{1X,52X,12) 

C 

C  AFTER  ALL  PERTINENT  INIXIRMATION  IS  ENTERH)  INTO  THE  OATA.XSX  FILE, 

C  THE  FILE  IS  RtVOUND  AND  THE  POINTER  MOVED  TO  READ  THE  DATA  NEEDED 
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C  TO  CALCULATE  THE  LEAST-SQUARES  POLYNOMIAL  FUNCTION. 

C 

READ (2, 2006) 

READ(2,3003) 

READ(2,3082) 

READ(2,3083) 

READ(2,3084> 

C 

C  TO  HAKE  THE  LEAST-SQUARES  POLYNOMIAL  CURVE  FIT  ROUTINE  AS 

C  FLEXIBLE  AS  POSSIBLE,  THE  USER  CAN  SELECT  A  VARIETY  OF 

C  FUNCTIONAL  DEFINITIONS  FOR  THE  VARIABLE  (X). 

C 

2013  TYPE  2006 

TYPE  4137,ICCBS(1) ,ICCBS(2) 

TYPE  2006 
TYPE  2  009 

2009  FORMATdX, 'THE  USER  IS  FREE  TO  SELECT  ANY  ONE  OF  THE  FOLLOWING'/ 
C  FUNCTIONAL  DEFINITIONS  OF  THE  VARIABLE  (X)./) 

TYPE  306 
TYPE  2010 

2010  FORMATdX,  '  tlJ  (X)  ’/ 


c 

(21 

1/ (X)  '/ 

c 

(31 

EXPONENTIAL 

(X)  '/ 

c 

(41 

EXPONENTIAL 

(-X)  '/ 

c 

(51 

EXPONENTIAL 

d/x) '/ 

c 

(61 

EXPONENTIAL 

(-1/X)  '/ 

c 

(71 

NATURAL  LOG 

(X)  '/ 

c 

(81 

LOG  BASE  TEN 

1  (X)  '/ 

c 

(91 

SIN  (X)  ’/) 

TYPE  306 
TYPE  2  011 

2011  FORMATdX, 'TO  SPECIFY  A  FUNCTIONAL  DEFINITION  OF  (X),  SELECT*/ 
C  THE  CORRESPONDING  NUMBER  INSIDE  THE  BRACKETS'/) 

TYPE  306 
TYPE  2004 
ACCEPT  2012, lAX 

2012  FORHATdll) 

TYPE  2006 

TYPE  4137,ICCBSd')  ,ICCOS(2) 

TYPE  2006 

IFtlAX.CT.O. ANO.1AX.lt. 10)  GO  TO  2014 
GO  TO  2013 
C 

C  DOUBLE  PRECISION  CALCUIATIONS  ARE  HADE  TO  RESOLVE  PROBLEMS 
C  ENCOUNTERED  WHEN  DEALING  WITH  DATA  THAT  RANGES  OVER  SEVERAL 

C  ORDERS  OF  MAGNITUDE. 

C 

C 

C  THE  FOttOMING  'DO  LOOP'  CALCUIATES  AN  BOCIVALEKT  VALUE 

C  FOR  THE  (X)  VARIABLES  BASED  ON  THE  .saSCTO)  FUNCTIONAL 
C  DEFINITION. 

C 

2014  DO  300  l-i,NP 
K-l 

R£AD(2,9002)ti,X(K)  ,y(|d 
iFdAx.Eo.i)  x(n«x(n 
IF(IAX.E0.3)  Xtl)>1.0/X(l> 
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X(I)*>DEXP(1.0/X(I)  ) 
X(I)»DEXP(-1.0/X(I) ) 


NAME 

MASK 

DATE 

TIME 


IF(IAX.EQ.3)  X(I)-DEXP(X(I) ) 
IFdAX.EQ.S)  X<I)»DIX)G10(X(I)> 
IF(IAX.EQ.7)  X(I)=DLOG(XlI) ) 
IF(IAX.EQ.9)  X(I)-DSIN(X(I) ) 
rF(IAX.EQ.4)  X(I)«»DEXP(-X(I) ) 
IFdAX.EQ.S) 

IFdAX.EQ.e) 

300  CONTINUE 

IF(REP.EQ.yES)  GO  TO  9777 
GO  TO  9888 
9777  READ(2,2006) 

READ(2,3046) 

READ(2,3047) 

READ(2,3048) 

READ(2,3049) 

READ(2,2006) 

READ<2,9006) 

READ(2,3094) 

READ(2,3120) 

IP(SELeCT(7)  .EQ.  '2')  GO  TO  3170 

READ (2, 3 051)  ICl.ICTPl 

READ(2,30S2)  IC2,ICTP2 

READ(2,3053)  IC3,ICTP3 

READ(2,30S4)  IC4,ICTP4 

READ(2,305S)  ICS.ICTPS 

READ (2, 3056)  1C6,ICTP6 

XCl-ICl 

XC2-IC2 

XC3-IC3 

XC4«rC4 

XC5-IC5 

XC6-IC6 

ICd)«XCl 

IC{2)«XC2 

lC(3)-iXC3 

IC(4)«XC4 

IC(S)-XC5 

:C(6)*XC6 

ICTPtD-ICTPl 

ICTP(2)-ICTP2 

ICTP(3>»ICTP3 

ICTp{4)»reTP4 

ICTP(5) -ICTPS 

lCTPt6)-ICTP6 

DO  9417  I-l,6 

ied)«icd)/(iooo.o«icTP(x) ) 

9417  CONTINUE 
GO  TO  3171 

3170  READt2,305l)ICl,lCTPl 
READ{2,31S6)IC2,1CTP2 
RBAD(2,3157) IC3.ICTP1 
READ ( 2 » 3 1 S a ) IC4 . ICTP4 
READ(2,31S9)IC5»lCrpS 
READ(2,3052) IC6,If^P6 
READ(2,3160)IC7.ICTP7 
RfiAD(2,5161)lca»ICTP<l 
READ(2,3162>IC9,ICTP9 


-APPENDIX  A— 


READ(2,3163) ICIO.ICTPIO 

READ(2,3053) ICll.ICTPll 

READ(2,3164) IC12,ICTP12 

READ(2,3165! IC13,ICTP13 

READ (2, 3166) IC14,ICTP14 

READ (2, 3 054) IC15,ICTP15 

READ<2,3055)IC16,ICTP16 

XC1«IC1 

XC2-IC2 

XC3-IC3 

XC4-IC4 

XC5-IC5 

XC6-IC6 

XC7-IC7 

XC8»IC8 

XC9-IC9 

XCIO-ICIO 

xcii-icn 

XC12-IC1-2 

xcn-ici3 

XC14..IC14 

XCIS-ICIS 

XC16-IC16 

ICTP(1)«ICTP1 

ICTP(2)"ICTP2 

ICTP(3)»ICTP3 

ICTP(4)«ICTP4 

ICTP>‘5)»ICTP5 

lCTPt6)-ICTP6 

ICTP(7) -ICTP? 

ICTP{8)-ICTP8 

ICTP(9)-ICTP9 

ICTP<10)»tCTP10 

ICTPdD-ICTPll 

ICTP(12)-ICTP12 

ieTPtl3>«ICrP13 

ICTP(14>«ICTP14 

ICTP(1S)«ICTP1S 

ICTP(16) •ICTP16 

IC<l)-SCl 

ic<2)«xe2 

ie<3)«XC3 

tCt4)«XC4 

lCtSl«XC5 

1056)  •xe6 

ICt7)-XC7 

tC(8)«XC8 

lCt9)«Xe9 

tC<10)«XC10 

ic(in-xcu 

IC(12)-XC12 

IC(131-XC11 

tC(U)«XC14 

IC51S)»XC1S 

tCd6)-XC16 

DO  3172  1*1. 1< 

ic({)«ic(i)/(iooo.o«tcTP<n) 
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3172  CONTINUE 
3171  CONTINUE 
REWIND  2 
9888  TYPE  2006 
306  FOBMATdX,'  '/) 

THE  USER  NOW  SELECTS  THE  ORDER  OF  THE  POLYNOMIAL  CURVE 
FIT  EQUATION.  THE  MAXIMUM  ORDER  IS  ONE  LESS  THAN  THE 
NUMBER  OF  (X)  DATA  VALUES. 

TYPE  2006 

TYPE  4137,ICCBS(1) .ICCBS(2) 

TYPE  2006 
TYPE  301 

301  FOBMATdX, ‘ENTER  THE  ORDER  OF  THE  DESIRED  POLYNOMIAL'/) 

TYPE  306 
TYPE  302 

302  FOBMATdX, '  (MAXIMUM  -  THE  NUMBER  OF  SCOIUM  CHLORIDE'/ 

C  CALIBRATION  CONCENTRATION  STANDARDS  -  1) '/) 

TYPE  306 
TYPE  2004 
ACCEPT  303,M 

303  FORMAT (12) 

TYPE  2006 

TYPE  4137,lCCBSd)  ,ICeBS(2) 

TYPE  2006 
DO  403  JK«1,14 
W(JK-l)«0,0 
S(JK-l)»0.0 
403  eOWTINUB 
W(0)aNP 
DO  401  I«l,NP 
DO  400  .I»1,2*M 
M(J)«W{J)+X(I)*»J 
5(J)«2(J>*V(I)''X(I)»»J 

400  COtWINUB 
J(0)-S(0)*Y{I) 

401  CONTINUE 
DO  !>00  1*1, M^l 
DO  SOO  3>1,M«1 
8(l,J)»N(J-2*l) 

D(I,M*2)»S(l-l) 

SOO  CONTINUE 

C  T!IB  FOLDOWINQ  SERIES  OF  NESTH)  DO  LOOPS  ACCOMPLISH  THE  FOLLONINCj 

C  THE  FIRST  IKCREHeNTS  THE  NESTED  DO  LOOP  SYSTEM  AND  THUS  SOLVES 

C  FOR  EACH  POLYNOMIAL  COEfFICIENT  ON  A  SUCCESSIVE  OASIS. 

C 

C  THE  SeCONL  INCREMENTS  THE  NESTQ)  DO  LOOP  SYSTEM  AND  THUS  DECREASES 

C  THE  MATRIX  SISE  OV  ONE  RCM  AND  COLUMN  EAdl  ITERAYON. 

C 

C  THE  THIRD  AND  FOURTH  TAKE  THE  ITH  RON  AND  DIVIDE  IT  BY  ITS  NTH 

C  COEFFICIENT  AND  THEN  MULTIPLIES  THE  ITH  MOM  OY  THE  NTH  COEFFICIENT 

C  OF  THE  (I«1)TH  MON,  FINALLY  IT  SUBTRACTS  THE  (I»I)TH  ROM  FROM  THE 
C  ITH  ROM. 

C 

C  THE  UST  TVO  DO  LOOPS  RESTRUCTURE  MATRIX  0  SO  THAT  IT  CONTAINS  THE 
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E  MATRIX  MINUS  THE  NTH  COLIWN  AND  THE  LAST  ROW. 

DO  50  N»1,M+1 
DO  30  11=1, M+I 
DO  30  Jl=l,H+2 
D(I1,J1)-B(I1,J1) 

30  CONTINUE 

DO  40  K=1,H 
IF(N-K)35,35,33 
33  L=1 

GO  TO  37 
35  L»2 

37  DO  31  I=1,H+1-K 
DO  31  J«1,M*3-K 

E(I,J)  =  (D(I,a)/’D(I.L)  )»D(l-t-l,L)-D<I  +  l,J) 

31  CONTINUE 

THE  COEFFICIENTS  FOR  THE  POLYNOMIAL  LEAST-SQUARES  CURVE  FIT  EQUATION 
ARE  CONTAINQ)  IN  THE  *AIN) *  ARRAY. 

DO  40  I»l,H+l-K 
DO  4Q.  J“i,M+2-K 
IF(J  .GE.  L)  GO  TO  45 
D<I,J)«E(I,J) 

G©  TO  40 

45  D{I,J)«E(I,J+1> 

110  CONTINUE 

AtN)-0(l,2)/D(l.U 

50  CQNflNUB 
WRlTBtJ, 20071 
WRITBd.SDH 

OUTPUT  THE  CALCULATB)  UNFORHATION  INTO  THE  CALCE.XX.X  FILE. 

51  F0WAT(1X,24X, 'THE  ORDER  OF  THE  OESIREP  POLYNONIAL*  *,12) 
WBJTEtJ.aOOS) 

WBITEtS,2006) 

IFdAX.BQ.l)  GO  TO  SJ 
irCIAX.^.S)  CO  TO  5> 

IFtlAX.KQ.J)  CO  TO  54 
irtIAX.eQ.4)  CO  TO  55 
IFUAX.BO.S)  GO  TO  S« 

SFdAS.EQ.S)  GO  TO  57 
Iff lAX.SOxT)  GO  TO  so 
;f flAX.EQ.e)  GO  TO  59 
If tIAX.BQ.9)  GO  TO  iO 

52  WRITBfJ.BU 

61  FO»HATtiX,J4X,*T«E  fOLYNORlAL  FtJ)*CTIOWL  OEnNITtOM  OF  THE 
C  VARIAOt.6  fX)  IS  IN  TESNS  Of  ♦  fX)  *//■<') 

SO  TO  89 
51  MRITEtl,62) 

62  FOH«AT(lX,24X, 'THE  POtVNOHIAL  FUNCTIONAL  DEFINITION  OF  THE 
C  VARIAfil.B  fX)  IS  IN  TERMS  OF:  (1/X) ’///) 

CO  to  89 
54  NRlTE(3,63) 

63  FOII«AT{lX,24X,’THE  I'OLVNCMIAL  F'tiMCTSONAL  DEFINITION  OF  THE 
C  VARIABLE  fXI  IS  IN  TERMS  OF:  ^EXPONENTIAL  tX))*///) 
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GO  TO  89 

55  WRITE(3,64) 

64  F0RMAT(1X,24X, 'THE  POLYNOMIAL 
C  VARIABLE  (X)  IS  IN  TEW4S  OF; 

GO  TO  89 

56  WRITE(3,65) 

65  F0RMAT(1X,24X, ‘THE  POLYNOMIAL 
C  VARIABLE  (X)  IS  IN  TERMS  OF; 

GO  TO  89 

57  WRITE(3,66) 

66  F0RHAT(1X,24X,'THE  POLYNOMIAL 
C  VARIABLE  (X)  IS  IN  TERMS  OF: 

GO  TO  89 

58  WRITE(3,67) 

67  FORMAT(iX,24X, 'THE  POLYNOMIAL 
C  VARIABLE  (X)  IS  IN  TERMS  OF: 

GO  TO  89 

59  WRITER, 68) 

68  F0RHAT<1X,24X, 'THE  POLYNOMIAL 


FUNCTIONAL  DEFINITION  O'"  THE 
(EXPONENTIAL  (-X) )'///) 


FUNCTIONAL  DEFINITION  OF  THE 
(EXPONENTIAL  (I/X))'///) 


FUNCTIONAL  DEFINITION  OF  THE 
(EXPONENTIAL  (-1/X> )'///) 


FUNCTIONAL  DEFINITION  OF  THE 
(NATURAL  LOG  (X)  '///) 


FUNCTIONAL  DEFINITION  OF  THE 


C  VARIABLE  (X)  IS  IN  TERMS  OF:  (LOG  BASE  TEN  (X)'///) 

GO  TO  89 
60  WRITB(3,69) 

69  F0RMAT(1X,24X,'THE  POLYNOMIAL  FUNCTIONAL  DEFINITION  OF  THE 
C  VARIABLE  (X)  IS  IN  TERMS  OF:  (SIN  (X))'///) 

GO  TO  89 

89  CONTINUE 

DO  90  I»l,M+l 
TYPE  70,I,A(I) 

71  FORMAT( IX. 24X, 'COEFFICIENT  NUMBER', 13, 'IPBIS.S) 

70  FORHATd.X, 'COEFFICIENT  NUMBER’ ,  J  3, '•  '.IPEIS.S) 
WR7TE(3,71)I,A<I) 

90  eOMTlNOE 
TYPi  2006 

TYPE  4137,lCCBS<l),ieC8S(2» 

TYPE  2006 
TYPE  2006 
16  TYPE  19 

19  FORMAttJS^'DO  YOU  WANT  ,A  LIST  OP  RESIDUALS  ?'/ 

C  (ANSfc'SB  YES  OR  NO>  SIX,*) 

ACCEPT  10S5,COMnEN 
TYPE  2036 

IPiCOMSEN  YES)  GO  fO  10 
IFtCOMMEN  .SO.  NO)  GO  TO  17 

.ME.  YSS  .AND.  COMMEN  .NE.  NO)  (»  TO  16 

TYPE  2006 

TYPE  4i37,lCCDS(i) ,ieeBS(2) 

TYS^'  2006 
10  «Pm<l,3007) 

KRm(J,23) 

23  fORMAtdX.24X. ‘THE  RESIDUALS  ARE  CALCULATED  BY  THE  FOLLOWING*) 
KRtTet3.)3) 

13  row  AT  (U,  2  *X, 'EQUATION:  tY(DATA)«y(EStINATED)  1 '(’/) 

4«)tTE(l,3  006) 

W81TE(3.24) 

24  fO»IAT<lS,iex,'SA«PU  NUMBER*  ,17X, ‘YOATA’ ,23S, 'YCESTIMATliD) ' , 
Ci9X,*R£SiOUAL'.') 
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THE  RESIDUALS  ARE  CALCULATED  USING  THE  POLYNOMIAL  CURVE 
FIT  EQUATION  AND  THE  ORIGINAL  DATA. 

SUMRES=0.0 
DO  15  K=1,NP 
THEORy=0.0 
DO  21  1=1, M+I 

IF(I  .EQ.  1  .AND.  X{K)  .EQ.  0.0)  GO  TO  20 
THEORY=A(I) *X(K) **(I-1) +THEORY 
GO  TO  21 

20  THEORY=A(I)*1.0**(I-1)+THEORY 

21  CONTINUE 

CALCULATE  THE  SOM  OF  THE  SQUARES  RESIDUAL. 

RESID=Y<I<)  -THEORY 
SRESD=n!.oir  *RESID 
SUHRES=.l  E'.D  f  SOMRES 
TYPE  10,K,  -.  .SID 

10  FORMAT (IX, 'SAMPLE  NUMHER  ' , 13 , 12X , ' RESIDUAL^  ',1PE15.5) 
WRITE(3,22)K,Y(K) , THEORY, RESID 

22  FORMAT(1X,23X,I2,20X,1G13.6,18X,1G13.6,16X,1G13.6) 

15  CONTINUE 

TYPE  2006 
TYPE  72,SUMRES 

72  FORMAT (IX, 'SUM  OF  THE  SQUARE  RESIDUALS=  ',1PE15.5) 

TYPE  2007 

)'miTE(3,2007) 

WRITE(3,73)SUHRSS 

73  F0RMAT(1X,24X, 'SUM  OF  THE  SQUARE  RESIDUALS-  ’,1PE15.5) 

17  CONTINUE 

CALCULATE  INDIVIDUAL  EXERCISE  PROTECTION  FACTORS 


NOTE:  SINCE  THE  LEAK  MEASURING  SENSITIVITY  OF  THE  SODIUM 
CHLORIDE  RQFT  INSTRUMENT  IS  ONE  PART  IN  TEN  TO  THE 
SIXTH,  ANY  EXERCISE  SCALED  INTEGRATOR  COUNT  VALUE 
YIELDING  A  PF  GREATER  THAN  1.0E-H06,  WILL  BE 
REPORTED  AS  l.OE+06.  REPORTING  A  MASK  LEAKAGE  GREATER 
THAN  1.0E-F06  WOULD  BE  ERRONEOUS.  ANY  EXERCISE  SCALED 
INTEGRATOR  COUNT  VALUE  YIELDING  A  PROTECTION  FACTOR 
GREATER  THAN  1.0  WILL  BE  REPORTED  AS  l.OE-EOO. 

TYPE  2006 

TYPE  4137,ICCBS(1) ,ICCBS{2) 

TYPE  2006 

IF(SELECT(7)  .EQ.  '2')  GO  TO  31/3 
IDL-6 
IDLP-7 
GO  TO  3174 

3173  CONTINUE 
IDL-16 
IDLP-17 

3174  CONTIN'  E 
PFEST>  J.O 

DO  9173  1=1, IDL 
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PF{I)=0.0 
9173  CONTINUE 
CAI,CPF  =  0.0 
DO  3057  1  =  1, IDL 

IFdCd)  .LE,XXAS(7)  )  GO  TO  9227 
IFdCd)  .GE.XXASd)  )  GO  TO  4136 
IFdAX.EQ.l)  PFEST  =  ICd) 

IFdAX.EQ.2)  PFEST=1.0/IC(I) 
rFdAX.EQ.3)  PFEST=DEXPdC(I)  ) 

IFdAX.EQ.4)  PFEST=DEXP(-ICd)  ) 

IFCIAX.EQ.S)  PFEST=DEXPd.O/IC(I)  '< 

IFdAX.EQ.6)  PFEST=DEXP(-1.0/ICd)  ) 

IFdAX.EQ.7)  PFEST=DIiOGdC(I)  ) 

IFdAX.EQ.8)  PFEST=DLOG10(ICd)  ) 

IFdAX.EQ.9)  PFEST=DSINdC<I)  ) 

CALCPF=0.0 
GO  TO  9228 

9227  CALCPF=0. 000001 
GO  TO  3058 

4136  CALCPF=0.0 
GO  TO  3058 

9228  DO  3058  J=1,M+1 
IF(J.EQ.1.AND.PFE5:T.EQ.0.0)  go  to  3059 
CALCPF=A(J) *PFEST**{J-1)+CALCPF 

GO  TO  3058 

3059  CALCPF=A(J)*1.0**(J-1)+CALCPF 
3058  CONTINUE 

C 

C  BY  DEFINITION,  THE  PROTECTION  FACTOR  (PF)  IS: 

C 

C  PF  =  (CHALLENGE  ATMOSPHERE  CONCENTRATION) / (MASK  LEAK  CONCENTRATION) 

C 

C 

IF(CALCPF.NE,0,0)  GO  TO  4135 

PFd)=0.0 

GO  TO  3057 

4135  PFd)=1.0/CALCPF 
PFd)»1.0/CALCPF 
3057  CONTINUF 
C 

C  CALCULATE  AN  OVERALL  ARITHMETIC  AVERAGE  PROTECTION  FACTOR  FOR  ALL 
C  EXERCISES. 

C 

KOUNT=0 
PFSUM=0,0 
no  3060  MT=1,IDL 
KOUNT=KOUNr  + 1 
PFSUM«PFS'''-.  i  PF(HT) 

3060  CONTINUE 

PF ( IDLP) =PFSUM/KOUWr 
C 
C 

C  CALCULATE  AN  OVERALL  TIME  WEIGHTED  AVERAGE  PROTECTION  FACTOR 
C  FOR  ALL  L'XERCI.SES. 

C 

C 

WPF=0.0 


84 


-APPENDIX  A- 


KKOUm'=0 
PPSUM=0.0 
DO  3194  IHT=1,IDL 
KKOUMr=KKOUNr  +  ICTP<IMT) 

PPSUM=PPSUH  +  (PF<IHT) *ICTP(IMT> ) 

3194  CONTINUE 

WPF  =  PPSUH/KKOOm' 

WRITE{3 ,3061) 

3061  FORMAT (IHl) 

WRITEO  ,3062) 

3062  FORMAT(6X, 'THE  DESCRIPTIVE  AND  PROTECTION  FACTOR  CALCULATIONS:'/) 
WRITE(3 ,9542) 

9542  FORMAT(6X, 'NOTE:  ANY  PROTECTION  FACTOR  THAT  IS  LISTED  AS'/ 

C'  l.OE+06  HAS  BEEN  ASSIGNEE  THIS  VALUE  BY  DEFAULT'/ 

C  BECAUSE  THE  SENSITIVITY  OF  THIS  RQFT  INSTRUMENT  IS'/ 

C  AT  MOST  ONE  PART  IN  TEN  TO  THE  SIXTH.  THE  INTEGRATOR'/ 

C  COUNT  VALUE  FOR  A  PARTICULAR  EXERCISE  IN  QUESTION’/ 

C  IS  MERELY  REPRESENTATIVE  OF  INTEGRATING  THE  ELECTRICAL'/ 

C  NOISE  AND  THE  TRUE  PROTECTION  FACTOR  IS  INDEED  LESS') 

WRITE (3, 9 543) 

9543  FORMAT {6X, 'THAN  l.OE+06.  ANY  EXERCISE  SCALED  INTEGRATOR'/ 

C  COUNT  VALUE  YIELDING  A  PROTECTION  FACTOR  GREATER'/ 

C  'THAN  1.0  WILL  BE  REPORTED  AS  O.OE-01.'//) 

VmiTE(3,3176) 

3176  FORMATdX,  '  '//) 

TYPE  3062 
TYPE  2006 
TYPE  9542 
TYPE  9  54’ 

TYPE  2000 

TYPE  413V, :CCBS(1),ICCBS(2) 

TYPE  2005 
WRITE (3, 3 046) NAME 
WRITE (3, 3 047) MASK 
WRITE(3 ,3048)DATE 
WRITE(3,3049)TIME 
WRITE(3,3176) 

TYPE  3 046, NAME 
TYPE  3 047, MASK 
TYPE  3 04 8, DATE 
TYPE  3049, TIME 
TYPE  2006 

TYPE  4137,ICCBS{1) ,ICCBS(2) 

TYPE  2006 
WRITE(3 ,3063) 

TYPE  3063 

3063  FORMAT(6X, 'EXERCISE' ,29X, 'PROTECTION  FACTOR'/) 

IF(SELECT{7)  .EQ.  '2')  GO  TO  3175 

WRITE(3 ,3064)  PF(1) 

TYPE  3064, PF(’) 

3064  FORMAT (6X, 'NORMAL  BREATHING  STRAIGHT  AHEAD ' ,6X , IPEl 2 . 1) 

WRITER  ,3065)  PF(2) 

TYPE  3065,PF(2) 

3065  FORMAT(6X, 'DEEP  BREATHING  STRAIGHT  AHEAD ', 8X , 1PE12 . 1) 

'.<'RITE(3 ,3066)  PF(3) 

TYPE  3066  , PF (3) 

3066  FORMAT(6X, 'TALKING ' ,30X, IPEl 2 . 1) 
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WRITE(3  ,3067)  PF(4) 

TYPE  3067, PF(4) 

3067  FORMAT(6X, 'SIDE-TO-SIDE  HEAD  MOVEMENTS'/ 

C  (DEEP  BREATHING) ' ,21X,1PE12.1) 

WRITE(3,3068)  PF(5) 

TYPE  3068, PF(5) 

3068  FORMAT (6X, 'UP-AND-DOWN  HEAD  MOVEMENTS'/ 

C  (DEEP  BREATHING) ' ,21X,1PE12.1) 

WRITE(3,3069)  PF(6) 

TYPE  3069, PF(6) 

3069  FORHAT(6X, 'FACIAL  GRIMACING ', 21X , 1PE12 . 1////) 

WRITE(3,3070)  PF(7) 

TYPE  3C70,PF(7) 

3070  FORMAT(lHO,5X, 'OVERALL  ARITHMETIC  AVERAGE  PROTECTION  FACTOR  '/ 

C  FOR  ALL  CATEGORIES  OF  EXERCISES  ACTUALLY  PERFORMED  =  ' 

C,1PE8.1) 

TYPE  2006 
WRITE(3,2006) 

TYPE  3195,WPF 
WRITE(3,3195)  WPF 

3195  FORMAT (1H0,5X, 'OVERALL  TIME  WEIGHTED  AVERAGE  PROTECTION  FACTOR  '/ 
C  FOR  ALL  CATEGORIES  OF  EXERCISES  ACTUALLY  PERFORMED  =  ' 

C,1PE8.1) 

GO  TO  3177 
3175  CONTINUE 

WRITE(3,3178)PF(1) 

TYPE  3178, PFd) 

3178  PORMAT(6X, 'NORMAL  BREATHING  STRAIGHT  AHEAD ', 8X, 1PE12 . 1) 
WRITE(3,3179)PF(2) 

TYPE  3179, PF(2) 

3179  FORMAT (6X, 'NORMAL  BREATHING  LEFT ' , 18X , 1PE12 .1) 

WRITE(3,3180)PP(3) 

TYPE  3180, PF  (3) 

3180  FORMAT (6X, 'NORMAL  BREATHING  RIGHT ', 17X, 1PE12 . 1) 

WRITE(3,3181)PF(4) 

TYPE  3181, PF(4) 

3181  FORMAT (6X, 'NORMAL  BREATHING  DOWN' ,18X , IPEI 2 . 1) 

WRITE(3,3182)PF(5) 

TYPE  3182, PF(5) 

3182  FORMAT (6X, 'NORMAL  BREATHING  UP ' ,20X , 1PE12 . 1) 

WRITE(3,3183)PF(6) 

TYPE  3183, PF  (6) 

3183  FORM  AT  ( 6  X, 'DEEP  BREATHING  STRAIGHT  AHEAD  '  ,  1  OX ,  IPEl  2 . 1) 
WRITE(3,3184)PF(7) 

TYPE  3184, PF(7) 

3184  FORMAT(6X, 'DEEP  BREATHING  LEFT ' , 2 OX , IPEl 2 . 1) 

WRITE(3,3185)PF(«) 

TYPE  3185, PF(8) 

3185  FORMAT(6X, 'DEEP  BREATHING  RIGHT ', 19X , 1PE12 . 1) 

WRITE(3,3186)PF(9) 

TYPE  3186, PF(9) 

3186  FORMAT (6X, 'DEEP  BREATHING  DOWN' ,20X , 1PE12 . 1) 

WRITE(3,3187)PF(10) 

TYPE  3187, PFdO) 

3187  FORMAT(6X, 'DEEP  BREATHING  UP ' ,22X , IPEl 2 . 1) 

WRITE(3,3188)PF(11) 

TYPE  3188, PFdl) 
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3188  FORHATtex, 'TALKING ' ,32X. 1PE12. 1) 

WRITE(3,3189)PP(12) 

TYPE  3 189, PF  (12) 

3189  FORMAT(6X, 'FACIAL  GRIMACING ' ,23X, 1PE12 . 1) 

WRITE(3,3190)PF(13) 

TYPE  3190, PF  (13) 

3190  FORHAT(6X, 'SIDE-TO-SIDE  HEAD  MOVEMENTS'/ 

C  (NORMAL  BREATHING) ',21X,1PE12.1) 

WRITE(3,3191)PF(14) 

TYPE  3191, PF  (14) 

3191  FORMAT(6X, 'UP-AND-DOWN  HEAD  MOVEMENTS'/ 

C  (NORMAL  BREATHING) ' ,21X,1PE12.1) 

WRITE(3,3192)PF(15) 

TYPE  3192, PF (15) 

3192  FORMAT(6X, 'SIDE-TO-SIDE  HEAD  MOVEMENTS'/ 

C  (DEEP  BREATHING) ',23X,1PE12.1) 

WRITE(3,3193)PF(16) 

TYPE  3193,PP(16) 

3193  FORMAT(6X, 'UP-AND-DOWN  HEAD  MOVEMENTS'/ 

C  (DEEP  BREATHING) ' ,23X,1PE12.1////) 

WRITE(3,3070)PF(IDLP) 

TYPE  3070,PF(IDLP) 

TYPE  2006 
WRITE(3,3176) 

TYPE  3195, WPF 
WRITE(3,3195)WPF 
3177  CONTINUE 
TYPE  2006 

TYPE  4137,ICCBS(1) ,ICCBS(2) 

TYPE  2006 

FORM  THE  ARRAYS  FOR  PLOTTING  PURPOSES  AND  LOAD  THEM  INTO  THE 
FILE  CALLED  GRPHX.XXX. 


FIND  THE  MINIMUM  AND  MAXIMUM  OF  THE  (X)  DATA  VALUES 

AND  FORM  A  DELTA  ELEMENT  SO  THAT  INTERMEDIATE  VALUES 

CAN  BE  CALCULATED  USING  THE  POLYNOMIAL  CURVE  FIT  EQUATION. 

TYPE  3061 
DO  2  I»1,NP 
GRPHKD-XXASd) 

2  CONTINUE 
Nl-NP-1 
DO  8  I»1,N1 
Il-I  +  1 
DO  8  J-I1,NP 

IF(XAS{J)  .GE.XASd)  )  GO  TO  8 
XSAVE-XAS ( I ) 

XAS(I)«XAS(J) 

XAS(J)-XSAVE 
8  CONTINUE 
XMIN-XASd) 

XNAX-XAS(NP) 

DEL-(XMAX-XMIN)/400.0 

XXX-XMIN-DEL 

DO  7  KA»1,401 
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XXX=XXX  +  DEL 
XX(KA)=XXX 
XXXX=XXX 

IF(IAX.EQ.l)  XEST*XXXX 
IFtIAX.EQ.2)  XEST-l.O/XXXX 
IF(IAX.EQ.3'  XEST»DEXP(XXXX) 

IF(IAX.EQ.4)  XEST-DEXP(-XXXX) 

IFdAX.EQ.S)  XEST-DEXPd.O/XXXX) 

I F (I AX . EQ . 6)  XEST-DEXP ( -1 . 0/XXXX) 

IFdAX.EQ.7)  XEST*DLOG(XXXX) 

IFdAX.EQ.8)  XEST-DLOGIO(XXXX) 

IFdAX.EQ.9)  XEST»DSIN(XXXX) 

BEGIN=0.0 
DO  6  KB=1,H+1 

IF(KB.EQ.l.AND.XEST.EQ.O.O)  GO  TO  5 
BEG1N=A(KB) *XEST** (KB-1) +BEGIN 
GO  TO  6 

5  BEGIN=A(KB)*1.0*’*(KB-1)+BEGIN 

6  CONTINUE 
Yy(KA)=BEGIN 

7  CONTINUE 
DO  4  1=1,401 
GRPH2d)“XXd) 

4  CONTINUE 
DO  3  1=1/401 
GRPH2d  +  401)=YYd) 

3  CONTINUE 
DO  1  isifNP 
GRPHld+NP)=Yd) 

1  CONTINUE 

OUTPUT  THE  INFORMATION  THAT  CAN  BE  USED  FOR  GENERATING  A 
CALIBRATION  PLOT  INTO  THE  PILE  CALLED  GRPHX.XXX. 

WRITEd  ,2007) 

WRITEd  ,9) 

9  FORMAT dX, 2 4X, 'THE  ARRAYS  FOR  PLOTTING'///) 

WRITEd, 2007) 

WRITEd, 3006) 

3006  F0RMAT(1X,24X, 'GRPHl  ARRAY  CONTAINS  THE  SODIUM  CHLORIDE 
C  CALIBRATION  STANDARD  DATA'///) 

WRITEd, 3004) 

3004  FORMATdX,24X, 'SAMPLE  NUMBER' ,  19X, ' VALUE ' /) 

DO  34  I-1,NP*2 

WRITE (1, 3005) I, GRPHl (I ) 

3005  FOBMATdX,28X,I3,22X,lG13.6) 

34  CONTINUE 

WRITEd, 2007) 

WRITEd, 3007) 

3007  F0RMAT(1X,24X, 'GRPH2  ARRAY  CONTAINS  THE  ESTIMATED  VALUES'///) 
WRITEd, 3004) 

DO  32  1=1,802 
WRITEd, 3005)1, GRPH2d) 

32  CONTINUE 
29  TYPE  2006 

TYPE  4137,ICCBSd)  ,ICCBS(2) 

TYPE  2006 
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TYPE  27 

27  FORHATdX, 'DO  YOU  WISH  TO  USE  THE  SAME  SODIUM  CHLORIDE  STANDARD'/ 
C  CALIBRATION  VOLTAGE  MEASUREMENTS  AND  EXERCISE  INTEGRATOR  COUNT'/ 
C  DATA  BUT  CALCULATE  A  DIFFERENT  DEGREE  OR  FORM  OF  THE'/ 

C  POLYNOMIAL  CURVE  FITTING  FUNCTION?  (ANSWER  YES  OR  NO)  ',2X,S) 
ACCEPT  1002, REP 
IF(REP.EQ.YES)  GO  TO  6000 
IF(REP.EQ.NO)  GO  TO  28 
IF(REP. NE.YES.AND.REP.NE.no)  GO  TO  29 

28  CONTINUE 
CLOSE <UNIT=2) 

CLOSE (UNIT=3) 

CLOSE <UNIT=1) 

9674  TYPE  2  007 

TYPE  4137,ICCBS<1) ,ICCBS(2) 

TYPE  2006 
TYPE  14 

14  FORMATdX, 'DO  YOU  WISH  TO  CALCUIATE  PROTECTION  FACTORS  FOR'/ 

C  ANY  OR  ALL  OF  THE  FOLLOWING  CONDITIONSs  1)  A  DIFFERENT'/ 

C  SUBJECT  2)  A  DIFFERENT  SET  OF  SODIUM  CHLORIDE  CALIBRATION'/ 

C  STANDARD  VOLTAGE  MEASUREMENTS  3)  A  DIFFERENT  SET  OF'/ 

C  EXERCISE  INTEGRATOR  COUNT  DATA.'/) 

TYPE  9999 

9999  FORMATdX,  '  (ANSWER  YES  OR  NO)',2X,S) 

ACCEPT  1002, COM 
IF(COM.EQ.YES)  GO  TO  5101 
IF(COH.EQ.NO)  GO  TO  9672 
IF(COM. NE.YES.AND.COH.NE.no)  GO  TO  9674 

9672  CONTINUE 
TYPE  2  007 

TYPE  4137,ICCBS(1) ,ICCBS(2) 

TYPE  2007 
TYPE  9599 

9599  FORHATdX, 12X, 'JOB  SUCCESSFULLY  COMPLETED',//////) 

STOP 

END 


SUBROUTINE  CLEAR  IS  CALLED  WHEN  ENTERRING  DATA  IN  THE 
SODIUM  CHLORIDE  CALIBRATION  CONCENTRATION  STANDARDS 
SECTION  AND  EXERCISE  INTEGRATOR  COUNT  SECTION.  BY 
CALLING  THIS  SUBROUTINE  THE  OPERATOR  CAN  MAKE  TtYO 
DATA  ENTRIES  ON  THE  SAME  LINE.  SUBROUTINE  CLEAR  ERASES 
THE  LINE  INWHICH  THE  FIRST  DATA  ENTRY  WAS  MADE  AND 
RETYPES  THAT  LINE,  INCLUDING  THE  FIRST  DATA  ENTRY; 

THIS  ALLOWS  A  SECOND  ENTRY  TO  BE  MADE  ON  THAT  LINE. 


SUBROUTINE  CLEAR(LINES) 
BYTE  A(3) 

Ad) -27 
A(2)-65 
A(3)-75 

IF  (LINES. EQ.O)  LINBS-1 
DO  1  I-l, LINES 
1  TYPE  4,A(1)  ,A(2)  ,A(1)  ,A(3) 
LINES-0 
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RETURN 

4  FORMAT  (1H+,4A1,$) 
END 


APPENDIX  B: 


DATA. XXX  File  Contents  for  Data  in  Table  7 
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TIIE  MUtlUKR  nf  SdOItItl  CIlU3RinE  CAMllRATIOH  COtlCEHTRATIOH  STAHDA1!DR=  7 


RAMPt.E  NtIMUri! 


vol.TAt’.r; 

ni;Ar.!iR[;ni’.!iT 

(Ill  voi.t:;) 

(  X  DATA  ) 


r.ODluti  ciii-nRit)E 

CAE  I  DRAT  10 11 
COllCEMTRATIOn 
(  Y  DATA  ) 


1 

2 

3 

4 

5 

6 
7 


3 .3ftOOO 
2.91^00 
2.31000 
1.50000 
0.545000 
0.165000 
0.105000 


1.00000 

0.100000 

O.lOOOOOE-01 

O.lOOOOOE-02 

O.lOOOOOE-03 

O.lOOOOOE-04 

O.1OQ00OE-O5 


MIIEtlXT  llAllEs  CAPTAIN  EDWARD  D.  KOI.EflAR,  JR. 
TYl'O  Of  HA.SI<:  USA  j  tU7  -  MEDIUM  -  NO  OIASSE.S 
DATE  TESTED  t  0  APR! I,  1980 
Tim:  TIXRTF.D  1  1330  HOURS 


EXfECtliE  lUTECRATOR  COUtfT  DATAj 


■■  EEIX'ISE 

lUTRORATOU  0711 OT 

TIME  PERIOD 
(III  SECONDS) 

l!''RMAt.  DnEATIUllC.  STRAUHIT  AHEAD 

3904 

10 

DEEP  HEEATlIIUt:  STRAIGHT  AHEAD 

4751 

10 

TAt.K  t  IIG 

SlDE-TO-iilDE  HEAD  MOVRHENT.S 

4P20 

10 

(Dfn-  HREAlllltlG) 

Up- AtID-nOWtl  HEAD  MOVEMENTS 

397C 

10 

(DEEP  HREAl'HIlIf:) 

401G 

10 

ea>-i,m,  ghimacuk; 

4937 

10 

^  sum4mwujm 


\ 

( 

i 


i 


f  smmm  s»m  bummuc  nmo 
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If 


■L'lit;  oupin;  or  tiu;  Di;r:inf:D  i‘ni.YiioinAt.=  c 

rur.  pot.YiUiiiTAr.  rutKTioM,M,  nKrjtiiTion  or  the  VAniM-i.i;  (x)  is  iii  terms  of:  (expohehtial  tx)) 


(•oKmciriiT  mninru 
rnEFi'H'irMT  tiuir.ci: 

roEFFiriisiT  '.luiinEU 
I'onmciEiiT  I'.unnER 
corFFiciritT  HUfiriEi: 
rnr.mcir.tiT  tuiiiiiEi: 
mr.FnciF.tr  tutiwrr. 


1=  --l.  l7‘J04F-05 

2=  -7.114r,7E-05 

3=  F.25r.C3E-0< 

4=  -2.f.nQ05E-05 

5=  lJI0ni‘.B-06 

6=  -K71S4U-.-07 

7=  7,.<)3'!C6c-09 


-'lit:  Ri:,;  lotiAi..';  .m;i;  c.M.cuiTiTrn  i'.y  the  fouu>uujg 

mUiV-'U'tt:  lY<tlATA)-Y(Er.Vtl!ATll»  1 


•i.r  tniiirrr 


VtUSTUIATEO) 


RESIDUAL 


I 

T 

i 


1 .A 00 AO 
tl.lQOOt'O 

o.iaonoftj:-oi 

o,»uf>ocor.-o2 

0.  lOlMlOOE-OJ 
O.ltU'OOOtSOil 


1.00000 
0. 10QOOO 
O.lQOOOOC-01 
0. JOftOOOE'OJ 
0.10Q00<5E-03 
0.10CI034E'-04 
0.l00413e-QS 


-0.403548^-07 
-o.iQomc-tt? 
-0.7n342B-0C 
-0.401‘)2!>E-Ot; 
-0.39a294E'0r 
-  0.34  060  2K-0r. 
-0.4325028'OC 


.'ui;:  tir  -."tin  stw.rt:  nisiuuA»..‘;“  2.oo«nfc-i5 


HttTHlIiM  DUS  “  mjSD 
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Ti'.i.  Di;rcni  protfctio!'  rACi’or.  cw,ci'U*.7io;'r: 


!Jr"r.:  A —  F-'C^OR  ’c:  uTI'TII 

}.o::+G'.  !■;.;•  eff"  affigfe:'  ?!'i?  VAi.nr  r.Y  Drr.-.:::;r 

■  ''r  *■'  qp"  r  j 'J*  jy  J QP  "•pjq  pnp*^  T ;"  P""  •  pip’  \ 


‘IC'F’’  O'.'F  ^APT  IP  TCP  TO  TPC  '"Cl  I'CCCI'/i'^CF. 

COPT  vM.ur  Fcn  .-■.  piw-Ticvu.?.  cxcRcicr.  i::  cjcctio:: 
pcncLv  p.cpp.EsnT.'.Tivc  or  iptlg-t.ticg  tmc  rccccnic/c, 

'.'.-'zcr.  AC.C-  THE  TRUF  PPOTECTIOI.  FACTOR  IS  niDEE.’'  CEfS 
tha;:  i.OE+of. .  m:y  exet-cisf  scaled  niTEGp.ATor. 

COOT  VALUE  YIELDING  A  PROTECTIOl!  FACTOR  GREATER 

tha:;  1.0  i;ill  ee  reported  as  o.oe-oi. 


SUPJECT  RAIIE:  CAPTAIN’  EDWARD  S.  KOLESAR,  JP . 
TYPE  or  MASK;  USA  ;  M7  -  MEDIUM  -  NO  GLf.SSE:’ 
DiATE  TESTED  ;  5  APRIL  19  80 

TIME  TESTED  ;  1330  HOURS 


EXERCISE  PROTECTlOtl  PACTOR 

NORMAL  r-REATIUHC  STRAIGHT  AHEAD  l.eE*0.I 

deer  EPEATHING  STRAIGHT  AHEAD  1.3E*0< 

TALKING  1.3E+0.’ 

.'•IDE-To-sioE  HEW)  MOVEttENTS 

(DEEP  SREI.TStlNCl  i.7E+!'- 

UP-AMD'tXiWN  HEAP  MOVEMENTS 

(DEEP  BREATSUNM)  J.7E-*S< 

FACIAL  GRIMACING 


overall  AJUTflMilTIC  AVERAGE  PROTECTIOt!  PACTOR 

FOR  ALL  CATEGORIES  OP  EXERCISES  ACTOALLV  PERPOWt©  •  l.SE+04 


OVERALL  TIMT  WEIGHTED  AVERAGE  PROTECTIO*.'  TACTOR 
I-OB  ALL  CATEGORIES  OF  EXERCISES  ACTGALLV  PSRFORRSS 


1.5E»04 


APPENDIX  0; 

GRPHX.XXX  File  Contents  for  Use  with  NACLGRAPH .FTN  Program 
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'i’f.’r  Arp.'.YP  FOR  ?LO'??r;!G 


GPPiil  ARa\Y  COHTAIRS  7HE  SODIUM  CHLORIDE  CALIBRATION  STAJIDARD  DATA 


S'..''.PLr  cuiiocr. 


VALUE 


1 


6 

7 

Q 

9 

10 

U 

n 

1 


j.iP.COC 

2.91500 

2.31000 

1.5COOO 

0.545000 

0.165000 

0.105000 

1.00000 

0.100000 

0.1000002-01 

O.lOOCOOB-02 

O.lOOOOOE-03 

0.100000E-0< 

O.lOOOOOE-05 


G.’'vPH2  array  cohtaihs  thr  estimated  values 


SAMPU:  N UNDER 


VALUE 


1 

2 

J 

4 

5 

6 
7 
(I 
9 

10 


0.105000 

0.113187 

0.12137S 

0.139562 

0.137750 

0.145936 

0.1S4U5 

0.162313 

0.170500 

0.17B607 


S' 
fe ' 
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11 

0 .1S6C75 

12 

0.195063 

13 

0.203250 

14 

0.211438 

15 

0.219625 

16 

0.227812 

17 

0.236000 

le 

0.244188 

19 

0.252375 

20 

0.260563 

21 

0.268750 

22 

0.276930 

23 

0.235125 

24 

0.293313 

25 

0.301500 

26 

0.309687 

27 

0.317875 

28 

0.326063 

29 

0.334250 

30 

0.342438 

31 

0.350625 

32 

0.358813 

33 

0.367000 

34 

0.375188 

35 

0.383375 

36 

0.391563 

37 

0.399750 

38 

0.407937 

39 

0.416125 

40 

0.424313 

41 

0.432500 

42 

0.440688 

43 

0.448875 

44 

0.457063 

45 

0.465250 

46 

0.473438 

47 

0.481625 

48 

0.489813 

49 

0.498000 

50 

0.506187 

51 

0.514375 

52 

0.522563 

53 

0.530750 

54 

0.S38938 

55 

0.547125 

56 

0.535313 

57 

0.563500 

58 

0.571688 

59 

0.579875 

60 

0.588063 

61 

0.596250 

62 

0.604438 

63 

0.612625 

64 

0.620813 

65 

0,629000 

66 

0.6371BB 

67 

0.645375 

66 

0.653563 

69 

0.661750 

70 

0.669937 

71 

0.6781  .IS 
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72 

73 

74 

75 

76 

77 

78 

79 

80 
31 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 


0.686312 
0 .694500 
0.702688 
0.710875 
0.719063 
0.727250 
0.735438 
0.743625 
0.753813 
0.760000 
0.768188 
0.776375 
0.784563 
0.792750 
0.800938 
0.809125 
0.817313 
0.825500 
0.833688 
0.841875 
0.850063 
0.858250 
0.866438 
0.874625 
0.882812 
0.891000 
0.899188 
0.907375 
0.915563 
0.923750 
0.931938 
0.940125 
0.948313 
0.956500 
0.964688 
0.972875 
0.981063 
0.989250 
0.997438 
1.00563 
1.01381 
1.02200 
1.03019 
1.03838 
1.04656 
1.05475 
1.06294 
1.07113 
1.07931 
1.08750 
1.09569 
1.10388 
1.11206 
1.12025 
1.12844 
1.13663 
1.14481 
1.15300 
1.16119 
1.16938 
1.17756 
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133 

1.18575 

134 

1.19394 

135 

1.20213 

136 

1.21031 

137 

1.21850 

138 

1.22669 

139 

1.23488 

140 

1.24306 

141 

1.25125 

142 

1.25944 

143 

1.26763 

144 

1.27581 

145 

1.28400 

146 

1.29219 

147 

1.30037 

148 

1.30856 

149 

1.31675 

150 

1.32494 

151 

1.33312 

152 

1.34131 

153 

1.34950 

154 

1.35769 

155 

1.36588 

156 

1.37406 

157 

1.38225 

158 

1.39044 

159 

1.39863 

160 

1.40681 

161 

1.41500 

162 

1.42319 

163 

1.43138 

164 

1.43956 

165 

1.44775 

166 

1.45594 

167 

1.46413 

168 

1.47231 

169 

1.48050 

170 

1.48869 

171 

1.49688 

172 

1.50506 

173 

1.51325 

174 

1.52144 

175 

1.52963 

176 

1.53781 

177 

1.54600 

178 

1.55419 

179 

1.56238 

180 

1.57056 

181 

1.57875 

182 

1.58694 

183 

1.59513 

184 

1.60331 

185 

1.61150 

186 

1.61969 

187 

1.62788 

188 

1.63606 

189 

1.64425 

190 

1.65244 

191 

1.66063 

192 

1.66881 

193 

1.67700 
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194 


1.68519 


195 

196 

197 
190 

199 

200 
201 
202 

203 

204 

205 

206 
207 
200 

209 

210 
211 
212 


1.69338 

1.70156 

1.70975 

1.71794 

1.72613 

1.73431 

1.74250 

1.75069 

1.75888 

1.76706 

1.77525 

1.78344 

1.79163 

1,79981 

1.80800 

1.81619 

1.82438 

1.832S6 


213 

214 

215 

216 
217 
21S 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 
23B 

239 

240 

241 

242 

243 

244 

245 

246 

247 
246 
249 
2$9 

251 

252 

253 

254 


1.84075 
1.84894 
1.85713 
1,86531 
1,87350 
1.88169 
i. 88988 
1.89806 
1.90625 
1.91444 
1.92263 
1.93081 
1.93900 
1.94719 
1.95538 
1.96356 
1.97175 
1.97994 
1.98813 
1.99631 
2.00450 
2.01269 
2.02087 
2.02906 
2.03725 
2.04S44 
2.0S363 
2.06181 
2.07000 
2.07819 
2.08637 
2.09456 
2.I027S 
2.11094 
2.11913 
2.12731 
2.13550 
2.14369 
2.15188 
2.16006 
2.16825 
2.17644 
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:55 

255 
25:^ 

256 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
232 

283 

284 
26S 
286 
287 
286 
209 

290 

291 

292 

293 

294 

295 

296 

297 
29S 

299 

300 

301 

302 

303 

304 
30$ 

306 

307 

308 

309 

310 

311 

312 


31S 


2.18463 

2.19281 

2.20100 

2.20919 

2.21738 


2.22556 

2.23375 

2.24194 

2.25013 

2.25831 

2.26650 

2.27469 

2.28288 

2.29106 

2.29925 

2.30744 

2.31563 

2.32381 

2.33200 

2.34019 

2.34838 

2.3S656 

2.36475 

2.37294 

2.38112 

2.38931 

2.397S0 

2.40S69 

2.41388 

2.42206 

2. 4302$ 

2.43844 

2.44662 

2.45481 

2.46300 

2.47119 

2.47938 

2.48756 

2.49575 

2.50394 

2.51213 

2.52031 

2.52850 

2.53669 

2.54488 

2.55306 

2.56125 

2.56944 

2.57763 

2.58561 

2.59400 

2.80219 

2.81038 

2.81856 

2.82675 

2.63494 

2.64313 

2.65131 

2.6S9S0 

2.66769 

2.67580 
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316 

2.68406 

317 

2.69225 

318 

2.70044 

319 

2.70863 

320 

2.71681 

321 

2.72500 

322 

2.73319 

323 

2.74138 

324 

2.74956 

325 

2.75775 

326 

2.76594 

327 

2.77413 

328 

2.78231 

329 

2.79050 

330 

2.79869 

331 

2.80687 

332 

2.81506 

333 

2.82325 

3  34 

2.83144 

335 

2.83963 

336 

2.84781 

337 

2.85600 

338 

2.86419 

339 

2.87238 

340 

2.880  56 

341 

2.88875 

342 

2.89694 

343 

2.90513 

344 

2.91331 

34S 

2.92150 

346 

2.92969 

347 

2.93788 

348 

2.94606 

349 

2.95425 

350 

2.96244 

351 

2.97063 

352 

2.97881 

353 

2.98700 

354 

2.99519 

355 

3.00338 

356 

3.01156 

357 

3.01975 

356 

3.02794 

359 

3.03613 

360 

3.04431 

361 

3.05250 

362 

3.06069 

363 

3.06886 

364 

3.07706 

365 

3.08525 

366 

3.09344 

367 

3.10163 

366 

3.10981 

369 

3.11800 

370 

3.12619 

371 

3.13438 

373 

3.14256 

373 

3.15075 

374 

3.15894 

375 

3.16713 

376 

3.17531 
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377 


3.18350 


378  3.1916S 

379  3.19988 


380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 

394 

395 


3.20806 

3.21625 

3.22444 

3.23263 

3.24081 

3.24900 

3.25719 

3.26538 

3.27356 

3.28175 

3.28994 

3.29813 

3.30631 

3.31450 

3.32269 

3.33088 


396 

397 

398 

399 

400 

401 

402 

403 

404 

405 

406 

4M 

98 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

436 

437 


3.33906 

3.34725 

3.35544 

3.36363 

3.37181 

3.38000 

0.100432E-05 

0.216730E-05 

O.335033E-0S 

0.455367B-05 

0.S77763E-05 

0.702248B-05 

0.8288S2E-05 

0.957605S-OS 

0.100B54E-04 

0.122168E-04 

0.13S706E-04 

0.149471E-04 

0.163466E-04 

0.177694E-04 

0.192160E-04 

0.30686SE-04 

0.22iai3e-04 

0.237007E'‘04 

0.252452E-04 

O.2681S0B-04 

0.2l4104e>04 

0.3003196-04 

0.3167976-04 

0.3335436-04 

0.3505606-04 

0.3678516-04 

0.3854216-04 

0.4032736-04 

0.4214116-04 

0.4398396-04 

0.4585606-04 

0.4775796-04 

0.4969006-04 

0.5165276-04 

0.5364636-04 

0.5567lflE-04 
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1 


I 


43a 

43? 

4  40 

441 

442 

443 

444 

445 

446 

447 

448 

449 
4  50 

451 

452 

453 

454 

455 
4  56 
4  57 
45S 

459 

460 

461 

462 

463 

464 

465 

466 

467 
44“  R 

469 

470 
4Tl 

472 

473 

474 

475 
4-6 
477 

471 
479 
4  8^3 

481 

482 
464 
464 

4  as 

486 

487 

488 

489 

490 

491 
493 

493 

494 

495 

496 

497 
4«C 


0.577283E-04 

0.598174E-04 

0.619392E-04 

0.640942E-04 

O.662827E-04 

0.685053E-04 

0.707623E-04 

0.730542E-04 

0.753816E-04 

O.777449E-04 

0.80144SE-04 

0.825810B-04 

0.850548E-04 

0.875665B-04 

0.901165n-04 

0.927055E-04 

0.9533386-04 

0.980021E-04 

0.100711E-03 

0.103461E-03 

0.106252E-Q3 

0.1090866-03 

0.1119626-03 

0.114881S-03 

0.117845B-03 

0,120053  ii-OS 

0.123906B-03 

0.127005B-03 

0.130150E-03 

0.133342B-03 

0.136582G-03 

0.139B70B-03 

0,1432078-03 

0.,\46594E-{n 

0.150031E-03 

0.153S19E-03 

0.1570S9e-03 

0.160692B-03 

0.1 64 29 08-03 

0. 1679988-03 

0.1717536*03 

0.175564B-&3 

0.179430E-03 

0.183 35 5E-03 

C.1I7337E-03 

0.19137IE-03 

0.195480B-03 

0.1996412-03 

0.303865E-03 

0.200151E-OJ 

0.312501E-03 

0.2169156-03 

0.22139SE-03 

0.225941E-03 

0.2305556-03 

0.1352378-03 

0.3399896-03 

0./448? 36-03 

0, 2497066-03 

0.2546746-03 

0.2597166-03 
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■;?9 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 

513 

514 

515 

516 

517 
5J3 

519 

520 

521 

522 

523 

524 

525 
$2$ 

527 

528 

529 

530 

531 

532 

533 

534 

535 
S3« 

537 

538 

539 

540 

541 

542 

543 

544 

545 

546 

547 

548 

549 

550 

551 

552 

553 

554 

555 

556 

557 

558 

559 


0.264834E--03 

0.27002eE-03 

0.275300E-03 

0.280652E-03 

0.286084E-03 

0.291598E-03 

0.297196E-03 

0.302878E-03 

0.308647C-03 

0.314504C-03 

0.320450E-03 

0.32648BE>03 

0.332617E-03 

0.338842E-03 

0.34S162E-U3 

0.3S3S80E-03 

0.3S8098E’03 

C.364718E-03 

0.371441E-03 

0.378269E-03 

''.3E5205E-03 

0.392251E-Q3 

0.399409E-03 


0.406680E>03 

0, 4140688-03 

0.421575E-03 

0.439203B-03 

0.4369546-03 

0.4448316-03 

0.4526386-03 

0,4609766-03 

0.4692466-03 

0.4776586-03 

0.48630BE-03 

O.t *490 16-03 

0.5037416-03 

0.5127316-03 

0.5218746-03 

0.5311746-03 

0.5406346-03 

0.5502586-03 

0.5600496-03 

0.5700136-03 

0.5801536-03 

0.5904726-03 

0.6009766-03 

0.6116696-03 

0.6225556-03 

0.6336406-03 

0.6449286-03 

0.6564246-03 

0.660)346-01 

0.6600636-03 

0.6922176-03 

0.7046016-03 

0.7172216-03 

0.7300646-03 

0.7431966-03 


0.7565636-03 

0.7701936-03 


0.7848936-0) 
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560 

561 

562 
56  3 

564 

565 

566 

567 

568 

569 

570 

571 

572 

573 

574 

575 

576 

577 

578 
579! 
580 

581 

582 

583 

584 

585 
SS6 

587 

588 
$89 
59Q 

591 

592 

593 

594 

595 
59« 

597 

598 

599 
«00 
601 
602 

603 

604 

605 

606 
607 
601 

609 

610 
611 
613 

613 

614 

615 

616 
617 
611 

619 

620 


0.79C268E-03 

0.912729E-03 

0.827481E-03 

0.a<2S33E-03 

0.857894E-03 

O.873571E-03 

0.889574E-03 

0.905912E-03 

0.922595E-03 

0.939631E-03 

Q.9S7032E-03 

0.974a08E-03 

0.992969E-03 

0.30n53E-02 

0.103049E-02 

0.104988B-02 

0.306970E-02 

0.108996E-02 

0.111060E'02 

0.il3187E-0: 

0.U53S5E-02 

Q.U7573B-Q2 

O.U9a43E-02 

0.X22165E-02 

0.U4543F.-02 

0.126976E-02 

0.129468E-02 

0.132019E-02 

0.136633E>Q2 

0.1373098-02 

0.140Q51E-02 

0.142861E-02 

0.145740E-02 

0.148692E-02 

0.1S1717E-02 

0.1548188-02 

0.157998B-02 

0.161259E-02 

0.164604E-02 

0.1680356-02 

0.1T1S54E-02 

O.mi666-02 

O.17««t36-02 

0.ie2676£-02 

0.1865816-02 

0.1905906-02 

0.1967066-02 

0.1989336-02 

0.2032736-02 

0.2077316-03 

0.3133106-03 

0.3170156-03 

0.2318506-03 

0.3368176-03 

0.3319336-03 

0.3371706-02 

0^3625656-03 

0.3481 106-03 


0.2S3 8136-02 


0.3596756-02 


0.2657056-83 
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021 

622 

623 

624 

625 

626 

627 

628 

629 

630 

631 
6  32 

633 

634 

635 

636 

637 

638 
635 

640 

641 

642 

643 

644 
64  5 

646 

647 

648 

649 

650 

651 

652 
65) 
6*4 

655 

656 

657 

658 
651 
660 
661 
662 

663 

664 

665 

666 
667 
661 

669 

670 

671 
67J 
67) 

674 

675 

676 

677 
671 

679 

680 
6«1 


0.271906E-02 

0.278284F.-02 

0.284846E-02 

0.291S96E-02 

0.298541E-02 

0.305687E-02 

0.313041E-02 

0.32060aE-02 

O.320398B-O2 

0.33641SE-02 

0.34466aE-02 

0.353164E-92 

0.361911E-02 

0.370918E-02 

0.38ni92E'02 

0.3B9742E“02 

0.399S77E-02 

0.409706E-02 

0.420129E-02 

0.430685E-C1 

0.441955E-02 

0.4S33S9E-02 

0.465108C-02 

0,477213E“02 

0.48968SE-02 

O.SO2S37E-02 

0.51S780E-02 

0.529427E-02 

0.S43492E-02 

0.5579876-02 

0.5729276-02 

0.588)266-02 

0.6041995-02 

0.62 05626-02 

0.6)74296-02 

0.6548186-02 

0.6727456-02 

0.6912296-02 


0.7102866-02 

0.7299)56-02 

0.7501976-02 

0.7710906-02 

0.7926)66-02 

0.8148556-02 

0.8)77696-02 

0.8614026-03 

0.8657766-02 

0.8109166-02 

0.9)68486-02 

0.96)SK6-02 

0.9911886-02 

0.1019656-01 

0.1049016-01 

0.1079)16-01 

0.1110566-01 

0.1142816-01 

0.1176016-01 

0.1210416-01 


0.1243846-01 

0.1283)96-01 


0.1)50126-01 
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682 

0.135905E-01 

6S3 

0.139923E-01 

684 

0,144070E-01 

685 

0.148351E-01 

686 

0.152769E-01 

687 

0.157330E-01 

688 

0.162038E-01 

689 

0.166899E-01 

690 

0.171917E-01 

691 

0.177099E-01 

692 

0.182449E-01 

693 

0.187974E-O1 

694 

0.193679E-01 

695 

0.199570E-01 

696 

0.205656E-01 

697 

0.2119  41E-01 

698 

0.218433E-01 

699 

0.225140F-01 

700 

0.232068E-01 

701 

0.239227E-01 

702 

0.246623E-01 

703 

0.254267E-01 

704 

0.262165E-01 

705 

0.270329E-01 

706 

0.278767E-01 

707 

0.287489E-01 

708 

0.296506E-01 

709 

0.305829E-01 

710 

0.315469E-01 

711 

0.32S438E-01 

712 

0.335748E-01 

713 

0.346412E-01 

714 

0.357443E-01 

715 

0.368856E-01 

716 

0.380665E-01 

717 

0.392885E-01 

713 

0.405533E-01 

719 

0.418624E-01 

720 

0.432177E-01 

721 

0.446209E-01 

722 

0.460739E-01 

723 

0.475787E-01 

724 

0,491375E-01 

725 

0.507523E-01 

726 

0.524255E-01 

121 

0.541594E-01 

728 

0.559565E-01 

729 

0.578195E-01 

730 

0.597510E-01 

731 

0.617540E-01 

732 

0.638314E-01 

733 

0.659864n-01 

734 

0.682223B~01 

735 

0.705426B-01 

736 

0.729509E-01 

737 

0.7545.10E-01 

738 

0.780469E-01 

739 

0.a07429B-01 

740 

0.835433E-01 

741 

0.864528E-01 

742 

0.894763B-01 
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743 

0.926189E-01 

744 

':-.9  58  860E-01 

745 

0.992332E-01 

746 

0.102816 

747 

0.106492 

•:’48 

0.110317 

749 

0.114297 

750 

0.118441 

751 

0.122755 

752 

0.127248 

753 

0.131929 

754 

0.136806 

755 

0.141889 

756 

0.147187 

757 

0.152711 

758 

0.158473 

759 

0.164483 

760 

0.170755 

761 

0.177300 

7  62 

0.184133 

763 

0.191267 

7  64 

0.198718 

765 

0.206502 

766 

0.214636 

767 

0.223137 

768 

0.232024 

769 

0.241316 

770 

0.251035 

771 

0.261202 

772 

0.271840 

773 

0.282975 

774 

0.294631 

775 

0.306835 

776 

0.319618 

777 

0.333007 

778 

0.347037 

775 

0.361740 

780 

0.377152 

781 

0.393310 

782 

0.410255 

783 

0.428028 

784 

0.4  46  6  73 

785 

0.466237 

786 

0.486770 

787 

0.508323 

788 

0.530952 

789 

0.554714 

790 

0.579671 

791 

0.605887 

792 

0.633431 

793 

0.662376 

794 

0.692797 

795 

0.724775 

796 

0.758395 

797 

0.793747 

798 

0.830926 

799 

0.870032 

800 

0.911171 

801 

0.954454 

802 

1.00000 
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DIMENSION  XYKSO)  ,Xy2(810)  ,TTL4(5)  ,NOE(2)  ,IPTTL(4) 

DIMENSION  IPCRV(2) ,KTYPE(2) ,TTL3(4) ,MrTL(4) ,TSIZE(4) ,RYTTL(4) 
BiTE  XTL(40)  ,yTL(40)  ,TTL(240)  ,DNAME(130)  ,PNAME(13) 

BYTE  YES,NO,ANS,SYMB 

BYTE  RED , GREEU , BLACK , PSCL , PTL , PCRV , ESC , HOME , CLR 

COMMON/GRPH/YTLRX,XTLRY 

DATA  YES/'YV 

DATA  NO/'NV 

DATA  RED/  'RV 

DATA  GREEN/ 'G'/ 

DATA  BLACK/ 'L  / 

DATA  TTL3/'-  AC','TOAL','  DAT', 'A'/ 

DATA  Ti%4/'«  ES', 'TIMA', 'TED  ' , 'VALU' , 'ES '/ 

ESC-27 
HOME-72 
CLRt«74 
SYMD* 1 
ISCL-0 

CALL  ERRSBT(64, . JRUE. , .TRUE. , .TRDE. , .FALSE. ,3 0) 

CALL  ERRSET  (29,. TRUE.,  .TRUE.  ,  .TR'IE. ,  .FALSE.  ,30) 

CALL  ERRSET(43,.TR0S., .TRUE., .FALSE., .FALSE, ,30) 

CALL  EHRSET(S9, .TRUE . , .FALSE. , .TRUL . , .FALSE. , 3 0) 

CALL  ERRSETOO,. TRUE.,  .FALSE.,  .TRUE.,  .FALSE. ,30) 

11  COmiNUE 

WRITE(5,S55)  ESC, HOME, BSC, CLR 
WRITE(5,415) 

1  READ(5,405,ERR-2)  IRW 
IP(IRW.EQ.1.0R.IRW.EQ.2.0R.IRW.EQ.3)  GO  .0  3 

2  WRITE(5,1S5; 

GO  TO  1 

3  IF(IRW.EQ.2.0R.IRW.EQ.3)  IRLUN-5 
IFdRW.Ey.l)  IRLON-1 

WRITE(5  65) 

10  'rSAD(5,55)  LENP,(PNAME(I) ,1-1,12) 

IF(LENP.GT.12)  GO  TO  20 
PNAME  (LENP-t-l) -0 
CALL  ERRrST(43,K) 

CALL  ASSIGN (2, PNAMB) 

CALL  ASSiaN(4, 'GHAPH1.TXT') 

IP(1RW.NE.3)  CALL  ASSIGN (1 , 'GRAPHl .COM ' ) 

CALL  ERRTST(43,K) 

IF(K.EQ.2)  GO  TO  21 

20  WRITE(S,115) 

GO  TO  10 

21  KRTTPd  SKI 

WRITE(4,305)  <P^ME(I)  ,I-1,LENP) 

30  ISCL-ISCL+l 

IXMM-0 
IYMM-O 

WRITE<5,S55)  ESC, ROME, ESC, CLR 
NDIQ-2 

IFdSCL.LT.lO)  NDlQ-1 

IP(IRW,EQ,2.0R,IRW.EQ.3)  WRITE(5,355)  ISCL 
IF(IRW.EQ.2.0R.IRW.SQ.3)  WRITE(5,463) 

40  READ(IRLUN,10£)  LEN,ANS 
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IF(LEN.EQ.O)  AN£-'N' 

IF(AK3.EQ.yES.0R.ANS.EQ.N0)  GO  TO  50 
WRITE(5,155) 

GO  TO  40 

50  IP(ANS.EQ.XES)  IXMM-1 

IF(IXMM.EQ.O)  GO  TO  110 
IF(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,485) 
60  READ(IRLUN,*,ERR-70)  XMIN 

GO  TO  80 

70  WRITE(5,155) 

GO  TO  60 


80 

90 


GO  TO  110 


IF<IRW.EQ.2.0R.1RW.EQ.3)  WRITE(5,495) 
READ(IRLUN,*,ERR-100)  XMAX 


100 


WRITE(5,155) 


GO  TO  90 

110  WRITE(4,245) 

WRITE(4,355)  ISCI. 

WRITB(4,46S) 

WRITB(4,26S)  ANS 
IP(IRW.EQ.2)  WRITE<1»585)  ANS 
IP (ANS.EQ.no)  go  TO  111 
WRITB(4,485) 

WRITE(4,305)  XMIN 
IP(IRW.EQ.2)  WRITE(1,435)  XMIN 
VIRITE(4,495) 

WRITE(4,385)  XMAX 
IP(IRW.EQ.2)  WRITB(1,435)  XMAX 

111  IP(IRW.BQ.2.0R.1RW.EQ.3)  WRITE(5,475) 

120  READ (IRLON,  105)  LBN, ANS 

IF(LBN.EQ.O)  ANS"'N' 
IPCANS.EQ.YBS.OR.ANS.BQ.NO)  go  to  130 
WRITE(5,1S5) 

GO  TO  120 

130  IP(ANS.EQ.YES)  IYMM-1 

IPdVMM.EQ.O)  GO  TO  170  . 

IP<lRW.Ea.2.0R.IRW.EQ.3)  WRITB(5,485) 

140  RBAD(IRDUM,*,BRR»150)  .YMIN 

GO  TO  160 

150  WRITE(5,155) 

GO  TO  140 

ISO  IP<IRW.EQ.2.0R.IRW.BQ.3)  WRlT£i(5,495) 

READdRDUN,*)  YMAX 

170  WR1TB(4,475) 

WRITB(4,26S)  ANS 
IFdRW.BQ.2)  WRITE(1,585)  ANS 
IF (ANS. BQ. NO)  GO  TO  171 
WRITE(4,4e5) 

WRITE (4, 385)  YMIN 
IPdRW.BQ.2)  WRITE(1,435)  YMIN 
WR1TB(4,495) 

HRITE(4,385)  YMAX 

lPdRW.EQ.2)  WRITB{1,435)  YMAX  . 

171  IPdRW,EQ.2.0R.IRW,EQ.3)  WR1TB(5,555) 
lPdRW.BQ.2.0R.ISW.EQ.3)  WR1TB(5,125) 


ESC  >  HOME,  ESC.  cut 
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350  IF(ANS.EQ.YES)  ISCX-1 

WRITE(4,245) 

WRITE(4,355)  ISCL 
WRITE(4,2S5) 

WRITE(4,195) 

WRITE (4, 275)  XLEN 
IF(IRW.EQ.2)  WRITE(1,425>  XLEN 
WRITE(4,205) 

WRITE(4,275)-  XBLINT 
IF(IRW.EQ.2)  WRITE(1,42S)  XBLINT 
WRITE(4,185) 

WRITE(4,265]  ANS 
IP(IRW.EQ.2>  WRITE(1,585)  ANS 

IF(IRW.EQ.2.0R.IRW.EQ.3)  WRITE{5,555)  ESC , HOME , ESC , CLR 
IF(IRW,EQ.2.0R.IRW.EQ.3)  WRITE(5,355)  ISCL 
IF(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,375> 
IF(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,215) 

360  READ (IRLUN,*,ERI^ 370)  YLEN 

IF(YLEN.GE.l..AMO.yLEN.LE.10.)  GO  TO  380 
370  WRITE(5,1S5) 

GO  TO  360 

380  IP(IRW.EQ.2.0R.IRW.EQ,3)  WRITE(5,205) 

390  REAO<IRLUN,*,ERR«400)  YBLINT 

IF<YBLINT.LE.0..OR,YBLINT.GT,yLEN)  GO  TO  400 
GO  TO  410 

400  WRITB(S,155) 

GO  TO  390 

410  IP<IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,185) 

420  READ (IRLUN, 105}  L£N,ANS 

IP(LBN.EQ.O)  ANS-'N' 

IP (ANS.BQ.YES.OR.ANS.EQ.no)  GO  TO  430 
WRITE(S,155) 

GO  TO  420 

430  IP (ANS. EQ. YES)  ISCY-1 

WRITE(4,24S) 

WRITE(4,355}  ISCL 
WRITB(4,375) 

HRITE(4,215) 

WRITE (4, 27 5)  YLEN 
IP(IRW.EQ.2)  WRITB(1,42S)  YLEN 
WRITE(4,20S) 

WRITE(4,275)  YBLINT 
IP(IRW.EQ.2}  WRITB(1,42S)  YBLINT 
WRITE(4,185) 

WRITE(4,265)  AMS 
IP(IRW.EQ.2}  WRITE(1,585)  ANS 

IP(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(S,555)  ESC, HOME, ESC, CLR 
IP(IRH.EQ.2.0R.IRW.EQ.3)  WRITE(5,355)  ISCL 
IP(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,95) 
IP(IRW.EQ.2.0R.IRW.EQ.3}  NRITE(5,225} 

440  RBAO(IRLUN,455,ERR-450}  LEN,IPSCL 

IP(LEM.EQ.0>  IP8CL-1 

IP(IPSCL.GT.0.AMD.IPSCL.LT.4)  GO  TO  460 
450  NRZTB(5,155) 

GO  TO  440 


) 
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IF(IRW.EQ.2.0R.IRW.EQ.3)  TOITE(S,365) 

180  READ(IRLON,4  05,ERR»19  0)  IC 

IF(IC.EQ.1.0R.1C.EQ.2.0R.IC.EQ.3)  GO  TO  200 
190  WRITE(5,155) 

GO  TO  180 

200  IP(IC.NE.l)  GO  TO  210 

NTYPE-0 
GO  TO  250 

210  IP(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,395) 

220  READ(IR]:.UN,4  05,ERR«23  0)  NCTCLE 

IP(NCYCLE.GE.I)  go  TO  240 
230  VJRITE(5,155) 

go  to  220 

240  IF(IC.EQ.2)  NTYPE-NCYCLE 

IP(IC.EQ.3)  NTYPE— NCYCLE 

250  WRITE(4,245) 

WRITE(4,125) 

TOITE(4,365) 

WRITE(4,285)  IC 
IP<IRW.EQ.2)  WRITE(1,405)  IC 
IP(IC.EQ.l)  GO  TO  251 
WRITE(4,395) 

WRITE<4,285)  NCYCLB 
IPdRW.EQ.a)  WRITEa,405)  NCYCLE 

251  IF<IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,165) 

260  ■  READ (IRLUN, 105)  LEN,AMS 

IF(DEN.EQ.O)  ANS-'N' 

IF(ANS. EQ.YES.OR.ANS.EQ.no)  GO  TO  270 
WRITE(5,155) 

GO  TO  260 

270  IP<ANS.EQ.yES)  IRT-5 

IP(ANS.BQ.NO)  IRT-1 
WRITE (4, 165) 

WR1TE(4,265)  ANS 
IP(IRW.EQ.2)  WRITE(1,585)  ANS 

IF(1RW.EQ.2.0R.IRW.EQ.3)  VIRITECS , 555)  ESC, HONE, ESC, CLR 
IP(IRW.EQ.2.0R.1RW.EQ.3)  WRITE(5,355)  ISCL 
IF(IRM.EQ.2.0R.IRW.EQ.3)  WRITE(5,255) 
IP<1RW.EQ.2.0R.IRW.EQ.3)  WRITB(5,195) 

280  READ < IRLUN, *,ERR-290)  XLEN 

IF(XLEN.GE.l.)  GO  TO  300 
290  WRITE(5,155) 

GO  TO  280 

300  IP(1RW.EQ.2.0R.IRW.EQ.3)  WRITE(5,205) 

310  READ ( IRLUN, *,ERR-320)  XBLiNT 

IF(XBLINT.LE.0..OR.XBLINT.GT.XLEN)  GO  TO  320 
GO  TO  330 

320  WRITE<5,155) 

GO  TO  310 

330  IF<IRW.EQ.2.0R.IRW.EU.3)  WRITE(5,185) 

340  READ (IRLUN, 105)  LEN,ANS 

IF(LEN.EQ.O)  ANS>'N' 

IF (ANS.EQ.YES.OR.ANS.EQ.no)  go  to  350 
WRITB(S,155) 

GO  TO  340 
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460  IF(IRW.EQ.2.0R.IRW,EQ.3)  WRITE(5,565) 

470  READ (1RLUN> 575, ERR-480)  LEN,ASIZE 

IF(LEN.EQ.O)  ASIZE-.l 

IF(ASIZE.I.T..07.0R.ASI2E.GT..5)  GO  TO  480 
GO  TO  490 

480  WRITE(5,155) 

GO  TO  470 

490  WRITE<4,245) 

WRITE (4, 35 5)  ISCL 
WRITE(4,95) 

WRITE(4,225) 

WRITE <4, 285)  IPSCL 
IF(IRW.EQ.2)  WRITE(1,405)  IPSCL 
WRITE(4,565) 

WRITE(4,275)  ASIZE 
IP<IRW.EQ.2)  WRITE(1,595)  ASIZE 

IF(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,555)  ESC , HOME, ESC, CLR 
IF(IRW.EQ.2.0R.ISW.EQ.3)  WRITE(5,355)  ISCL 
IF(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,445) 
IF(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,315) 

500  READ  (IRLON,  325)  )1XTL,XTL 

IF(NXTL.GT.40)  GO  TO  510 
IP(IRW.EQ.2.AND.NXTL.EQ.O)  WRITE(1,585) 

IP(NXTL.EQ.O)  GO  TO  589  ’ 

GO  TO  520 

510  IP(IRW,EQ.2.0R.1RW.EQ.3)  WRITE(5,1S) 

GO  TO  500 

520  IP(IRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,565) 

530  RBAD(IRLCN,57S,ERR-540)  LEM,XSIZE 

IF(LEK.EQ.O)  XSIZE-.2 

IF(XS1ZE.LT..07.0R.XSIZE.GT..5)  GO  TO  540 
GO  TO  550 

540  WRITE(5,155) 

GO  TO  530 

550  IF(IRW.EQ.2.0R.IRH.EQ.3)  WRITB(5,225) 

560  READ(IRLUN,455,ERR-57  0)  L£l),IPXTL 

IF(LEN.EQ.O)  IFXTL-1 
IFdPXTL.GT.O.AMO. IPXTL.lt. 4)  GO  TO  580 
570  WRITE(5,155) 

GO  TO  560 

580  '  WRITE(4,245) 

WRITE(4,445) 

WRTTPJA 

WRITE(4)305)  (XTL(I) ,I-1,NXTL) 

IF<IRW.EQ,2)  WP.ITE(1,585)  (XTL (I)  , >1  ,NXTL) 

WRITE(4,S65) 

WRITE<4,275)  XSIZE 
IF(IRW,EQ.2)  WR1TE(1,595)  XSIZE 
WRITE(4,225) 

WRZTE(4,285)  IFXTL 
IF(IRW.EQ.2)  WRITE(1,405)  IPXTL 

589  IF(IRH.EQ.2.0R.IRW.EQ.3)  WRITE(S,335) 

590  READ  (IRLOM,  325)  irrFL,m 
IF(NyiL.QT.40)  GO  TO  600 
IF(IRW.EQ.2.AND.Mm.EQ.0)  WRITE(1,585) 
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IP(Nm.EQ.O)  GO  TO  676 
GO  TO  610 

600  IF(IRW.EQ..2.0R.IBW.EQ.3)  WRITE(5,15) 

GO  TO  590 

610  IF(1RW.EQ.2.0R.IRW.EQ.3)  WRITE(5,565) 

620  READ (IRLON, 575, ERR-63  0)  LEN,YSIZE 

IP(LEN.EQ.O)  VSIZE-.2 

IP(YSIZE.LT..07.0R.YSI2E.GT..S)  GOTO  630 
GO  TO  640 

630  WRITE(5,155) 

GO  TO  620 

640  IP(1RW.EQ.2.0R.IRW.EQ.3)  WRITE(5,225) 

650  READ(IRLUN,455,ERR«66  0)  LEN,  IPYTL 

IF(LEN.EQ.O)  IPOTL*1 

IPdPYTD.GT.O.AND. IPYT1.lt. 4)  GO  TO  670 
660  WRITE(5,15S) 

GO  TO  650 

670  WRITE(4,335) 

WRITE(4,305)  (YTL(I) ,I«1,NYTL) 

IF(IRW.EQ.2)  WRITE{1,585)  (YTL(I)  ,I-1,NYTL) 

WRITE(4,565) 

WRITE(4,275)  YSIZE 
IPdRW.EQ.a)  WRITE(1,595)  YSIZE 
WRITE(4,225) 

WR1TE(4,28S)  IPYTL 
IPdRW.EQ.2)  WRITE(1,405)  ipm 

676  lPdRW.EQ.2.0R.lRW.EQ.3)  WRITE(5,555)  ESC, HOME, ESC, CIR 

WRITE(4,24S) 

DO  780  1-1,4 

IPdRW.EQ.2.0R.IRW.EQ.3)  MRITB(5,345)  I 

ISTRT-d-l)*60+l 

lEND-ISTRT+59 

680  READ (IRLON, 325)  NTTL (1) , (TTL(J) ,J-ISTRT,IEND) 

IF(NTTLd)  .GT.50)  GO  TO  690 
lFdRW.EQ.2.AND.NTn(I).EQ.O)  WRITE(1,S85) 

IP(NTTLd)  .EQ.O)  GO  TO  780 
GO  TO  720 

690  IPdRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,15) 

GO  TO  680 

720  IPdRW.EQ.2.0R.lRW.EQ.3)  WRITE(5,235) 

730  READ  (IRLON, 575, BRR-74  0)  LEN,TSIZEd) 

IP(LEN.EQ.O)  TSIZE(l)-.2 

IF(TSIZEd)  .LT..07.0R.TSI2B(I).GT..5)  GO  TO  740 
GO  TO  750 

740  WRITE (5, IS S) 

GO  TO  730 

750  IPdRW.EQ.2.0R.IRW.EQ.3)  WRITE(5,175) 

760  READdRLUN,455,SR»-77  0)  LEN,IPTTL(I> 

IF<LEN.EQ.O)  IPTTL(I)-l 

IPdPTTLd)  .QT.O.AND.IWTL(I)  .LT.4)  GO  TO  771 

770  WRITE(5,155) 

GO  TO  760 

771  WRtTB(4,345)  I 
JBND-ISTRT+NTTL  <1 ) -I 

IP(WrTLd)  .GT.O)  WRITE(4  ,  295)  (TTL(J)  ,J-1STRT,JEND) 
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IF(IRW.EQ.2)  WRITE(1,S85)  (TTL( J) ,J*ISTRT, JEND) 
WRITE{4,235) 

WRITE(4,275)  TSIZE(I) 

IF(IRW.EQ.2)  WRITE(1,S95)  TSIZE(I) 

WRITE(4,175) 

WRITE(4,285)  IPTTL(I) 

IF(IRW.EQ.2)  WRITE(1,405)  IFTTLC) 

780  CXJNTINOE 

WRITE(5,555)  ESC, HOME, ESC,CIiR 

NDIGl-2 

NDIG2-2 

IFdSCL.LT.lO)  NDIGl-1 
IP(ICRV.LT.IO)  ND1G2-1 
WRITE<5,3S5)  ISCL 
WRITE(5,45)  ISCL 
IPNS-(ISCL-l)  *12+1 
IPNE-ISCL*12 

790  READ(S,S5)  LEND, (DNRME(I) ,1-IFMS,IFNE) 

LENDE-IPNS+LEND 
DNAME(LENDE)-0 
CALL  BRRTST(29,X) 

OPEN  ( UNIT-3, NAME-DNAMEdFNS)  , READONLY, TYPE-' OLD ', ERIN 800) 
CALL  ERRrST(29,K) 

IP<K.EQ.2)  GO  TO  810 

800  WRITE(5,155) 

GO  TO  790 

810  WRITB(4,245) 

WRITB<4,35S)  ISCL 
WRITB(4,4S)  ISCL 

HRITE(4,30S)  (DNAMEd)  ,I-IFNS,IFNE) 

DO  820  1-1,24 
RBAD<3,S5)  LEN,T 

820  CONTINUE 

NOE(1)-0 

I-l 

830  RBAD(3,5)  LEN,XYld) 

IF(LEN.EQ.O)  GO  TO  840 
NOB(l)-NOB(l)+l 
l-I+l 
GO  TO  830 

840  NOE(l)-NOE(l)-l 

IF(NOB(l) .LE.IOOO)  GO  TO  850 
IFdRW.EQ.2.0R.IRW.EQ.3)  WRITB(5,75) 

GO  TO  9999 

850  DO  860  1-1,11 

READ (3, 55)  LBN,T 

860  CONTINUE 

NOE(2)-l 
I-l 

870  READ (3, 5, 0)0-880)  L£N,XY2d) 

NOE(2)-NOE(2)+l 

I-I+l 

GO  TO  870 

880  •  NOE(2)-NOB(2)-l 

REWIND  3 
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DO  890  1-1,24 
READ (3, 55)  LEN,T 
890  CONTINUE 

DO  900  I»l,NOE(l) 

READ(3,*)  LEN.XYKI) 

900  CONTINUE 

DO  910  1-1,12 
READ (3, 55)  LEN,T 
910  CONTINUE 

DO  920  I-l,NOE(2) 

READ(3,*)  I,EN,XY2(I) 

920  CONTINUE 

CLOSE  (UNIT-3) 

N0E(1)-N0E(1)/2 

NOE<2)-NOE(2)/2 

XMINT-1.B27 

YMINT-1.E27 

XNAXT-1  .E-27 

YMAXT-l.E-27 

DO  930  I>l,NOE(l) 

IP(XY1(I).LT.XM1NT)  XMINT-XVKI) 

IP(XY1(I)  .GT.XMAXT)  XNAXT-XYKI) 

IP(Xyi(NOE(l)+I)  .LT.YHINT)  YHINT-XYKNOBd)  ♦!) 
IP(XYl(NOB(l)+I)  .OT.YMWCT)  YNAXI-XYKNOBd) +1) 

930  CONTINUE 

DO  940  I-l,NOE(2) 

IP(XY2{1).LT.XHINT)  XHINT-XY2(I) 

IP{XY2{I).QT.XMAXT)  XHAXT-XY2(I) 

IP(XY2(NOE(2)+I> .LT.YHINT)  YNINT-XY2(NOB  2>*1 
lP(XY2(NOB{2)+I).(IT.YHAJtT>  YNAXT-XY2(NOE{2)  ♦!) 

940  COtrriNUB 

IPdXMN.EQ.l)  GO  TO  950 

XNIN-XNINT 

XMAX-XNAXT 

950  IPdYNN.EQ.l)  GO  TO  960 

YMIN-YNINT 
YNAX-YNAXT 

960  tmXTB(4,245) 

DO  990  ICRV-1,2 

IPdCRV,BQ.l.AH0.dRI».BQ.2.0R.IRN.BQ,3) )  KRXTB(5,13S) 
IPdCRV.BQ.2.AND.(XRW.BQ.a.0R.IRN.E0.3> )  NRXTB{5,145) 
IP(XRN.EQ.2.0R.XRW.EQ.3)  WR1TB<5,225) 

970  REAO(XRL0N,455.BRi^980)  LBN, XPOlVtXCRV) 

IP(LEN.BQ.O)  XPCRV{XCRV)-l  _ 

IP{XPCRV(XCRV)  .OT.O.AND.IPCRVdCRV)  .LT.4'  GO  TO  981 

980  NRXTE(5,15S) 

GO  TO  970 

981  IF(XCRV.EQ.l)  NRXTE(4,135) 

XP(XCRV,EQ.2)  WRXTB(4.145) 

MRXTE(4,229) 

>IRITE(4,205)  XPCRVdCRV) 

XP(XRH.BQ.2)  NRXTBd,409)  XPOtVCXCRV) 

990  CONTINUE 

IFdRN.BQ.2.0R.IRN.EQ.3)  WRXTE(S,S55)  SSC,tlONE»£8C,CLR 
XQXBB-XBLXHT/.25 


m 
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yQIBB»mrNT/.25 

CALL  BGNSrP(OTYPE/-2, CHIT. 0.1) 
^^(^SCX.EQ.Q)  GO  TO  1000 

10^^^  ®°**T?!”‘XMlN,XMAX,Xt£N,X0IB8) 
„  IFtlSCY.EQ.O)  GO  TO  1010 

att  BOPTMM(YHIN,yMAX,YLEN,y01BB) 

NWPSmiPSCtJ 
DO  1020  JPI-1,18T 
Pa«PWAT(JPI-l) 

TPA«< . 283 1853  E0*Fj/5 . 
AX«.008*COS(TPA) 

AY«.008*SIN(TPA) 
rP(XRT.EQ.l)  AX-0.0 
IPdRT.EQ.l)  A»-0.0 
CAtt  PWT(AX,Ay,-3) 

j/a.+i.s 

Rmti- (YLEH- (RyTL*ysiX8) ) /2  *1 
RnTt{4).yLBH*1.2 

00  1040  r-i,4 
00  1040  JPt«l,iRT 
rj»PU)AT{JPI-l) 

H^*i:i****MM*N/s. 

AX».OOI*COS(TPA) 

AY».00l»*51ll(trA» 
rP(IRT.EQ.l)  AX-0.0 
IFflRT.EO.i)  AJf«0.0 

ilu  Sigi?},gJ''  ““ 

00  1050  jpuittnt 
fO-PlOATtJPt-1) 
m«<.  213 1853  C0«pjr/S. 
W».00l»COS{TfA) 

AY*i.808*SIN<TrA) 

AX-0.0 

traRT.tQ.l)  AY-0.0 
PU}t<AX«A¥»<‘3) 

00  1080  JM-1,2»t 
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FJ»FL0AT(JPI-1) 

TFA-6.2a31853E0*FJ/5. 

AX«.008*COS<TFA) 

AY»,008*SIN(TFA) 

IFdRT.EQ.l)  AX-0.0 
IF(IRT.EQ.l)  AY-0.0 
CALt  PU)T{AX,AY,-3) 

XF(NXTL.EQ.O)  GO  TO  1060 

CALL  BQITTL ( XTL , NXTL , XS2  ZE  >  RXXTL , RYXTL , 1 . , 0 . ) 

1060  COOTINUE 

1061  RXTTL3-XLEN^3. 

CALL  NEWPCNdPXTL) 

DO  1070  JPI-l.lRT 
FJ-PLOAT(JPI-l) 

TPA-6. 2831853 B0*FJ/S. 

AX-.00e*COS(TPA) 

Ay-.00e*SIN(TPA) 

IFdRT.EQ.l)  AX-0.0 
IPdRT.BQ.l)  AY-0.0 
CALL  PU>T(AX.AY.-3) 

CALL  B(»Ttl.(TTL3.1S,.14,lUCm3.4,0,l.,0,> 

CALL  B<»TTL(m4.20..14,RXTTL3.3.7.1.,0.) 

CALL  SYNBOL(i(XTTL3,4.07,.14.1.0.0,>l) 

CALL  SYNBOL(RXTTL3.3.77..14.24>0.0«-1) 

1070  CONTIHUe 
CALL  RCNKNdKViVd)) 

DO  IGBQ  JPI-l,tXT 
M«PLOATtJPt-l) 

TPA-6. 213 llS3eO*PJ/S. 

AX».008*COS(TPA) 

AY-.00l*SX»<TPA) 

IFdRT.EQ.l)  AX-0.0 
IFdRT.EQ.l)  AY-0.0 
CALL  FU)T(AX«AY,*3) 

CALL  BPU}Td8CL.NOE<l)»O.XYl.Xn(ltOEU)«l).1.0 
1  l.SV)a,.14) 

CALL  BOVFLTdSCL) 

10«0  COttriNUi 
CALL  MtWtiRd) 

CALL  WWFENdFOtVd)) 

00  1010  JFI-UIRT 
rj-PtOAT(JFt-l) 

TFA-6 .2i31IS3e9«FJ/S. 

AX-.OOBKOSITFA) 

AY-. 001 OSIN ITFA) 

IFdRT.EQ.l)  AX-0.0 
IFdRT.EQ.l)  AY-0.0 
CALL  FLOTtAX«AY«>3) 

CALL  BFU>TdSCL.tl0eU)  ,l.XY2»XY3tltdt(2)d)  .1.0 

1  0«0..0.) 

CALL  BOVFLY<18CL) 
lORO  CONTMOt 
CALL  NIWFtRd) 

CALL  ailDFLT(O) 

tFdRM.fiO.l.Oa.llM.tO.3)  CALL  CLOSt  d) 
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CALL  CLOSE  (2) 

CALL  CLOSE  (4) 

WRITE<5,S55)  ESC, HOME, ESC, CLR 
WRITE (5,35) 

1100  REAO(S,105)  LEM,ANS 

IF<LEN,EQ.O)  ANS-'N' 

IP(AMS>EQ.yES.OR,ANS.£Q.}»)  GO  TO  1110 
WRITE(5,155) 

GO  TO  1100 
1110  COWriNOB 

IFCANS.EQ.TES)  GO  TO  11 
CALL  NBMPfiN(l) 

CALL  BNDPLT(l) 

WRXTE(5,5SS)  ESC,{10t(E,£SC,CLR 
WRITE(S,2S) 

9999  CALL  EXIT 

5  FOimAT(Q,S5X, F14.il) 

IS  POfnATt'  You  auat  enter  60  characters  or  leas.') 

25  F0IWAT(///////////»29X, 'Succaaaful  coimlation. ' ,//////////) 

35  FORMAT (///////////,'  Create  another  graph  tY/NJ?  ',«) 

45  FORMATC/,'  Input  flic  for  Graph  ‘,I<N1)1G>,'  tfileapec)?  *,4) 

55  FORNAT(Q,13Al) 

65  FORHAT(///,‘  Output  file  (filcspcc)?  ',$) 

75  FORHAT(//, *  Four  data  flic  contains  aorc  than  1000  X/Y* 

1  *  eoordlnatcc.*,//, *  Prograa  atop.  No  output  produced.*) 

15  POiMAT(I4) 

95  FORNATC^t  X  and  Y  axes  apcclflcatlont)  *) 

105  P0RNAT(Q«1A1) 

US  PORNATt*  That  flic  docs  not  exist.*) 

125  PORMATt//,*  Plot  specif  lea  tlonst  *> 

U9  FORNAT(///  * ‘Actual  data**  curve.*) 

145  PORNAT(/,'  **CatlBatcd  values'*  curve.*) 

155  FORNATC  Whatl?  *,S) 

165  POIMATC/,*  Ret  race  option  ty/N)7  *,l) 

175  PORNATt  *  Ml  SLACR  t2)  RO)  U)  GRfiCN*,/, 

1  *  Color?  *,i) 

US  PORNATt/,*  Scale  option  tY/N)?  *,*) 

195  PORNATt/,*  X*Axla  length  In  Inches?  *,*) 

205  FOWATt/,*  Dlstanoe  betaeen  blips  in  Inches?  *,S) 

215  PORNATt/, •  y-Aals  length  in  inches?  *,#) 

225  PORNATt/,*  (1)  UACR  (2)  RED  Ui  GJieEH*,,*, 

1  *  Color?  *,») 

US  PORNATt*  Chsraeter  SUE  .‘n  Inchts  t.07>.S)?  *,*) 

245  PORNATt/,*  , 

^  I ■i»w,a»a»6»ea«eeea»<a^ee«ie«e^eei^ea‘eaeeeMeea^<a^aiyasefcee-«ai^aa^*»«a^eaea^  ^ 

2S5  PORNATt** I  X  axis  speclf Icctioosi ' ,/) 

265  PORNATt •♦•,1A1) 

275  POlNAr{'**,PlS.?J 

265  rOfNAtt***,!!!) 

295  POfNAYt  *  ' ,C0A1) 

305  PORNATt ***,40A1> 

315  POlHATt*  X>Axis  Tltla  tl>40  chart. )?  *,S) 

325  PORNATtQ,60Ai) 

335  PORNATt/,’  YHtxia  Title  tl-40  ehcts.)?  ',S) 

345  fOWArt*  Graph  tlUe  (line  *,Xl.'!  11*60  ebara.l?  ') 
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355  FOI«AT('  Graph  niaaber  '.KNDIO.S) 

365  FOBMAT(/f'  tl}  Linear  (21  5'eai-log  (31  Log-log* 

1  ,/,*  Type  of  scale?  '|S) 

375  FOPMAT('+:  Y  axis  specifications;',/) 

385  FOI»!AT'*+',E15.7) 

395  FOmATl/,'  tJunber  of  cycles  OOl?  ',$) 

405  FOBHATdl) 

415  FORMAT! '  (11  Read  responses  froo  GRAFH1.COM 

1  *  and  do  not  updar-e  GRAFH1.COM.  * ,//, 

2  '  (21  Read  responses  ^  'a  keyboard  and',/, 

3  '  update  GRAJBl. CUM. ',//, 

4  '  (31  Read  respons"’  froa  keyboard  a;^',/, 

5  '  00  not  upda.  GRAFHl.COH. ' ,//, '  Option?  *,$) 

425  ?ORMAT(F15.7) 

435  FORMAT(E15.7) 

445  F0H!AT(  *+  Title  s:  cifications: ' ,/) 

455  fo:mat(q,ii) 

4^'  FOBMAK/,'  X-Axis  mirinuB/waximum  overide  option  (Y/Nl?  ',$) 

475  FORMAT!/,'  Y-Axis  ninisuoi/aaxlBuin  overide  option  (Y/Nl  ?  ',$) 

485  FOWATC  Mini  a.®  value?  ',?) 

495  FO'tIAT!  '  Maxinum  vilue?  ',$) 

555  FORMAT!'  ',4Al) 

365  FORMAT!/,'  Character  size  in  inches  (.07-. 51?  *,5) 

575  FORMAT!Q,P5.4) 

585  FORMAT!60A1) 

595  F0IMAT!P5.3) 

END 

SDBROOTINE  BOPTMM(XMIN,XMAX,AXLEN,QIBB) 


C  BOPTMM  will  determine  a  scale  for  plotting  values  between  XMIN  and  XMAX 
C  such  that  each  l/20th  inch  tic  mark  has  an  easily  READ  decimal  value. 

C  If  QIBB  specifies  tic  marks  on  the  axis,  '0.'  must  fall  on  a  blip  even 
C  if  off  the  plot  page.  AXLEN  is  truncated  to  the  nearest  1/2 0th  inch  and 
C  XMIN  and  XMAX  are  changed  to  the  endpoints  of  the  axis  to  be  drawn. 

C 

C  Arguments: 

C 

C  XMIN  »  Minimum  value  in  array  to  be  plotted. 

C 

C  XMAX  •  Maximum  value  in  array  to  be  plotted. 

C 

C  AXLEN  ■  Axis  length  in  inches. 

C 

C  QIBB  B  Number  of  quarter  inches  between  tic  marks  on  axis. 

C 

C  If  QIBB  is  positive  (>01  -  fit  scale  to  XMIN  and  XMAX 

C  using  tic  mark  interval. 

C  If  QIBB  is  negative  (<01  -  permit  a  2.5%  overscale  on 

either  end  using  tic  marks. 

C  _ 

C - 


REAL*8  XL0W,XEI, POWER 


128 


•I-  - 


—APPENDIX  E— 


rXHl  -IXHI  /lo 
FXMAX-FXMAX/l 0. 
lEXP-IEXP+l 
GO  TO  120 

130  IF  <I2,E5.0)  GO  TO  140 

IX-IABS(IXU3»;-10*{IXLOW/10)  ) 

I?  (IX.EQ.O)  GO  TO  140 
IF  (IXLOW.LV.O)  IX-IO-IX 
IL-IXHI-IX 

IF  (FXHWC.GT.IDGO  to  140 

IXHI-IL  i 

IXLOW-IXIjOW-IX 

GO  TO  120  !' 

140  POWER-10.Da**IBXP  ! 

XLOW-IXLOH*POHER 

XHl  -IXHI  ‘POWER  j; 

IF  tPANGB.ET.O.)  GO  TO  ISO  i 

XMIN-XLOW  h 

XMAX-XHI 
RBTORN 

150  XMIN-XHI 

XKM(-XLOH  P 

RETORN 

5  FOBMATC  Invalid  BOPTMM  arguments: • ,//, 

1  '  Minimum  value  •  •,P14.7,/, 

2  '  Maximum  value  •  ',F14.7,/,  ^ 

3  '  Axis  length  -  ',P14.7,/,  \ 

4  '  Blip  intervals  •  •,F14.7,//,'  Program  stop.') 

END 
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12*0 

RANGE- XMAX-XM IN 
IF  (RANGE. NE.O.)  GO  TO  10 

;P(1RH.£Q.2)  WRIT£(S,S)  XMIN,XHAX,AXL£N,QIBB 
CALL  EXIT 

10  R«ABS( RANGE) 

IF  (RANGE. GT.O.)  GO  TO  20 

X-XMIN 

XMIN-XHAX 

XMAX*X 

20  AXL20«INT{AXLEN*20.) 

AXLEN-AXL20/20. 

aNP20>>R/AXL20 

IF  (CIBB.GT.O.)  go  to  30 

X-.025*R 

XMAX-XMAX-X 

XMIN-XMIN+X 

R-XMAX-XMIN 

tJNP20»R/AXL20 

30  TWPBLP«INT(ABS(QIBB>  )*5 

IF  (TWPBLF.EQ.O.)  go  to  40 
IF  (AXL20/TNPBLP.GT.1.0001)  GO  TO  50 
IF  (R.LE.ABS(XMAX+XMIN))  GO  TO  50 
40  TWPBLP-1. 

IZ«1 

50  lEXP— 1 

IF  (ONP20.GE.1.)  IEXP-9 
60  X-0, 

n-10.**IEXP 
70  X-X+1 . 

PONP20-X’»Y1 

IF  (DNP20.GT.PONP20)  GO  TO  70 
IF  (X.NE.l.)  GO  TO  80 
IP  (ONP20.EQ.PUNP20)  GO  TO  80 
IEXP-IEXP-1 
GO  TO  60 
80  IX-X 

IAXL20-AXL20 
ITWPBL-TWPBLP 
PXMAX-XMAX/Yl 

90  IUNPBL«ITWPBL*IX 

Y-XMIN/(IUNPBL*Y1) 
lY-Y 

IP  (Y.EC  lY)  GO  TO  100 
IF  (Y.LT.O.)  IY-IY-1 
100  IXL0W-IY*1UNPBL 

IXHI-IXLO; 

IF  (PXMAX.L."  IXHI)  GO  TO  110 

IX-IX+1 

GO  TO  90 

110  IL-IDNPBL*INT ( (PL0AT( IXHI) -FXMAX) / (2  *IDNPBt) ) 

IXLOW-IXLOW-IL 
IXHI  -IXHI  -IL 

120  IF  (10*(IXLOW/I0) .NS.IXLOW.OR.10*(IXHI/10) .NE.IXRI)  GO  TO  130 

IXLOH-IXLOW/10 


APPENDIX  F: 


CALCOMP  Generated  Plot  (Semilogarithmic)  of  the  GRPHX.XXX  Data 
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K£bPlKAIIIK  OUAHIIIATIVE  F11  lEbllllb 

ItlblKULllUtlb  li)l<  llblNIi  Hit  COMl'UIhK  iLKMlUALb  li;  UbAFbArt/VN  TD  FtAlFbb  lllL  IIAIA  COLLELTEn  Oi  THE  SALT  FOG 
llJblKUllLlJl 

KhV;  -  COMPUTER  GENERATED  INFORMATION  (CRT  SCREEN) 

-  Piogiam  Uie>t  Geni>iated  Inf^oimation  [t.nti'ied  via  the  kegboaid;  in^o^imation  dtiptaijed  o«  the  CRT 
6c>ieeyil 

-  Sequential  Step  Numbers;  not  displayed  on  CRT  screen 


—  Coininents  to  help  the  program  user 

LKI  SiKEEN  UlSHLAY 

>  _ 

2)  >  HULO  _ 

3)  ACCOUNT  OR  NAME:  _ 

<)  ACCOUNT  OR  NAME:  XXXXXX_ 

5)  PASSWORD:  _ 

6)  PASSWORD:  XXXX  _ 

.  (SERIES  OF  COMPUTER  SYSTEM  MESSAGES) 

ENTER  TERMINAL  TRACKING  NUMBER  (WHITE  (AG  ON 
RIGHT  FRONT)?  „ 


COniiENT 

CRT  display  status  normally  found  on  an  Idle 
terminal. 

Type  In  'HKl.LO';  depress  "Return"  key  on  keyboard. 
Computer  response. 

Type  In  your  last  name;  depress  "Return"  key  on 
keyboard. 

Computer  rusponse. 

Type  In  your  password;  password  characters  are  not 
displnynd  on  CRT  screen;  depress  “Return"  key  on 
kuybourd. 

No  response  required  on  your  part. 

laist  messUKe  ul'  the  Btuup. 
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LKl  bLKttN  UlbPLMY 


7)  tilTEH  terminal  tUIMBEH  (WHITE  TAG  0)  RIGHT 
FRONT)'!'  _ 


(SERIES  OF  COMPUTER  SYSTEM  MESSAGES) 


>_ 

8)  >  Run  Olt3i[30S.4.mCLKq,FT 


(PROGRAM  STATEMENTS) 


CU'liltlll 

tutor  two  dlijLts;  ilciiross  "Return"  key  on  keyboard. 
No  response  required  on  your  part. 

Computer  is  ready. 

Type  in  "Run  liKdi  I 'Jl)S,4  INACLIUJFT;  depress  tlie 
"Return"  key  on  tiio  keyboard. 

Ho  response  required  on  your  part 


9)  ENTER  THE  FOLLOWING!  (KJl  FOR  THE  FIRST  DATA  Last  of  proijraiii  statements. 

SET;  UU)!  FOR  THE  SECOND  DATA  SET;  (XJ3  FOR  THE 

THIRD  DATA  SET,  ETC* 

10)  entry  =  _  Enter  a  three  di^il  nur.iber  to  name  this  set  of 

.  "  data.  Depress  tlie  "Return"  key  on  keyboard. 


11) 


(PROGRAM  statements) 


No  response  required  on  your  (lart. 


12)  enter  THE  FOLLOWING!  UUl  FOR  THE  FIRST  RESIDUAL  “1  pro(jram  sl.ilemcnts. 

set;  (X)/I  for  the  SECOND  RESIDUAL  SET;  UU5  FOR 
THE  T.HIRD,  ETC. 

12)  entry  =  XXX  _ 


Enter  a  three  Jl,'. U  mmhur  to  name  this  set  of 
residuals;  depress  the  "Return"  key  on  the 
keyboard. 
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LKl  DlbPLAY  LUHtlU 

14)  . 

.  (PHOGRAM  statements)  response  reipilreR  on  your  part. 


15)  enter  the  following  UQl  FOR  THE  FIRST  GRAPH  Last  of  projiram  statements. 

SET;  UU2  FOR  THE  SECOND  GRAPH  SET;  003  FOR  THE 

THIRD.  ETC- 

16)  tNTRY  =  _  Enter  a  throe  iltglt  number  to  name  this  graph  set; 

depress  the  "Return"  key  on  the  keyboard. 

ll'I  ENTER  THE  NUMBER  OF  SODIUM  CHLORIDE  CALIBRATION  There  are  seven  sodium  chloride  calibration  con- 
CONCEHTRATION  STANDARDS-  centratlon  standards.  Enter  the  number  '7'; 

depress  the  "Return"  key  on  the  keyboard. 

ENTRY  =  ^  _ 

18)  DEPRESS  RETURN  KEY  AFTER  A  VOLTAGE  MEASUREMENT  No  response  required  on  your  part. 

19)  ENTER  THE  DATA  POINTS 


SAMPLE  NUMBER 

1 

2 

5 

‘I 

3 

D 

/ 


voltage 

MEASUREMENT 
(X  t»TA) 

_XXX_ 

XXX 

XXX 

“XXX'' 

XXX 

“XXX" 

XXX~ 


SODIUM  CHLORIDE 
CONCENTRATION 
(y  data) 

1.0 

~o.r 

~0.0T 

"O.OOT 

~O.OOOT 

~O.OOOOT 

‘O.QOOOOT 


Enter  a  voltage  measurement,  then  depress  the 
“Return"  key  on  the  keyboard.  Type  the  correspond¬ 
ing  sodium  ehlnriile  eallbratlon  concentration  then 
depress  the  "Return"  Key  on  the  keyboard.  Repeat 
this  procedure  until  nil  data  points  have  been 
entered. 
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l.KI  UI^IIAY 

20)  SUBJECT  UAtlE;  XXXXX_ 

21)  TYPE  CP  MASK:  XXXXX_ 

22)  DATE  TESTED;  _XXXXX_ 

23)  TIME  TESTED;  _XXXXX_ 

24)  THE  USER  IS  FREE  TO  SELECT  ONE  OF  TWO  GROUPS 
OF  EXERCISE  PROTOCOLS. 

•  (PROGRAM  STATEMENTS) 

25)  TO  SPECIFY,  THE  EXERCISE  PROTOCOL  (IROUP  OF 

Interest,  type  either;  Group  1  or  (jroup  1 

ENTRY  =  0<Oup  X  _ 

26)  . 


LUItItNI 


Enter  siiBieel’s  n.iiiii-  (up  to  30  characters  long); 
depress  the  '‘KetFirn"  key  on  tlie  keyboard. 

Enter  mask  noiiieiii:  l.iture  (up  to  30  characters  long), 
depress  "l{eturn"  key  on  keyboard. 

Enter  date  subjei'l  was  tested;  depress  the  "Return” 
key  on  tlie  keyhoai'il. 

Enter  tlie  time  ol  day  subject  was  tested  (for 
Example;  1430  hours);  depress  "Return"  key  on 
keyboard. 

First  of  a  series  of  program  statements  explaining 
the  two  groups  of  exercises.  No  response  required 
on  your  part. 


Last  of  program  statements. 

Enter  either  ’(.'roup  I’  or  'Croup  2'';  depress  the 
"Return"  key  on  the  keyboard. 


The  program  statements  explain  how  to  enter 
Integrator  data.  No  response  on  your  part. 


(PROGRAM  STATEMENTS) 
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I.Kl  :)«.KLLI1  lUotlAY 


LU'VltIJl 


2!)  EXEkCI5fc  cotjlir  UATA: 

integratok  time  period 

EXERCISE  COUNT  (IN  SECONDS) 

NORMAL  UREATHING 

STRAIGHT  AHEAD  XXXXXX  _XX_ 


Earner  liitL'^r.itiir  Ise  count  liuta  as  previously 

spectlieil  In  llio  program  scateiaencs;  depress  Che 
"KeCurn"  key  on  the  keyboard  after  each  Integrator 
count  value.  Eaiter  the  time  period  for  the  exer¬ 
cise  using  no  decimal  points;  depress  Che  "Keturn~ 
key  on  Che  keyboard.  Repeat  until  data  for  each 
exercise  l>as  been  input. 


(PROGRAM  statements) 


The  prograiii  slatei.ienCs  explain  the  functional 
definitions  of  the  variable  (X);  no  response  is 
required  on  your  piirt. 


29)  TO  SPECIFY  A  FUNCTIONAL  UEFINiriON  OF  (X) 
SELECT  THE  CORRESPONDING  NUMBER  INSIDE  DIE 
BRACKETS 

ENTRY  =  3 


30)  ENTER  THE  ORDER  OF  THE  DESIRED  POLYNOMIAL 
(MAXIMUM  =  THE  NUMBER  OF  SODIUM  Ct«.0RIDE 
CALIBRATION  CONCENTRATION  STANDARDS  *1) 

ENTRY  =  6 


31) 


(program  STATEMENTS) 


Last  of  prograi.i  statei.ients. 


Enter  the  nni.iber  '3';  depress  the  "Return”  key  on 
Che  keyboard. 

Enter  the  numhor  'b';  depress  Che  "Return”  key  on 
tile  keyboard. 


A  list  of  coef  f  ii:  lent  numbers  Is  printed  on  Che  CRT 
screen;  no  response  required  on  your  part. 


32)  Do  YOU  WANT  A  LIST  OF  RESIDUALS? 
(answer  Yes  ok  no)  Ye-i 


Enter  'Yes’,  unless  otherwise  instructed;  depress 
the  "Return"  key  on  the  keyboard. 
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CUlltlJl 

A  list  ut  residuals  Is  displayed  on  the  CKT  screen. 
No  response  rei|ulred  on  your  part.  If  you  want  to 
stop  the  flow  of  InforMUtion  on  the  CKT  screen, 
depress  simultaneously  the  ’CTRL*  and  'S'  keys  on 
the  keyboard.  To  resume  tlie  flow  of  Information  on 
the  CKT  screen,  depress  simultaneously  the  'CKTL' 
and  'Q'  keys  on  tlie  keyboard. 

The  descriptive  and  protection  factor  calculations 
are  displayed  on  CKT  screen;  no  response  required 
on  your  part. 


35)  Uo  YOU  WISH  TO  USE  THE  SAME  SOOItll  CHLORIDE  Type  In  'Yes'  or  'No'  dependtns  upon  your  desires; 

STANDARD  CALIBRATION  VOLTAGE  MEASUREMENTS  AND  depress  the  Return  key  on  the  keyboard. 

EXERCISE  INTEGRATOR  COUNT  DATA,  BUT  CALCULATE  A  -  /'***  fesult  will  be;  36) 

DIFFERENT  DEGREE  OR  FORM  OF  THE  POLYNOMIAL  CURVE  “  ^es  the  result  will  be:  29) 

FITTING  FUNCTION?  (ANSWER  YES  OR  liO) 

36)  Uo  YOU  WISH  TO  CALCULATE  PROTECTION  FACTORS  FOR 
ANY  OR  ALL  OF  THE  FOLLOWING  CONDITIONS: 
i)  A  DIFFERENT  SUBJECT  2)  A  DIFFERENT  SET  OF 
SODIUM  CHLORIDE  CALIBRATION  STANDARD  VOLTAGE 
MEASUREMENTS  t>)  A  01 FFERENT '  SET  OF  EXERCISE 
INTEGRATOR  COIKIT  DATA?  (ANSWER  YES  OR  IJo) 

37)  >  PRINT  DATA. XXX  If  you  want  a  cenipucdr  print  out  on  paper  of  the 

—  ~  data  that  you  Input,  type  'PRINT  DATA. XXX',  where 

XXX  Is  the  number  you  entered  from  the  keyboard  In 
step  9);  depress  the  "Return"  key  on  keyboard. 


Type  In:  'Yes’  or  'No'  depending  upon  your 
desires;  depress  the  "Return"  key  on  the  keyboard. 
_  If  you  tyi>u  In  'No'  Che  result  Is  37) 

—  If  you  type  In  'Yes'  the  result  la  9) 
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LKl  bLKhU'J  IJISHLWY 


LOf'.iLIJl 


38)  >  PKIWr  CALCX.XKX 


39)  >  PRIWr  GUPHX.XXX 


40)  >  eve 


(f  yon  wont  i  loi ipni i>r  print  out  on  {vipor  of  tlic 
result!!,  type  CAI.C  X.  XXX ' ,  where  tlie  first  X 

of  CAl.ttX.XXX  is  ;isslj;ned  the  niimhcr  I  hy  tile  con- 
puter  ,iiMi  1  111- rei.ientisl  by  I  cnch  time  .1  form  of  the 
polynomial  euree  littint’  fiinetlon  is  ralciilated, 
and  XXX  is  t  lie  noinher  you  entered  from  the  keyboard 
In  step  1 3);  depress  the  ■’Return"  key  oo  the  key¬ 
board. 

If  you  want  a  eonputer  print  out  on  paper  of  tlie 
graphing  data,  type  'PRINT  (StPIlX.  XXX '  ,  tdiere  tlie 
first  X  of  CHi’IIX,  XXX  Is  assigned  tlie  number  I  by 
the  computer  and  Incremented  by  I  each  time  a  form 
of  tlie  polynomial  curve  fitting  function  is  calcu¬ 
lated,  and  XXX  is  the  number  you  entered  from  the 
keyboard  in  step  Hi);  depress  tlie  "Return"  l:cy  on 
the  keyboard. 

If  you  are  finlshi>d  with  the  computer,  type  In: 
'BYE*  from  tlie  keyboard;  depress  tlie  "Return"  key 
on  keyboard  and  you  will  be  automat  teal  1  y  lugged 
off  tile  computer  terminal.  Computer  printed 
results  can  be  picked  up  In  computer  center,  Bldg 
150,  approximately  2  hours  after  you  log  off  the 
computer  terminal. 


IP 
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2-11  HLUIfllKi  KIIUIINL 

llJi>ll(llLlll)llb  H.K  UbIIJI.  THt  anPlllLk  iLmlll/lS  III  Ui)AFS/n/VIJ  TO  HJ)I  A  (iKAHII  IN  UAJ  UlMEHSKlUS. 

KLY:  -  COMPUTER  generated  INFORMATION  (CRT  SCREEN) 

-  PTOfliam  (Uei  Gtntnated  Infotmation  {f»te.ne.d  via-  the  kenhnnid;  iaio>imation  di-ipfaned  on  the  CPT 


icieeni 

-  Sequential  Step  Numbers;  not  displayed  on 

_  Cunnents  to  help  the  pro(;r;ici  user 

LHl  SLKLhN  DlbUAY  . 

n  >  _ 

^)  >  HELLO  _ 

3)  account  or  NAME:  _ 

^1)  ACCOUNT  OR  NAME:  XXXXXX_ 

5)  PASSWORD: 

6)  PASSV/OHD:  XXXX  _ 

« 

.  (SERIES  OF  COMPUTER  SYSTEM  MESSAGES) 


CRT  screen 


COTILIIl 

CRT  display  status  normally  found  on  an  Idle 
terminal. 

Type  In  'IIELLO';  depress  “Return”  key  on  keyboard. 
Computer  response. 

Type  In  your  last  name;  depress  "Return”  key  on 
keyboard. 

Computer  response. 

Type  In  your  passwurd;  password  characters  are  not 
displayed  on  CRT  screen;  depress  “Return”  key  on 
keyboard. 

No  response  required  on  your  part. 


ENTER  TERMINAL  TRACKING  NUMBER  (WHITE  TAG  ON  RIGHT  FRONT)?  _  l-ast  iiessnc.c  of  the  group. 
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Ikl  i>CKl-Ltl  DISPLAY  ClWILIJl 


7)  LNILK  TERMINAL  TRACKING  NUMBER  (WHITE 
TAG  ON  RIGHT  FRONT)?  XX_ 


(series  of  computer  system  messages) 


8)  Run  VR3i[30S.4MCLGRAPH 


9)  Mil  READ  RESPONSES  FROM  GRAPmI  .COM'aND  DO 
NOT  UPDATE  GRAPh1.COM. 

121  READ  RESPONSES  FROM  KEYBOARD  AND  UPDATE 
graph1.com. 

1 i I  HEAD  RESPONSES  FROM  KEYBOARD  AND  DO  NOT 
UPDATE  GRAPh1.COM. 

(iPllUN?  !<_ 


10)  OUTPUT  file  IfilespecIV  xxxxxx.xxx 


U)  GHAPH  HUMBER  N 

X-AXIS  MI Ml  mum/max I  MUM  OVEHRIDE  OPTION 

lY/HlV  V 


Enter  two  dl^itsi;  licpruKS  "Return"  key  on  keyboard. 


No  response  required  on  your  part. 


Computer  is  ready. 

Type  In  "Run  Dll'):  IdOS.AlNACLGRAPli;  depress  the 
“Return"  key  on  keyboard 

Computer  response  explaining  the  three  graphing 
options  of  the  program.  No  response  required  on 
your  part. 


Enter  the  number  '2'  or'T'i  depress  the 
“Return”  key  on  the  keyboard.  If  you  enter 
the  number  or  the  result  la  10).  If  you 

enter  the  number  the  result  Is  SB). 

Enter  the  name  u|  the  output  file  In  the  form 
'XXXXXX.XXX' ;  the  file  mime  may  have  up  to  nine 
eharaclers  before  the  period  ana  must  have  three 
tharacters  alter  the  period.  For  example, 
'PLOT. Dill '  i  depress  the  "Return”  key  nn  the 
keyboard. 

Enter  the  letter  'Y';  depress  the  “Return"  key  on 
tlie  keyboard. 
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LKl  M.KllN  lll'jIlAV 

12)  MINIMUM  VALUEV 

13)  MAXIMUM  VALUE? 

14)  Y-AXIS  MINIMUM/MAXIMIM  OVERRIDE  OPTION 

lY/NI? 

15)  MiHlMUH  VALUE?  XXX_ 

161  MAXIMUM  VALUE?  XXX_ 

171  PLOT  SPECIFICATIONS: 

111  LINEAR  |?1  SEMI-LOG  lil  LOG-LOG 
TYPE  OF  SCALE? 

18)  number  of  cycles  1>0|?  7_ 

19)  RETRACE  OPTION  lY/Nl?  Y_ 

201  GRAPH  NUMBER  N:  X’AXIS  SPECIF ICATtONS: 
X-AXIS  LENGTH  IN  INCHES?  XX_ 

211  DISTANCE  BETWEEN  SLIPS  IN  INCHES?  XX_ 


CUItltNl 


Enter  the  smallest  value  on  the  X-Axis;  depress  the 
"Return"  key  un  the  keyboard. 

Enter  the  largest  value  on  the  X-axls.  Depress  the 
“Return"  key  on  the  keyboard. 

Enter  the  letter  'Y';  depress  the  "Return"  key  on 

the  keyboard. 

Enter  the  smallest  value  on  the  Y-axis;  depress  the 
"Return"  key  on  the  keyboard. 

Enter  the  largest  value  on  the  Y-axis;  depress  the 
"Return"  key  on  the  keyboard. 

Enter  the  number  '2'  unless  otherwise  instructed; 

depress  the  "Return"  key  on  the  keyboard. 

Enter  the  number  unless  otherwise  instructed; 

depress  the  "Return"  key  on  tlie  keyboard. 

Enter  the  letter  'Y';  depress  the  "Return"  key  on 

the  keyboard. 

Enter  the  length  of  the  X-axls  in  inches 
number  not  m  word;  depress  the  "Re' 
keyboard. 

Enter  the  dlstanve  In  in'' 

X-axis.  Enter  .t  no- 
"Return"  key  un  tlie 
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LKi  y.KlLli  DibliAY 


SCALE  UPriON  IV/Ul?  V. 


^J)  niMPH  tlllfIBCK  IJ;  Y-AXIS  SPECIFICAf  IlMIS: 
Y-AXIS  LENC.TH  in  inches?  XX  _ 


<?4)  DISTANCE  HETWEEN  BLIPS  IN  INCHES?  XK_ 


<^b)  SCALE  OPTION  lY/IJl?  W 


26)  ORAPH  NllfIBER  IJ:  X  AND  Y  AXIS  SPECIFICATIONS; 
ID  BLACK  I?  I  RED  I  31  C.REEN 
Color?  L 

^7)  CHARACTER  SIZE  IN  INCHES  1 -U/  *  -bl?  •XX_ 

^H)  r.RAPM  NUMBER  IJ  TITLE  SPECIFICATIONS; 

X-AXIS  TITLE  ll-‘IU)?  XXXXX_ 

?y)  CHARACTER  SIZE  IN  INCHES  l.U/-.bl?  XX_ 

30)  ID  liLACK  I?  I  KEI)  131  (iREEN 
(.OLOR?  L 


Dtlll.iJl 


Cuter  the  letler  'N'  unless  utherwlse  Instructed; 
depress  till'  "Ueturn"  key  on  the  keyboard. 

Enter  the  lenj;tli  ul  the  Y-axls  In  inches.  hunter  a 
number  not  a  wui'il;  depress  the  "Return"  key  on 
keyboard. 

Enter  tiu-  dlsl.iiui-  in  inches  between  points  on  the 
Y-axIs.  Enter  .i  number  not  a  word;  depress  Che 
"Return"  key  i>n  keyboard. 

Enter  tlio  lettei  ’I."  unless  otherwise  Instructed; 
Repress  the  "Return  key  on  the  keyboard. 

Enter  the  niiiiber  Depress  the  "Return"  key  on 

tile  keyboard. 

Enter  the  size'  ol  the  number  to  be  displayed  along 
the  X  .ind  Y  axis.  Enter  a  number  between  0.07  and 
0.5  In  Inches.  Depress  the  "Return"  key  on  the 
keyboard. 

Type  In  the  title  lor  the  x-axis.  Limit  the  number 
of  charaeters  to  AO;  depress  the  "Return"  key  on 
the  keyboard. 

Enter  the  size  of  the  characters  for  the  X-axls 
title  In  Inches.  Type  a  number  between  0,07  and 
0,5;  depress  the  "Return"  key  on  the  keyboard. 

Enter  tlie  number  depress  the  "Return"  key  on 

the  keyboard. 
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LiU  ^LKLLU  IJlbUrtY 


■it)  Y-AXIS  TITLE  I  l-itU  CHARS- IV  XXXXX_ 

■tV)  CHARACTER  SIZE  IN  INCHES  l-O/'-blV  XX_ 

33)  111  Black  IVl  kEn.  131  Green 

Color?  t_ 

34)  GRAPH  TITLE  (LINE  1)  ll-bO  CHARCS- |V  XXXXX 

3!>)  CHARACTER  SIZE  IN  INCHES  ) -.U/' -5 IV  XX_ 

36)  111  Black  I2l  Ked  13)  (jreen 

LaoR?  t_ 

37)  ck«PH  TITLE  (line  2)  ll-bU  CHARS- IV  XXXXX_ 


3«)  CHARACTER  SIZE  IN  INCHES  l-O/'-blV  XX 


ClYVTLIil 


Type  "In  tin-  1 1  lo  for  tlie  Y-axIs.  Limit  the  number 
of  chiir.ieters  ti>  4(1;  depress  the  "Return"  key  on 
the  keyboard. 

Enter  the  sl/e  ol  the  characters  for  the  Y-axis 
title  in  inches.  I'.iiter  a  number  between  0.07  and 
0.5;  depress  tlie  “Return"  key  on  the  keyboard. 

Enter  the  niinber  '1';  depress  the  "Return"  key  on 
the  keyboard. 

Type  in  the  first  line  of  the  graph  title.  Limit 
the  number  of  characters  to  60;  depress  the 
"Return"  key  on  the  keyboard. 

Enter  the  size  "of  the  characters  in  the  first  line 
of  the  gr.aph  title  in  Inches.  Enter  a  number 

between  0.07  and  0.5;  depress  the  "Return"  key  on 
the  keyboard. 

Enter  the  number  depress  the  "Return"  key  on 

the  keyboard. 

If  a  line  of  the  graph  title  Is  to  be  centered 

under  a  preceding  line.  Include  an  appropriate 
number  of  blank  spaces  before  typing  the  title. 
Type  in  'the  second  line  of  the  grapli  title.  Limit 
the  number  of  characters  to  60;  depress  the 
"Return”  key  on  the  keyboard. 

Enter  the  size  of  the  characters  In  the  second  line 

of  the  graph  title  in  inches.  Enter  a  number 

between  0.07  and  0.5;  depress  the  “Return"  key  on 
the  keyboard. 


cf 


Cf} 

CO 
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CKI  SLkLtll  DlbHLAY  CUVILIll 


39)  11)  Hack  l)!l  Ked  131  bREEH 
Color?  L 

'10)  GRAPH  TITLE  (LINE  i)  11-bU  CHARS.  I?  XXXXX_ 


41)  CHARACTER  SIZE  IN  INCHES  l.U/-.bl?  XX_ 


42)  111  Black  121  Ked  13)  (ireen 
Color?  L 

4'3)  graph  title  (line  9)  ll-bO  chars)?  xxxxx_ 


44)  CHARACTER  SIZE  IN  INCHES  1.0/-.bl?  XX_ 


4b)  111  Black  |21  Bed  |31  (ireeii 

COLOR?  L 

46)  GRAPH  NUMBER  II 

INPUT  FILE  FOR  GRAPH  IJ  I  F I LE:;i’!  (; |? 

XXXXXX.XXX 


Enter  the  niinher  depress  the  "Return"  key  nn 

the  keyboard. 

Type  In  the  tliird  lino  of  the  (traph  title.  Limit 
the  number  of  eluiracters  to  60;  depress  the 

"Return"  key  on  the  keyboard. 

Enter  the  size  of  the  characters  In  the  third  line 
of  the  yraph  title  In  Inches.  Enter  a  number 
between  0.0?  and  li.S;  depress  the  "Return"  key  on 
the  keyboard. 

Enter  tbe  number  depress  the  "Return”  key  on 

the  keyboard. 

Type  In  the  fourth  line  of  the  graph  title.  Limit 
the  number  of  characters  to  60;  depress  the 

"Return"  key  on  the  keyboard. 

Enter  the  size  of  the  characters  In  the  fourth 
line  of  the  graph  title  In  Inches.  Enter  a  number 
between  0.07  and  0.5;  depress  the  "Return"  key  on 
the  keyboard. 

Enter  the  number  '1';  depress  the  "Return"  koy  on 

the  keybo.ird. 

Enter  the  name  of  the  li.put  file  In  the  form 

'XXXXXX.XXX';  the  file  name  may  have  up  to  nine 
characters  before  the  period  and  must  have  three 

characters  after  the  period.  for  example, 
'CRI'lll.OOl  ' ;  depress  the  "Return"  key  on  the 
keyboard. 
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CKl  SlKtlM  DbKAY 


'»7)  'actual  nATA*  curve. 

Ill  Black  I'/?!  Ben  151  liREEtL 

(.OLOR?  L 

48)  ‘ESTlt\AIED  values'  CURVE. 

Ill  Black  I'^l  Bed  151  Breeh 
ColorV  I_ 

49)  CREATE  ANOTHER  GRAPH  lY/Ul? 


90)  SUCCESSFUL  CmtFLETlON 


91) 


52)  )«k)t>«SO.CMP_ 


53)  >*I11SK  DRX  Isl:  X. 


54)  >*UIC  m,n  (si:  m.XK^ 


55)  >*FILeNAME.EXT  Isll  XXXJ(XX.XXX_ 


cmiLni 


Enter  the  number  'I';  deproett  the  "Return"  Key  on 
the  keybourd. 


Enter  the  number  *1';  .iepress  the  "Return"  key  on 
the  keyboard. 


Enter  either  the  le' ter  'V*  or  'N*;  depress  the 
“Return"  key  on  the  keyboerd. 

Coaputer  response.  U  no  other  graph  is  to  be 
crested.  No  response  required  on  your  part.  If 
another  graph  Is  to  be  treated,  the  coaputer  util 
return  to  step  •>(  the  directions. 

Coaputer  is  reuly  to  plot  the  graph  on  the  Calcoap 
plotter. 

Enter  •^tPHSU.CttO* ;  depress  tl»e  "Return"  key  un  the 
keyboard. 

Enter  th"  disk  uiw'oer  of  the  account  that  you 
Skgn.'d  on  the  torulnal  ultht  depress  the  "Return* 
hey  on  ttu;  kvyha.trd. 

E.  ir  the  UlC  nuuhur  of  tl»e  account  that  you  signed 
on  the  ternlt.al  with;  use  the  fora  XXX, XXi  depress 
the  "'•etuni"  kuy  m\  the  keyboard. 

Enter  the  mine  uf  llie  output  file  tor  the  CRAEH. 
Use  the  torn  VXXKX.'CX. XXX'.  This  util  be  the  saae 
naae  as  In  step  lu)  nt  the  directions. 
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LKI  SLKtt.ti  IMSUrtY 

^>b)  *  Plain/lineo  Is  I:  Plain 

•  (SEHIES  OF  COMPUTER  STAFEMEMts) 

sa)  >8«e_ 

">'*)  CWIIPIIT  FtVF  iFlLlSPedV 

«0>  fiRAPti  NllMtttR  H 

INPlU  FILE  for  rOJAPSi  N  iFaESPtcIV  WWf.JtKX, 


LiniUii 


Kiitur  till'  w.)rit  'I’l.iln*  tiiilL-sii  otherwise  liistrueced; 
iiei>ress  the  "Ketiirii"  key  on  the  keyboardt 

No  resi)onse  iviiolreil  im  your  |>>irt. 


host  ot  lOMiputer  'it .tteuents;  OKAPtI  has  been  created 
and  stored. 

II  you  are  tlol-ilsed  uslnit  the  coaputer,  type  tn 
•Hye*  t  rom  the  keyhi>ard;  depress  the  “Return"  key 
on  tile  keyho.iril  and  you  will  he  autoaatically 
loitlied  ott  ils'  loupaier  teralnaU 

The  OkiVHII  eaa  Iw  pleked  up  in  the  coaputer  center, 
Bldli  ISO,  tn  approelDMIoly  1  hours. 

tHrectluns  tor  uptlun  I 

Knier  the  naiae  ni  the  output  lUe  tn  the  fora 

‘XXy.XXX.XXS* ;  tlw  11  le  naae  aay  have  up  to  nine 

charaetefs  tnuore  the  period  sad  «wst  liave  tliree 

eharaeleu  aiiet  the  period.  For  enaopte, 
•uRAPill.ooi ' ;  .leprestt  the  “return"  key  on  itu? 
keyhitard. 

Knter  the  naae  of  the  Input  ttle  tn  the  loro 

'X.XXXXX.XXX’ j  the  Itte  oaae  »My  have  up  to  nine 

ehararters  Indove  the  period  and  must  twve  three 

eharaelefs  alter  Ih..  perlodt  depress  the  “keturn" 
key  Ott  tile  keyboard. 
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LKl  bLItttJ  UlbPLAY 


LUniHJl 


6U  CHEATE  ANOTHEH  GHAPH  IY/(J1V 


6tl  SUCCESSFUL  COMPLETION 


E^ntor  the  lotter  'Y*  or  'N';  depress  the  "Return” 
key  on  the  keyboard.  If  you  enter  the  letter  'Y* 
the  computer  goes  to  direction  9);  If  you  enter  the 
letter  'tv*  the  coinputer  (toes  to  62). 

Computer  response;  no  response  required  on  your 
part.  Computer  now  ruturi\s  to  direction  $1). 


ABBREVIATIONS,  ACRONYMS.  AND  SYMBOLS 


ASCC 

C 

Ca 

COE 

cm 

CU 

OB 

dc 

OTL 

Ein 

Eos 

Eout 

e* 

re 

» 

‘s 

1C 

(.d. 

ITT 

K  or  k 
kP« 

UO 

vA 

uE 

urn 

ws 

i>V 

lA 

•In 


Air  Standardization  Coordinating  Committee 
capacitor 

ambient  challenge  atmosphere  concentration 
Chemical  Defense  Establishment  (United  Kingdom) 
centimeter 

sampled  leakage  concentration 

chemical  Marfare' 

deep  breathing 

direct  current 

diode-transistor  logic 

input  voltage 

offset  voltage 

output  voltage 

exponential  of  t 

facial  grimacing 

current;  or,  the  i^^  exercise 

input  bias  current 

Integrator  count 

inside  diameter 

current  leakage  through  the  integrator  capacitor 

tniemational  Telephone  atsd  Telegraph 

kilo  (1000) 

kilooPascals 

light  emitting  diode 

logarithm  to  the  base  ten 

microampere 

microfarad 

micrometer 

microsecond 

microvolt 

mlltiampere 

minute 

(Coat'd,  on  facing  page) 
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